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Abstract—Movements of the tiger shark 
(Galeocerdo cuvier) by life stage are 
largely unknown yet are necessary 
to determine essential fish habitat 
and sustainable fishery management 
practices. In an effort to elucidate dis- 
tribution and movements of the tiger 
shark, we analyzed fishery-dependent 
and -independent tag (sample size 
[n]=10,516) and recapture (n=762) 
records for sharks caught in the North 
Atlantic Ocean during 1963-2018. Sea- 
sonal distribution of the tiger shark was 
examined by life stage—young of the 
year, juveniles (immature), and adults 
(mature)—and distribution patterns 
were used to identify potential nursery 
habitat. Tiger sharks were caught over 
a wide area from the Grand Banks of 
Newfoundland, Canada, south to Bra- 
zil and from coastal to offshore waters 
and into the eastern North Atlantic 
Ocean. Seasonal north-south move- 
ments were observed in all life stages, 
and 14 immature sharks were found to 
have migrated from the western to the 
eastern North Atlantic Ocean. A broad 
nursery area and a potential birthing 
area were identified on the continental 
shelf between Florida and Georgia on 
the basis of the repeated presence of 
neonates in summer across years and of 
the recapture of multiple tagged young- 
of-the-year sharks from the same loca- 
tion over a period of at least 2 years. 
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The tiger shark (Galeocerdo cuvier) is a 
wide-ranging, globally distributed pred- 
ator found in tropical and warm tem- 
perate waters (Bigelow and Schroeder, 
1948). Tiger sharks make extensive 
migrations in many parts of the world, 
and the range of this species is consid- 
ered both coastal and oceanic (Heithaus 
et al., 2007; Ferreira et al., 2015). In the 
western North Atlantic Ocean (WNA), 
the tiger shark is known to undertake 
seasonal inshore—offshore movements 
in the Gulf of Mexico (GOM) and sea- 
sonal north-south movements along 
the East Coast of the United States 
(Lea et al., 2015; Ajemian et al., 2020). 
Evidence also exists for basin-wide 
migrations in the North Atlantic Ocean 
(Kohler et al., 1998; Kohler and Turner, 
2019). Tiger sharks have been found in 
water temperatures ranging from 4.0°C 
to 34°C (Afonso and Hazin, 2015; Ham- 
merschlag et al., 2022), but in the WNA 
they seem to prefer water temperatures 
between 26°C and 28°C, regardless of 
season (Hammerschlag et al., 2022). 
Because of its broad distribution, this 
species occurs in the waters of many 
countries and crosses numerous man- 
agement boundaries. Therefore, under- 
standing the movements of tiger sharks 
in general, and by life stage, is import- 
ant to international management. 


Tiger sharks are caught in commer- 
cial and recreational fisheries through- 
out the North Atlantic Ocean (Natanson 
et al., 1999; Domingo et al., 2016). The 
tiger shark is currently listed on the 
IUCN Red List of Threatened Species 
as near threatened globally (Ferreira 
and Simpfendorfer, 2019). However, in 
the North Atlantic Ocean, shark vulner- 
ability to pelagic longline gear has been 
estimated to be low to moderate on the 
basis of an ecological risk assessment 
based on the species’ productivity and 
susceptibility to the gear (Cortés et al., 
2015). Furthermore, increases in tiger 
shark abundance in the WNA have 
been documented since the early 1990s, 
a trend attributable to the implementa- 
tion of the National Marine Fisheries 
Service (NMFS) Shark Management 
Plan in 1993 (Carlson et al., 2012; 
Peterson et al., 2017). 

Life history studies of tiger sharks 
have expanded in recent years. Results 
of the most recent study on growth 
dynamics, in which ages validated with 
bomb radiocarbon techniques were 
used, indicate a longevity of 27-29 
years and an age at maturity of 10 years 
for both sexes (Kneebone et al., 2008). 
Updated reproductive parameters cor- 
roborated these estimates, indicating 
a size and age at median maturity of 
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258.9 cm in fork length (FL) and 9.5 years for males and 
of 261.4 cm FL and 11.6 years for females (Natanson et al., 
2023). Tiger sharks are well known as generalist predators 
(Lowe et al., 1996), and they have a high trophic level indica- 
tive of a top predator (Cortés, 1999). They are considered an 
important component of the trophic structure for reef sys- 
tems and other coastal habitats (Heithaus, 2001; Ruppert 
et al., 2013), and results of studies indicate that ontogenetic 
shifts in their diet may partly result from sharks of different 
sizes or sharks in different life stages occupying different 
habitats (Lowe et al., 1996; Aines et al., 2018). 

Movement and migration information on the tiger shark 
is limited to discrete studies in which high-technology tag- 
ging methods were used with limited numbers of satellite 
or acoustic tags. The results of these studies indicate that 
seasonal migrations to higher latitudes are common and 
that, although there are resident individuals, migrations 
over long distances are also undertaken (Sulikowski et al., 
2016; McClain et al., 2022). In some studies, longer migra- 
tions have been attributed to males (Lea et al., 2015), but 
more recent data indicate that they are undertaken by both 
sexes (McClain et al., 2022) and are likely related to forag- 
ing, especially in relation to turtle nesting sites (Fitzpatrick 
et al., 2012), temperature (Papastamatiou et al., 2013; Lea 
et al., 2018; Hammerschlag et al., 2022), or important hab- 
itats related to life stage (e.g., birthing and nursery areas) 
(Meyer et al., 2018). Data from these types of studies, how- 
ever, are limited by the number of tagged individuals and a 
lack of representation by all life stages, and such data may 
not provide a complete picture of movements for the species 
(Meyer et al., 2018; Sequeira et al., 2019). 

In this study, we used updated reproductive parameters 
(Natanson et al., 2023) in conjunction with conventional 
tagging data to examine the demographic characteristics of 
the population of tiger sharks in the North Atlantic Ocean. 
Specifically, we looked at potential variation in spatial and 
seasonal distributions and migrations by sex and life stage, 
young of the year (YOY), juveniles (immature), and adults 
(mature). We also examined locations of mature females 
and neonate individuals (<66.2*' cm FL) from March 
through August by using estimated time and size of birth 
(Castro, 2011) to identify possible areas of parturition and 
by using the distribution of YOY to identify potential nurs- 
ery areas for the tiger shark. Information on movements by 
life stage can be used in management to understand essen- 
tial fish habitat and to help define protected areas or other 
approaches to protect vulnerable portions of the population. 


Materials and methods 
Mark and recapture 
Tag and recapture data from the NMFS Cooperative 


Shark Tagging Program (CSTP), for the period 1963-2018, 
were used with published maturity estimates (Natanson 


1 An asterisk (*) after a size value denotes a fork length that has 
been converted from a total length. 
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et al., 2023) to investigate distribution and movements 
of tiger sharks by sex and life stage. Tiger sharks were 
tagged by cooperating commercial and recreational fish- 
ermen, scientists, and fisheries observers. Tagging meth- 
ods and quality control have remained consistent since 
the inception of the CSTP and are detailed in Casey and 
Kohler (1992) and Kohler and Turner (2019). Distances 
between release and recapture locations are expressed 
in straight-line nautical miles (nmi). Tag and recapture 
data were reported in FL, total length (TL), or weight 
(WT) and as measured or estimated. Lengths are pre- 
sented in centimeters, and weights are provided in kilo- 
grams. All tagging and recapture data were converted to 
F' Lor, (i.e., from the tip of the snout to the fork in the tail, 
over the body [OTB]) by using the I-MARK web appli- 
cation (Kohler and Turner, 2019) and conversions from 
Kohler et al. (1996). Published straight-line TL (TLg7p) 
data were converted to FLo7pp (these converted values are 
indicated with asterisks) by using the conversion from 
Natanson et al. (2022): 


TLorp = 11.90 ae 1.180'Lorp), 
coefficient of determination [r7]=0.996, n=605. 


Unless otherwise noted, all lengths are presented in 


Demographic structure 


To investigate the spatial distribution of the population, 
tiger shark catch records of tagging and recapture events 
were divided into subsets by sex and life stage, where 
possible. Classification into stages was done by length. 
Sharks <106 cm FL were considered YOY (Natanson et al., 
1999), sharks >106 cm FL and less than the median size 
at maturity by sex (258.9 and 261.4 cm FL, for males and 
females, respectively; Natanson et al., 2023) were classi- 
fied as juveniles (immature), and sharks greater than or 
equal to the median size at maturity were classified as 
adults (mature). For sharks for which no sex was recorded, 
maturity status was determined by using the more con- 
servative, larger, female size at maturity. Additionally, 
locations at which sharks <66.2 cm FL (neonates) were 
caught from May through August were examined as possi- 
ble locations of birthing areas on the basis of an extension 
of the timing of parturition from Castro (2011) (June and 
July) to account for early and late parturition. The length 
at birth was based on the range (57.7*—66.2* cm FL) pro- 
vided by Castro (2011) with the assumption that any free- 
living shark <66.2 cm FL had just been born. A chi-square 
test was used to determine if there was a significant dif- 
ference between the sample sizes of females and males in 
the tagging data. 

To visualize potential variation in the spatial distribu- 
tion of the life stages of tiger sharks, catch locations by 
sex, life stage, and season were mapped by using R (vers. 
4.1.3; R Core Team, 2022). Seasons were defined as win- 
ter (December—February), spring (March—May), summer 
(June-August), and fall (September—November). The 
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distribution of catch records for neonates, adult females, 
and YOY were examined for patterns related to birthing 
and nursery areas. 


Results 
Tagging and recapture 


Tiger sharks (n=10,516) were tagged in coastal waters 
throughout the continental shelf in the WNA from the 
Grand Banks of Newfoundland, Canada, south to Brazil 
and in offshore waters (past the 200-m depth contour) to 
40°W (Table 1, Fig. 1). The majority of tags (n=9211) were 
deployed along the East Coast of the United States from 
Cape Hatteras, North Carolina, down and into the U.S. 
waters of the GOM. Sharks tagged were primarily cap- 
tured on longline gear (80.8%) during commercial (32.7%) 
or scientific (83.2%) fishing operations or by fisheries 
observers (17.0%), followed by recreational anglers using 
rod and reel (17.0%). All industries fished in the same 
general locations, although spatial effort varied between 
them (Table 1, Fig. 1). Recreational anglers, in general, 
fished closer to shore, and commercial and scientific 
operations fished both inshore and offshore, particularly 
between northeastern Florida and north of Cape Hatteras 
(Table 1). Commercial effort was primarily focused on the 
east coast of Florida, and the bulk of the scientific effort 
was evenly distributed along the coast from Cape Hatteras 
to the Florida Keys. Observer coverage was highest from 
Cape Hatteras to the border of Florida and Georgia, and 
although recreational fishing was highest north of Cape 
Hatteras, it was fairly evenly distributed south of the cape 


and into the GOM but very limited in the Caribbean Sea 
and the Bahamas (Table 1). 

Sharks were captured in all months of the year. Lengths 
of tagged tiger sharks encompassed the entire size range 
of the species (from birth to 391.0* cm FL; Bigelow and 
Schroeder, 1948), with a range of 40.1-457.2 cm FL for 
both sexes combined and a mean of 120.2 cm FL (stan- 
dard deviation [SD] 51.1) and 134.7 cm FL (SD 60.4) for 
males and females, respectively (Fig. 2A). The lengths of 
5 individuals were greater than the published maximum 
size, but because these lengths were estimated, none can 
be considered a definitive new maximum size; 3 of these 
sharks were females, and 2 of them were of unknown sex. 
Although caution should be exercised in taking these val- 
ues as a new maximum size for tiger sharks, the sharks 
would be classified as mature, and, for that reason, these 
lengths were included in our analyses. Of those tagged 
sharks for which sex was reported, the proportions of 
females (57.3%) and males (42.7%) were significantly dif- 
ferent from the expected ratio of 1:1 (X*(1, n=9419)=200.14, 
P<0.05). 

A total of 762 tiger shark recapture events were 
reported between 1967 and 2018 throughout the North 
Atlantic Ocean (Fig. 3), with an overall recapture rate of 
7.2%. Of these, 35 individuals were recaptured multiple 
times (3 sharks were recaptured twice, and 2 sharks 
were recaptured 3 times). For 13 recaptured sharks, no 
information was turned in at tagging and their tag data 
could not be used for movement analysis. In all but 
5 cases, sex was reported the same at both tagging and 
recapture. Sex was reported as unknown at recapture for 
2 of those exceptions, and sex was reported unknown at 
tagging for the other 3 cases. Similar to the proportions 


Table 1 


The numbers of tagged tiger sharks (Galeocerdo cuvier) from the western North Atlantic Ocean caught between 1963 and 2018 
by area tagged for each sex, life stage, and tagger industry. Also provided is the percentage of tagged fish that were mature for 
each sex. Tagged sharks were classified into 3 life stages: young of the year (YOY), juveniles (immature [IMM]), or adults (mature 
[MAT]). The tagger industries include commercial (COM) and scientific (SCI) fishing operations, fisheries observers (OBS), and 
recreational anglers (REC). The latitudinal borders of the tagging areas are as follows (in decimal degrees): north of Cape Hatteras 
(35.2500—45.5833°N), from Cape Hatteras to Florida (35.2499-30.7167°N), east coast of Florida (30.7168—25.3400°N), west coast of 
Florida (25.3401—30.2600°N), rest of the Gulf of Mexico (GOM, 302601—21.0333°N), Caribbean Sea (15.5667—21.3000°N), and the 
Bahamas and adjacent waters (13.0000—25.0833°N). The total number of sharks tagged includes those of unknown sex or maturity 
and 2 fish caught outside the location parameters. 

Males 


Females Tagger industry 


Total no. % % 


Area tagged YOY IMM MAT MAT YOY IMM MAT MAT COM SCI OBS REC 
North of Cape Hatteras 1318 141. 392 23 ai 81 209 22 7.0 160 164) = 2 EST 2 BHO 
Cape Hatteras to Florida 3066 Tal. 727 73 4.6 506 39 3.2 298 1424 1062 282 
Florida East Coast 4692 1253 1135 95 3.8 762 35 LL]? .,.2008. 1422.-. 273. ~ 185 
Florida West Coast 677 OA Bao 23 6.7 143 8 32 80 too” “270 at? 
Rest of GOM 419 ga 82 4.0 87 10 8.1 25 109 AQ? * 22 
Caribbean Sea 143 41 20.0 28 44 64 

Bahamas and adjacent waters 229 1B EZ) 19:3 36 29 194 


Total 10,516 2810 4.9 1771 , 3444 3486 


148 


Fishery Bulletin 121(4) 


Figure 1 


Map showing the distribution of all tiger sharks (Galeocerdo cuvier) tagged in the western North 
Atlantic Ocean by participants in the Cooperative Shark Tagging Program of the National Marine 
Fisheries Service between 1963 and 2018. Circles indicate locations where sharks were tagged. 
The dashed black line indicates the 200-m contour. GOM=Gulf of Mexico; CS=Caribbean Sea. 


of tagged sharks, recaptured females (56.2%) outnum- 
bered recaptured males (43.8%) for sharks for which sex 
was reported at recapture (n=714). Sizes of recaptured 
tiger sharks ranged from 66.0 to 373.0 cm FL with a 
mean of 142.8 cm FL (SD 50.4) and 154.7 cm FL (SD 
56.6) for males and females, respectively (Fig. 2B), nearly 
representing the entire size range of this species. Com- 
mercial fishermen (69.3%) recaptured the majority of 
tiger sharks, followed by recreational anglers (18.0%), 
biologists (5.8%), fisheries observers (4.5%), unknown 
(2.4%), and educators (<1%). The majority of sharks 
recaptured were caught on longline gear (75.7%) during 
commercial fishing operations. 

The range of overall time at liberty for tiger sharks was 
0.0-13.4 years with a mean of 1.0 year (SD 1.4). Most 
sharks were at liberty less than a year (71.8%). Time at lib- 
erty varied depending on life stage, with mature males hav- 
ing a higher mean time at liberty (5.7 years [SD 3.9]) than 
the means for males at all other maturity stages (Table 2). 
Distance between tagging and recapture locations ranged 
from 0 to 3643 nmi with a mean of 243.2 nmi (SD 449.5). 
Although mean distance between tagging and recapture 
locations was higher for mature sharks of both sexes, imma- 
ture fish had higher maximum distances (Table 2). 

Recapture locations were concentrated close to the coast 
between Long Island, New York, and Florida, into the U.S. 
portion of the GOM. Additional recapture events occurred 
along the coast of Mexico in the GOM into the Caribbean 
Sea and south to Brazil. North of Long Island, sharks 
were recaptured farther offshore. The northernmost recap- 
ture was on the Grand Banks of Newfoundland, and the 


easternmost recapture was in the coastal waters of Guinea- 
Bissau, a country in West Africa (Fig. 3). 

Tiger sharks transited throughout the WNA, including 
the GOM and Caribbean Sea (Fig. 3). A total of 37 sharks 
crossed between these water bodies. The majority of these 
individuals moved from the main body of the Atlantic Ocean 
into either the GOM (n=11) or the Caribbean Sea (n=14). 
Movement was also seen from the GOM into both the main 
body of the Atlantic Ocean (n=5) and the Caribbean Sea 
(n=5) and from the Caribbean Sea into both the main body 
of the Atlantic Ocean (n=1) and the GOM (n=1). 

A total of 116 sharks, of both sexes, traveled more than 
500 nmi; of these, 115 individuals had length data at tag- 
ging and 109 individuals had length data at both tagging 
and recapture. The majority of tiger sharks that travelled 
over 500 nmi were immature (82%, 84% including YOY), 
reflecting the prevalence of this life stage in the tagging 
data (Table 2). Even though adults account for only 5% of 
the total recapture events, this life stage had the greatest 
percentage of recaptured sharks that travelled over 500 
nmi. (Table 2). Of those sharks that travelled >500 nmi, 23 
individuals were tagged when YOY, 88 sharks were tagged 
when immature, and 4 fish were tagged when mature. Six 
fish were immature at tagging and mature at recapture. 
Two female tiger sharks recaptured in 1995 had traveled 
from the East Coast (off Frying Pan Shoals in North 
Carolina and off St. Augustine, Florida) across the middle 
of the Atlantic Ocean (45°W). One of these females, imma- 
ture at both tagging and recapture, traveled to 41.5°W and 
was recaptured after 571 d at liberty. The other female 
was a YOY at tagging and still immature when recaptured 
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Figure 2 


Size-frequency histograms for (A) tagged and (B) recaptured tiger sharks (Galeocerdo cuvier) 
from the North Atlantic Ocean for the period between 1963 and 2018. 


off the coast of Guinea-Bissau at 17.5°W; the transit of this 
shark remains the longest distance traveled to date for a 
tiger shark in the CSTP (8648 nmi; Kohler and Turner, 
2019) (Fig. 4). Since 1995, an additional 12 individuals 
that had been tagged in the coastal waters of the WNA 
have been reported to have been recaptured in the eastern 
North Atlantic (ENA). These sharks were all immature at 
tagging and recapture and consisted of 6 females, 5 males, 
and 1 individual of unknown sex (Fig. 4). 

All of the 35 tiger sharks that were recaptured multi- 
ple times were immature at tagging and final recapture 
(Table 3). Twenty-two (63%) of these sharks were origi- 
nally tagged as YOY. Multiple mark-recapture records 
indicate that individual tiger sharks tagged and recap- 
tured as YOY, and up to 124 cm FL, were recaptured in 


the vicinity of their tagging location for up to 3 months 
and/or returned to the area where they were tagged for 
up to 2 years (Table 3). The continental shelf off northern 
Florida was used by YOY each month of the year. Several 
of the immature sharks also made long-distance move- 
ments; for example, a 91-cm-FL male was tagged off 
northern Florida in November, recaptured there the fol- 
lowing May, and recaptured again in December nearly 
1400 nmi away in the Caribbean Sea off Venezuela 
(Table 3). A 134-cm-FL female was tagged and recap- 
tured off North Carolina in September of the same year 
and was recaptured again 2 years later east of the 
Mid-Atlantic Ridge (>2300 nmi away). Only 1 shark was 
recaptured multiple times in the GOM, a 156-cm-FL 
female at tagging that used the northeastern GOM over 
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Figure 3 


Maps showing the locations where tiger sharks (Galeocerdo cuvier), by sex, were tagged and recaptured 
in the North Atlantic Ocean between 1967 and 2018. The dashed black line indicates the 200-m contour. 


GOM=Gulf of Mexico; CS=Caribbean Sea. 


multiple years (Table 3). The largest tagged tiger sharks 
with multiple recapture events were a 214-cm-FL female 
and a 223-cm-FL male (Table 3). 


Demographic structure 


Females with length information (n=5797), based on 
combined tag and recapture data, ranged in size from 
48.0 to 457.2 cm FL. Of these, 284 females were classi- 
fied as mature. Males with length information (n=4336) 
ranged in size from 41.8 to 373.0 cm FL. Of these, 
153 males were considered mature. Sharks of unknown 
sex with length data (n=1132) ranged in size from 40.1 
to 457.2 cm FL. Of these, 107 individuals were classified 
as mature. All of these samples were used in combined 
sex analyses. 


Life stage—immature Across both sexes, 95.2% of the spec- 
imens (n=11,221, tagged and recaptured combined) were 
classified as ajuvenile tiger shark (n=9654, including YOY); 
of these, immature females (56.9%) outnumbered imma- 
ture males (43.1%) overall and in all seasons. Individuals of 
both sexes were found in the same general areas, although 
2 immature males were caught off Brazil, the southern- 
most extent of where our data were collected, 311.6 nmi 
farther south than any female was caught (Figs. 5 and 6). 


Immature animals of both sexes were caught in all sea- 
sons from Cape Hatteras into the GOM and down in the 
Bahamas and the Caribbean Sea (Suppl. Figs. 1A and 2A) 
(online only). In winter and spring, immature individuals of 
both sexes were found concentrated in these areas, but in 
summer and fall, immature sharks ranged farther north 
and east, with females traveling up to Nova Scotia, Canada, 
and males moving up to Massachusetts. 


Life stage—mature Adult females (n=284) outnumbered 
adult males (n=151) in all seasons, except in summer, 
although the largest difference was in spring. Mature 
fish of both sexes overlapped the range of immature fish; 
however, mature females had a slightly different seasonal 
distribution than mature males (Suppl. Figs. 1B and 2B) 
(online only). In fall, mature males were found up to Cape 
Hatteras, and mature females were found as far north 
as Long Island and slightly offshore. In winter, mature 
females were found as far north as North Carolina, and 
mature males were found only in the GOM and south. 
In spring and summer, mature fish of both sexes shared 
approximately the same range. 


Life stage—young of the year The majority of YOY 
(n=4606) were encountered between the U.S. waters of 
the GOM and Long Island, although they were found as 
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Table 2 


Summary of tag and recapture data by life stage for tiger sharks (Galeocerdo cuvier) from the western North Atlantic Ocean caught 
between 1963 and 2018. Tagged sharks were classified into 3 life stages: young of the year (YOY), juveniles (immature), or adults 
(mature). Mean distances traveled and times at liberty are provided with standard deviations (SDs) and maximum values. Other 
values include the number and percentage of tagged fish that traveled >500 nautical miles (nmi) and the number of tagged fish 
that crossed from the western to the eastern North Atlantic Ocean. The asterisk (*) indicates that this shark was tagged as a YOY 
and recaptured as a juvenile. 


No. of No. of Distance travelled (nmi) No. Jo No. that Time at liberty (years) 
sharks recapture a arcana aS 1k GS A >500 >500 crossed aga ET Gel et ee 
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Figure 4 


Map showing the locations where immature tiger sharks (Galeocerdo cuvier) that traveled east of 
longitude 60°W were tagged and recaptured in the North Atlantic Ocean between 1963 and 2018. 
The dashed black line indicates the 200-m contour. GOM=Gulf of Mexico; CS=Caribbean Sea. 
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Table 3 


List of tagged tiger sharks (Galeocerdo cuvier) that were recaptured multiple times, by tagging location and month and by size in fork length (FL), sex, and estimated age at 
tagging calculated by using conversions from Kneebone et al. (2008). Records of sharks tagged and caught in the North Atlantic Ocean during 1963-2018 were used in this 
study. The letter X marks the month the shark was tagged and the following months in which it was recaptured. If there is no superscript character after the X, the tagging 
and recapture events occurred in the same location. A superscript numeral 2 denotes that the shark was recaptured in the same month as it was tagged or that it was recap- 
tured twice in the same month. If the superscript 2 follows a superscript letter, both recapture events in that month occurred at the same location. The superscript letters NV, 
V,J, M,Z, C,F,S, and R denote the following recapture locations, respectively: off North Carolina, Virginia, New Jersey, Maryland, Venezuela in the Caribbean Sea, northern 
Cuba, Florida’s western panhandle, Florida’s southern Atlantic coast, and east of the Mid-Atlantic Ridge. The asterisk (*) identifies a tiger shark that was recaptured 6 years 
following the tagging year in June, approximately 260 nautical miles north-northwest of Bermuda. 
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Figure 5 


Maps showing the distribution of female tiger sharks (Galeocerdo cuvier) tagged and recaptured in the western North 
Atlantic Ocean between 1974 and 2018, by life stage, young of the year (YOY), juveniles (immature), and adults 
(mature). The dashed black line indicates the 200-m contour. GOM=Gulf of Mexico; CS=Caribbean Sea. 


far north as the Grand Banks of Newfoundland (n=1) and 
as far east as Georges Bank (n=3) (Figs. 5 and 6). In 
spring, YOY were tightly clustered from the east coast of 
Florida up to Cape Hatteras and from off the west coast 
of Florida to the Florida Panhandle (Suppl. Figs. 1C 
and 2C) (online only). In summer, YOY were up the coast 
from Cape Hatteras, as far north as Massachusetts, and 
were found south in the GOM, but by fall, they appeared 
to be retreating south and, by winter, the majority were 
between Cape Hatteras and the tip of Florida, similar to 
the distribution of YOY in spring (Suppl. Figs. 1C and 2C) 
(online only). The majority of YOY (82.0%, n=3777) were 
caught between 29°N and 34°N, with a concentration 
(65.6%, n=3023) at 29-31°N. In the GOM, YOY (n=273) 
were found primarily along the coast between 80.3°W 
and 97.0°W, although one was captured in Mexico, and in 
winter, individuals of both sexes were found farther off- 
shore (Figs. 5 and 6, Suppl. Figs. 1C and 2C [online only]). 


Neonates and mature females In May—August, 225 neo- 
nates were caught from off eastern Florida up through New 
York. Although neonates were caught in all 4 months, the 
majority were caught in July (55.5%) and August (32.9%), 
and 86.2% of them were caught between Florida and Geor- 
gia (between 29°N and 32°N) (Suppl. Fig. 3) (online only). 


Mature females were primarily caught in May (50.4%) 
and were less concentrated in distribution than neonates; 
the majority of the mature female sharks in May (66.1%) 
were evenly distributed between 29.6°N and 36.8°N 
(Suppl. Fig. 4) (online only). Of those females caught in May, 
24.2% were caught in the area of highest neonate abun- 
dance, and 43.5% of those females caught in June were 
caught in that area. In July and August, the portion of 
mature females caught in the area of neonate abundance 
decreased (to 21.4% and 20.0%, respectively). In these lat- 
ter months, mature females were more spread out and dis- 
tributed as far north as 40.9°N. No neonates were found 
in the northeast Caribbean Sea. Not enough samples were 
obtained from the GOM for neonate inquiry. 


Discussion 


The results of our analysis of 55 years of tag and recapture 
data advance our knowledge of the distribution and habi- 
tat use by life stage of tiger sharks in the North Atlantic 
Ocean. Conventional tagging relies on cooperation from 
commercial and recreational fishing industries, whether 
they participate as taggers or simply capture tagged 
sharks. Therefore, results include fishery-dependent data, 
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Figure 6 


Maps showing the distribution of male tiger sharks (Galeocerdo cuvier) tagged and recaptured in the western North 
Atlantic Ocean between 1963 and 2018, by life stage, young of the year (YOY), juveniles (immature), and adults 
(mature). The dashed black line indicates the 200-m contour. GOM=Gulf of Mexico; CS=Caribbean Sea. 


and the number of tagging and recapture events in an 
area may reflect fishing effort rather than fish distribu- 
tion. Additionally, the CSTP is based in the United States, 
likely leading to a higher proportion of recapture events 
occurring in U.S. waters of the Atlantic Ocean, the GOM, 
and the Caribbean Sea. In this study, the tagging and 
recapture data that were analyzed are a compilation of 
fishery-dependent (from commercial and recreational fish- 
eries) and fishery-independent (from research trips and 
surveys) sources. Differences in type and location of fish- 
ery, gear selectivity, and fishing method can influence the 
size and distribution of the individuals caught. For exam- 
ple, most of the fishing for all gears and industries occurred 
close to shore, leading to a higher likelihood of encounter- 
ing immature sharks (Fig. 1) (Natanson et al., 2023). These 
same limitations have been noted in other studies on 
migratory sharks (Mucientes et al., 2009; Heupel et al., 
2018; Natanson et al., 2020). Because the large amounts of 
data on such sharks are inaccessible when using other 
methods, the use of fishery-dependent methods will 
remain an important source of data. 

Because fishery-dependent data could bias results if 
fishing occurred in areas of concentrated juveniles, or of 
single-sex groups, and because the inclusion of fishery- 
independent data from scientific surveys and comparisons 
to published data from the use of fishery-independent 


methods (i.e., electronic tagging) have their own inherent 
biases, it is necessary to use multiple sampling methods 
and gear types to mitigate much of such biases. It is also 
important to note that some of the fishery-independent 
data were obtained by using sampling of fixed stations, 
possibly contributing some bias due to site selection that 
may not result in sampling across the entire habitat 
(McClelland and Sass, 2012). The ability to use data from 
both fishery-dependent and -independent sources allowed 
us to have a broader geographic range and, as a result, 
a broader biological range (the entire size range for each 
sex) for examining this species (Natanson et al., 2023). 
We observed substantial geographic overlap in the 
distribution of individuals of all life stages and of both 
sexes of the tiger shark. In particular, all life stages 
and both sexes spent considerable time between Cape 
Hatteras and northeastern Florida. The majority of our 
data are from immature sharks caught at depths less 
than 200 m. Driggers et al. (2008, fig. 3) found YOY and 
juvenile sharks in this same area as well. Our data also 
indicate that tiger sharks have some degree of seasonal 
north-south migration, a finding that corresponds with 
results of other studies that indicate similar migrations 
in the Atlantic Ocean (Lea et al., 2015; Ajemian et al., 
2020) and that indicate that tiger sharks are capable of 
both making extensive migrations and having periods of 
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residency (Sulikowski et al., 2016; McClain et al., 2022). 
Mature tiger sharks of both sexes were primarily caught 
in the GOM in winter, off the east coast of Florida and off 
North Carolina in spring, and ranged from the GOM to 
New England in summer. In fall, most catches of adults 
were off the coast from the east coast of Florida to North 
Carolina, although adult females were found as far north 
as Long Island. The observation of adult tiger sharks sea- 
sonally in the GOM is consistent with reports of mature 
tiger sharks in the GOM in spring and summer (Springer 
1938, 1940; Baughman and Springer, 1950; Clark and 
von Schmidt, 1965; Dodrill, 1977). 

Timing and location of mating cannot be determined 
with our data. Mature males and females overlap through- 
out the majority of their range in all seasons, except for 
in winter when overlap is limited and a low sample size 
for mature males in winter may bias this result. Results 
of previous studies indicate that a mating season occurs in 
spring in the GOM and the Caribbean Sea (Baughman and 
Springer, 1950; Clark and von Schmidt, 1965; Rivera-L6opez, 
1970). The Bahamas, in the vicinity of Tiger Beach, have 
been proposed as a mating area for tiger sharks in the 
Atlantic Ocean, on the basis of mating scars found on 
females in winter (Sulikowski et al., 2016). Lea et al. (2015) 
tracked mature males and found they routinely returned 
to this area in winter and suggested mating as a possible 
driver. Our data do indicate overlap of mature males and 
females in this region of the Bahamas during winter. Addi- 
tionally, our data indicate some overlap in the GOM off the 
panhandle and west coast of Florida near Apalachicola and 
Englewood, respectively, in the spring, but no overlap in 
the Caribbean Sea during spring. 

Pregnant females with embryos in various stages of 
development have been observed off the west coast of 
Florida in spring and early summer (Springer, 1938, 1940; 
Baughman and Springer, 1950; Clark and von Schmidt, 
1965; Dodrill, 1977) and off southwest Puerto Rico in the 
Caribbean Sea (Rivera-Lépez, 1970). Data from the CSTP 
indicate the presence of YOY in this area; however, there 
is no evidence of neonates, indicating that these YOY may 
move into these areas from the GOM or from the main 
body of the Atlantic Ocean after parturition. Sulikowski 
et al. (2016) suggested that the area off Tiger Beach may 
be a gestation area, given the consistent number of large 
females with an extended residency and observations of 
pregnant females. Smukall et al. (2022) reported captures 
of near-term and recently postpartum females, verified by 
ultrasonography, near the Bimini Islands in the Bahamas 
during late spring and summer. Additionally, Smukall et al. 
(2022) reported frequent catches of young-of-the-year tiger 
sharks, including neonate-size individuals, but acoustic 
data indicate limited residency, potentially meaning that 
YOY moved out of the area soon after parturition. 

Because near-term pregnant females are found in the 
GOM and Atlantic Ocean, it is probable that they birth, 
and possibly gestate, in their respective water bodies, 
although we do not have direct data from the GOM to cor- 
roborate it for that area. Our data indicate that the high- 
est concentration of neonates between March and August 
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overlap the locations of mature females in the waters off 
the East Coast of the United States at this time (Springer, 
1938; Clark and von Schmidt, 1965). The presence of post- 
partum females between New Jersey and southern New 
England during June-August indicates a northward 
movement of mature females after birthing (Natanson 
etale. 2023) 

Seasonal north-south migrations were undertaken by 
tiger sharks in the YOY, juvenile, and adult stages. The YOY 
concentrate between the east coast of Florida and Cape 
Hatteras throughout the year, expanding north in summer 
and fall (Suppl. Figs. 1C and 2C) (online only). Natanson 
et al. (1999) reported concentrated catches of neonate and 
YOY tiger sharks on the continental shelf between 29°N 
and 33°N. Driggers et al. (2008), using data on catch per 
unit of effort of neonates, expanded this region to 27—35°N, 
with the area of 31—33°N as the most important. The major- 
ity of YOY (82.0%) in this study were caught between 29°N 
and 34°N, with a concentration between 29°N and 31°N 
(65.6%). In contrast, Driggers et al. (2008) found the high- 
est concentration of YOY tiger sharks at 32°N in waters off 
the East Coast. The differences in these areas of high con- 
centration for YOY between our study and Driggers et al. 
(2008) may be a result of timing and type of sampling. Data 
from this study were collected year-round and are a combi- 
nation of 55 years of fishery-dependent and -independent 
(fixed-station and random-stratified surveys and targeted 
research sampling) data, and the data from Driggers et al. 
(2008) are temporally discrete (August-September) and 
based on 13 years of a single fishery-independent survey 
(with a random-stratified design). Additionally, the concen- 
trations may shift depending on optimal environmental 
conditions for each stage in a given year. 

Shark nursery areas can be detected by using the pres- 
ence of gravid females with near-term embryos, neonates, 
and small juveniles for viviparous species (Castro, 1993). 
Although presence data can be used to detect nursery 
habitat, isolated events do not provide information on the 
importance of this habitat in supporting juvenile shark 
populations (McCandless, 2007). Monitoring through 
long-term fishery-independent surveys, conventional 
mark-recapture work, and acoustic telemetry studies can 
help develop a better picture of how nursery habitat is 
used (McCandless, 2007). With fishery managers in mind, 
Heupel et al. (2007) developed a quantitative definition 
for identifying important shark nursery areas by using 
3 criteria for neonates and YOY: 1) these juveniles are 
more commonly encountered in the area than other areas, 
2) they have a tendency to remain in or return to the area 
for extended periods, and 3) the area is repeatedly used 
across years. 

On the basis of these criteria, the tagging and recapture 
data used in our study indicate that the area from off the 
east coast of Florida north of Cape Canaveral to off Jekyll 
Island in Georgia (between 29°N and 31°N) provides nurs- 
ery habitat for tiger sharks. Data from multiple recapture 
events indicate that this area has concentrations of neo- 
nate- and YOY-size sharks, that these sharks spend sev- 
eral months there following parturition that occurs within 
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the area or following migration to the area, and that these 
sharks repeatedly use this area across years. Expanding 
this region to Cape Island, South Carolina (33°N), incor- 
porates data from Driggers et al. (2008) indicating the 
importance of this region for both neonates and YOY. Our 
data on neonates, YOY, and multiple recaptured imma- 
ture tiger sharks, combined with data from Driggers et al. 
(2008) and McClain et al. (2022), support the existence of 
separate, broad birthing areas and nursery habitats in 
the GOM and in waters off the East Coast (Heupel et al., 
2007; Driggers et al., 2008). Knowledge of defined areas 
of concentration informs management decisions regarding 
refinement of essential fish habitat and spatial manage- 
ment and conservation goals to protect the vulnerable life 
stages of this species (Heupel et al., 2007, 2018; Kinney 
and Simpfendorfer, 2009). The area of concentration for 
YOY in this study encompasses the range where neonates 
were concentrated. 

This report is the first to document that immature tiger 
sharks undertook long-distance movements from the 
western to the eastern Atlantic Ocean. These data indicate 
movements of tiger sharks that are much broader and of 
longer distances and taken at smaller sizes than previ- 
ously known for this species. Although results from other 
studies indicate that adult males and females and sub- 
adult (no longer a juvenile but not yet an adult) females 
undertake long-distance movements (Hammerschlag 
et al., 2012; Lea et al., 2015; Ajemian et al., 2020), the out- 
comes of many studies indicate that juvenile tiger sharks 
cannot undertake long-distance migrations because of 
inefficient swimming in the more anguilliform body of 
the young tiger shark (Lea et al., 2015, 2018; Aines et al., 
2018) or because of the effect of environmental factors on 
body size (i.e., small sharks are less tolerant of low sea- 
surface temperatures) (Lea et al., 2018). 

Our data indicate that YOY use coastal waters of the 
United States fairly continuously from Louisiana to Long 
Island, depending on the season, and are found offshore 
and north to New England in summer. Individual tiger 
sharks of both sexes tagged as YOY traveled >500 nmi, 
indicating that this life stage is capable of movements of 
longer distances than previously thought. For example, 
one YOY, was tagged in February and recaptured 6 months 
later southeast of Newfoundland, having traveled nearly 
1300 nmi. Another was tagged in November and recap- 
tured after 1 year in the Caribbean Sea off Venezuela, 
having traveled nearly 1400 nmi. By the juvenile stage, 
the range of tiger sharks becomes Atlantic basin-wide. The 
coastal range extends from New England down and into 
the Caribbean Sea. The north-south migration along the 
Atlantic coast of the northeastern United States taken by 
immature fish is complicated by offshore and transatlantic 
movements. These sharks move north in summer past New 
England to the east; some remain in coastal waters, and 
others venture offshore presumably with the Gulf Stream 
and as far north as the Grand Banks of Newfoundland. 
This northward offshore movement with the Gulf Stream 
has been documented for adult males (Lea et al., 2015, 
2018) and for adult and subadult females (Hammerschlag 
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et al., 2012) by using satellite tracking data, and results 
of our study indicate that immature tiger sharks of both 
sexes also follow this pattern. Additionally, our data indi- 
cate movement of immature tiger sharks, including YOY, 
into and out of the Caribbean Sea from both the main body 
of the Atlantic Ocean and the GOM. 

It has been suggested that ontogenetic movements into 
more oceanic habitats are possibly related to foraging (Lowe 
et al., 1996; Lea et al., 2015), and more recently McClain 
et al. (2022) suggested a reproductive component to migra- 
tion. On the basis of their electronic tagging data, Lea et al. 
(2018) further suggested an ontogenetic shift in foraging 
strategy and diet that leads to longer migrations by adult 
males. Although a foraging component may certainly drive 
these migrations, our data indicate that immature tiger 
sharks make extensive migrations throughout the Atlantic 
Ocean, refuting the ontogenetic component of the basis for 
longer distance migration. Long migrations are evident 
even for young juveniles; for example, the transatlantic 
migration of over 3600 nmi of a tiger shark tagged as a 
YOY and recaptured 2 years later. 

In conjunction with the data we present on juveniles 
crossing the Atlantic Ocean, and data from other studies, 
we suggest that a northeasterly track could be the route 
that these sharks take to the eastern Atlantic Ocean. 
Following the major current system in the North Atlan- 
tic Ocean, this track starts with the Gulf Stream, mov- 
ing north and east with the North Atlantic Drift and then 
south to the Canary Current before following the North 
Equatorial Current back to the WNA as has been previ- 
ously postulated for the migration of the blue shark (Casey, 
1985). Afonso et al. (2017) suggested an equatorial route 
for juvenile tiger sharks linking the eastern and west- 
ern Atlantic Ocean. Additionally, Domingo et al. (2016), 
using observer data, documented that tiger sharks in the 
North Atlantic Ocean were more broadly distributed than 
previously thought. A limited number of mature sharks 
of both sexes (n=10; 2 males, 7 females, and 1 individual 
of unknown sex) in this study made movements greater 
than 500 nmi within the WNA, although none crossed the 
North Atlantic basin. Although mature individuals have 
not yet been documented to travel across the North Atlan- 
tic basin, they have been repeatedly shown to have the 
ability to travel the long distances required for this move- 
ment (Hammerschlag et al., 2012; Ferreira et al., 2015; 
Lea et al., 2015). The lower percentage of mature sharks 
in our sample may also be a factor in this result. Much 
of the data are limited to immature sharks crossing the 
North Atlantic basin, but these movements and distribu- 
tions strongly indicate that there is only one population in 
the North Atlantic Ocean. 

Conventional tags are significantly less costly than elec- 
tronic tags; therefore, more individuals of a species can 
be tagged with them at various times, locations, and life 
stages. Additionally, long-term studies with increased 
numbers of recapture events can enhance the body of 
knowledge. Sample size is important in interpreting 
results from tagging studies, and an increase in sample 
size allows for a larger range and scope of questions that 
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can be answered (Sequeira et al., 2019). When examining 
tagging data for population trends and movements, an 
increased number of tags makes it possible to see vari- 
ation in movements between individuals, life stages, or 
sexes (Sequeira et al., 2019); it is therefore important to 
tag individuals of all sizes and sexes to obtain a complete 
overview of the movements of a species. Additionally, the 
duration of tag deployment and the location of release of 
tagged fish can influence movements for certain shark 
species (Rooker et al., 2019). It is apparent from studies of 
a variety of species, including the tiger shark, that individ- 
uals can be resident or migratory for given periods of time 
and stages of maturity. 

Despite the number of tiger sharks tagged in the CSTP, 
it took 34 years and 5832 tiger sharks to be tagged before 
the first tiger shark tagged in the WNA was recaptured 
in the ENA. It is certainly doubtful that we have seen 
all possible scenarios, particularly with mature sharks. 
These recapture events are important not only in that the 
tiger shark can cross the North Atlantic Ocean, therefore 
indicating potential mixing within the population, but 
in that immature sharks are capable of this movement. 
Pooling all tagging data on the tiger shark in the WNA 
creates a better picture of the movements of this species 
than examining those studies individually. For example, 
Lea et al. (2015, 2018) stated that only mature males 
make long-distance migrations and that movements into 
and out of the Caribbean Sea were not observed; how- 
ever, with the addition of Hammerschlag et al. (2012), 
Rooker et al. (2019), and Ajemian et al. (2020), it is appar- 
ent that females also travel long distances and exchange 
between these water bodies. Combining data from con- 
ventional and electronic tagging studies can provide a far 
more robust, comprehensive picture of the movements of 
a species and of the drivers for those movements than 
the representation formed when using data from just one 
type of tagging, leading to advancements in management 
measures (Hammerschlag et al., 2022) to better support 
sustainable fisheries. 


Conclusions 


We examined seasonal movement patterns of the tiger 
shark on the basis of sex and life stage and have added 
to the literature describing seasonal north-south move- 
ments of the tiger shark. Additionally, our data indicate 
that immature tiger sharks tagged in the WNA make long- 
distance migrations to the ENA. We propose that there is a 
potential route that runs up the coast of North America in 
the WNA and crosses over to the ENA, revealing potential 
mixing within the population not indicated by data from 
electronic tagging to date. We also identified a broad nurs- 
ery area and potential birthing grounds off the East Coast 
of the United States. Tagged tiger sharks crossed multiple 
international jurisdictions, indicating that management 
measures need to have a basin-wide scale. The results of 
this project highlight the need for the use of more long- 
term, conventional tagging along with electronic tagging 
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to fully understand the movements of shark species and to 
adequately implement management strategies. 


Resumen 


Los movimientos del tiburén tigre (Galeocerdo cuvier) 
durante las etapas de su vida, son en gran medida descono- 
cidos. Sin embargo, es necesario conocerlos para determinar 
su habitat esencial y las practicas sostenibles de manejo 
pesquero. En un esfuerzo por dilucidar la distribuci6n y los 
movimientos del tiburén tigre, analizamos datos de marcaje 
dependientes e independientes de la pesqueria (tamano de 
la muestra [n]=10,516) y registros de recaptura (n=762) 
de tiburones tigre capturados en el océano Atlantico norte 
entre 1963 y 2018. Se examin6o la distribucién estacional 
del tiburén tigre por etapa de vida—juveniles del ano, juve- 
niles (inmaduros) y adultos (maduros)—y se utilizaron 
patrones de distribucién para identificar posibles habitats 
de crianza. Se capturaron tiburones tigre en una amplia 
zona desde los Grandes Bancos de Terranova (Canada) 
hacia el sur, en Brasil y desde la costa hasta alta mar y 
hacia el este del océano Atlantico Norte. Se observaron mov- 
imientos estacionales de norte a sur en todas las etapas de 
la vida, y se descubrié que 14 tiburones inmaduros habian 
migrado del océano Atlantico Norte occidental al oriental. 
Se identificé una amplia zona de crianza y una posible zona 
de nacimiento en la plataforma continental entre Florida y 
Georgia con base en recurrencia de neonatos en verano a lo 
largo de los afios y en la recaptura de multiples tiburones 
tigre juveniles del aio marcados en el mismo lugar durante 
un periodo de al menos 2 anos. 
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Abstract—The Atlantic sturgeon (Aci- 
penser oxyrinchus oxyrinchus) has been 
listed under the Endangered Species 
Act since 2012, with the Chesapeake 
Bay distinct population segment listed 
as endangered. We tracked the tim- 
ing of occupancy and movement in 
Chesapeake Bay of adult Atlantic stur- 
geon from the York River to best iden- 
tify when Atlantic sturgeon are likely 
to interact with anthropogenic threats. 
We monitored 84 adult (40 male and 
44 female) Atlantic sturgeon from 
August 2013 through January 2020 by 
using acoustic telemetry. Both spawn- 
ing and non-spawning fish regularly 
utilized Chesapeake Bay, with females 
and males arriving as early as 27 Feb- 
ruary (when the mean water tempera- 
ture was 7.7°C) and 4 March (6.5°C) and 
departing as late as 24 January (6.3°C) 
and 27 January (6.5°C), respectively. 
Peak occupation of the bay by Atlan- 
tic sturgeon occurred from 1 April to 
31 August and again from 15 October to 
1 December. Females of above average 
size (>1880 mm in fork length) spent 
significantly longer in the bay (>118 d) 
than smaller females and all males 
before spawning; therefore, the females 
capable of producing the most eggs were 
disproportionately exposed to anthro- 
pogenic threats. Although changes in 
arrival and departure dates are not sta- 
tistically significant, during the 7 years 
of this study, both males and females 
generally arrived earlier and departed 
later each year than the previous year, 
increasing residency in Chesapeake 
Bay by a month in that time. 
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The National Marine Fisheries Service 
listed 5 distinct population segments 
(DPSs) of Atlantic sturgeon (Acip- 
enser oxyrinchus oxyrinchus) under 
the Endangered Species Act in 2012 
(Federal Register, 2012a 2012b). Since 
this listing, new spawning populations 
have been discovered (Hager et al., 
2014; Balazik et al., 2017; Farrae et al., 
2017; Savoy et al., 2017). It now appears 
that some rivers support 2 genetically 
different spawning populations and 
that rivers identified within the geo- 
graphic area of a DPS may have a sea- 
sonal spawning group of sturgeon that 
are more closely related to sturgeon 
spawning in other rivers outside of the 
area of that DPS (White et al., 2021). 
Of the 5 DPSs in the United States, the 
Gulf of Maine DPS was listed as threat- 
ened under the Endangered Species 
Act, and the New York Bight, Chesa- 
peake Bay, Carolina, and South Atlan- 
tic DPSs were listed as endangered. At 
the time of listing, the Chesapeake Bay 
DPS had one known spawning popula- 
tion: the population in the James River, 
which at the time had been confirmed 
to spawn after summer as temperatures 
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U.S. Department of the Navy 
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Washington Navy Yard 
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decline (Balazik et al., 2012a). Like the 
York River, the James River is located in 
Virginia. Since the listing in 2012, adult 
Atlantic sturgeon have been confirmed 
to spawn in the York River (Hager et al., 
2014); spawning is likely occurring in 
the Nanticoke River, located on the 
Delmarva Peninsula, and possibly in the 
Rappahannock River, located in east- 
ern Virginia (ASMFC, 2017). Results of 
genetic analyses of the populations in 
the James and York Rivers indicate that 
these 2 populations are significantly dif- 
ferentiated (White et al., 2021). 

Atlantic sturgeon are long-lived, wide- 
ranging, late-maturing, iteroparous, 
anadromous fish that spawn intermit- 
tently (Smith, 1985; Bemis and Kynard, 
1997; Dadswell, 2006; NMFS, 2007; Peter- 
son et al., 2008; Hager et al., 2020). Adult 
and migratory juvenile Atlantic sturgeon 
spend significant portions of each year 
along the East Coast for feeding, growing, 
overwintering, and migrating (Dovel and 
Berggren, 1983; Smith, 1985; Bain, 1997; 
Stevenson, 1997; Wirgin et al., 2015). 
These individuals undertake seasonal 
migrations to utilize inshore and offshore 
habitat (Ingram et al., 2019). 
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In Chesapeake Bay, there have been no peer-reviewed 
reports of the timing or windows of occupancy of Atlantic 
sturgeon. The objective of this study was to identify that 
timing and the movement of adult Atlantic sturgeon within 
Chesapeake Bay. We evaluated the migrations of males 
and females separately in terms of seasons, years, calen- 
dar days, photoperiod, and temperatures. Understand- 
ing the periods of occupancy is essential for the National 
Marine Fisheries Service to protect and ultimately recover 
this endangered species. 


Materials and methods 


We identified the entrance to Chesapeake Bay as an invis- 
ible line from Fort Story, Virginia, to Fisherman’s Island 
National Wildlife Refuge near Cape Charles, Virginia, 
because this location was the most restricted point at the 
mouth of the bay (Fig. 1). We determined that a fish had 
left Chesapeake Bay if it was detected on receivers beyond 
the boundary denoted by this line or if it entered the 
mouth of the James, York, or Rappahannock River (Fig. 1). 
At locations where passive acoustic receivers (VR2W, 
Innovasea Systems Inc.', Boston, MA) were deployed to 


1 Mention of trade names or commercial companies is for identi- 
fication purposes only and does not imply endorsement by the 
National Marine Fisheries Service, NOAA. 
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detect fish implanted with acoustic tags, water tempera- 
tures were collected with HOBO U12-015 data loggers 
(Onset Computer Corp., Bourne, MA) zip tied to the acous- 
tic receivers or from the First Landing (near Cape Henry, 
Virginia) and York Spit (at the mouth of the York River 
near Perrin, Virginia) monitoring locations of the NOAA 
Chesapeake Bay Interpretive Buoy System (data avail- 
able from website, accessed January 2022). 

Atlantic sturgeon were captured, implanted with acous- 
tic transmitters (or tags), and released while on spawning 
runs in the Pamunkey and Mattaponi Rivers, the 2 main 
tributaries to the York River. For our analysis, individuals 
were considered to be on spawning runs if they met a 
2-part test: 1) the fish moved a minimum of 20 km above 
the saltwater interface, and 2) it spent at least 2 weeks in 
fresh water (Kahn et al., 2019; Hager et al., 2020; Kahn 
et al., 2021). All tagged Atlantic sturgeon were in the adult 
size range (1250-2272 mm in fork length [FL]; Grunwald 
et al., 2008; Waldman et al., 2019; Kahn et al., 2019, 2021) 
and were captured in upstream freshwater locations on 
confirmed spawning grounds (Kahn et al., 2019; Hager 
et al., 2020). 

Atlantic sturgeon were captured by using gill nets 
made of mesh with stretch measures of 23-36 cm from 
late July through mid-October during 2013-2019 (see 
Kahn et al., 2019) in the Pamunkey and Mattaponi Rivers. 
Our research methods followed the permit requirements 
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Figure 1 


Map of the mouth of Chesapeake Bay showing the locations of passive acoustic receivers and the boundaries (white 
lines) identified to define the extent of Chesapeake Bay for documenting arrival, departure, and residency of tagged 
Atlantic sturgeon (Acipenser oxyrinchus oxyrinchus) in the bay. Fish were monitored with acoustic telemetry from 


August 2013 through January 2020. 
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and protocols of the Endangered Species Act (Kahn and 
Mohead, 2010). 

Results of the genetic analysis of a previous study (White 
et al., 2021) indicate that fish from the York River are 
distinct from fish from other populations along the East 
Coast. However, some of the fish implanted with trans- 
mitters by our team and used in this study were fish natal 
to the James, Hudson, and Savannah Rivers, as well as to 
Albemarle Sound (Kazyak et al., 2021). Therefore, the occu- 
pancy timing presented herein is weighted toward adults 
of the population of the York River but is representative of 
movements of fish natal to 3 different DPSs. 

The sex of individuals was determined by applying pres- 
sure to the ventral surface, moving from the anterior to 
posterior ends, ending at the vent. This technique typically 
produced milt from males (>90% of the time), and females 
on occasion produced eggs (<20% of the time). Gravid 
females were usually confirmed to be female when trans- 
mitters were implanted. Sex was routinely confirmed upon 
recapture (Kahn et al., 2021). Ultimately, sex was confirmed 
for every fish tagged and recaptured in this study. 

During this study, 84 adult Atlantic sturgeon (40 males 
and 44 females) were implanted with V16P-4H, V16P-6x, 
or V16-6x acoustic transmitters (Innovasea Systems Inc.), 
with each tag weighing no more than 17.3 g. Transmitters 
produced a 69 kHz signal every 70-150 s and had a min- 
imum life span of 6 years (11 transmitters) and a maxi- 
mum life span of 10 years (73 transmitters). Transmitters 
were placed into incisions of 3-4 cm that were made most 
often between the 3rd and 4th ventral scutes anterior to 
the anal fins. The incisions were closed by using Vicry] dis- 
solvable sutures (Ethicon Inc., Raritan, NJ). 

The implanted transmitters were passively detected 
from August 2013 through January 2020. No females were 
tagged with acoustic transmitters in 2013; therefore, emi- 
erations in 2013 and immigrations in 2014 were recorded 
for only males. Through January 2020, 100 acoustic receiv- 
ers in the York River system, Chesapeake Bay, and Atlantic 
Ocean remained in place (Fig. 1). They were serviced and 
their data was downloaded either monthly or every other 
month depending on their location. Receivers in rivers, 
Chesapeake Bay, and nearshore waters of the Atlantic 
Ocean were downward-facing and deployed within 6 m of 
the surface, and offshore receivers were anchored to the 
seafloor and recovered with an acoustic release. 

When monitoring movements of Atlantic sturgeon, we 
assumed acoustic transmitters were in fish that were 
behaving normally. We also performed range tests by 
using stationary transmitters and receivers at variable 
distances in Chesapeake Bay near the mouth of the James 
River. We used transmitters of the same size for the mon- 
itoring in this study and for the range tests to minimize 
error, but detection distances can also be affected by envi- 
ronmental conditions, such as tides and sea state (Mathies 
et al., 2014). A single stationary transmitter with a 10-s 
ping rate was deployed, and receivers were deployed every 
100 m away from this transmitter. When receivers were 
placed at a depth of approximately 10 m, the distance at 
which over 90% of signals from the test transmitter were 


detected under optimal conditions was 1.3 km; in contrast, 
for receivers deployed at depths less than 2 m under rough 
conditions, this detection rate extended out only 200 m. 
Under various conditions over a 1-week period, receivers 
within 700 m of the transmitter detected at least 90% of 
transmissions. For the purposes of understanding when 
Atlantic sturgeon were present in Chesapeake Bay, we 
defined presence as multiple detections on 1 receiver or 
detections on 2 neighboring receivers on the same day. 
Therefore, any variability in detection distance would lead 
to underestimation of the timing of arrival and departure. 
Generally, we expected to detect fish within 700 m of a 
receiver. 

Atlantic sturgeon belonging to the spawning population 
in the York River spawn only during the fall (August-— 
October) of each year (Kahn et al., 2019; Hager et al., 
2020). The times of first arrival to Chesapeake Bay from 
the Atlantic Ocean, of departure from Chesapeake Bay 
into tributary river systems, of entry into Chesapeake Bay 
from those same systems, and of departure from Chesa- 
peake Bay to the Atlantic Ocean were recorded for each 
transmitter implanted in a fish as ordinal dates and as 
photoperiod in hours, minutes, and seconds from sunrise 
to sunset. We report the timing of these movements for 
each sex but considered mean, median, and timing for 
different proportions of transmitters (10%, 20%, etc.). In 
any given year, the general pattern of migration followed 
a normal distribution once it began; therefore, we report 
the earliest arrival and latest departure because 1) not 
all fish in the population are tagged with acoustic trans- 
mitters and not all tagged fish were acoustically detected, 
meaning that the timing is an approximation of actual 
first arrival or last departure, and 2) the purpose of this 
analysis was to identify periods of presence and absence 
from Chesapeake Bay. However, results from analyses of 
statistical differences are presented as mean values. We 
recorded the temperature from the nearest source (HOBO 
data logger, acoustic receiver, or buoy). 

Most individuals left Chesapeake Bay in the summer 
by initiating a spawning run. Non-spawning fish gener- 
ally remained near the mouth of the bay and were con- 
sistently detected throughout the summer, and spawning 
fish moved farther into the bay and had periods during 
which they were not detected as they moved through the 
areas without receiver coverage (Fig. 1). We conservatively 
calculated residency times in the bay from the date of first 
detection to the date of last detection, such that data for 
an animal detected leaving the mouth of the York River 
in the fall and not detected again were aggregated with 
data for other tagged fish leaving the river and used to 
calculate an average date of entry into the bay from the 
river, but we did not use that single data point to estimate 
bay residency. However, the dates of first or last detec- 
tion, when used to calculate residency, were not necessar- 
ily representative of the beginning or end of occupancy 
because receiver coverage was incomplete. These periods 
of presence in Chesapeake Bay were analyzed by sex, by 
length and sex (because adult Atlantic sturgeon are sexu- 
ally dimorphic), and by year. 
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Each Atlantic sturgeon was classified as either above or 
below the average length for its sex. Measurements were 
made by using a 2.44-m PVC T-board with a 3-m tailor’s 
tape measure. The mean lengths of males and females 
were 1518 and 1880 mm FL (Kahn et al., 2019). Because 
fish generally grew through time, when adults were recap- 
tured, they occasionally had moved from below average to 
above average length. For example, Fish 13-015, a male, 
was captured 6 times during the 7 years of this study (in 
2013-2017 and 2019), ranging in size from 1490 to 1548 
mm FL. This individual was categorized as below average 
the first 2 years and above average the last 5 years. In 
the event that a fish was not captured in the year of the 
analysis, its most recent measured length was used for 
categorization. 

Differences in arrival and departure times (ordinal 
dates) and temperature were assessed by using analysis 
of variance in statistical software R (vers. 4.2.2; R Core 
Team, 2022). The significance level was 0.05. When sig- 
nificant differences were identified, Tukey’s honestly sig- 
nificant difference test was used to identify significantly 
different variables. Data were organized in a table identi- 
fying sex, size category, year, bay entry, river entry, river 
departure and arrival, bay departure, photoperiod for each 
day of arrival and departure, and days spent in the bay 
before and after each spawning period for each tagged 
individual. Because photoperiod and day are correlated, 
only day of the year is discussed for ease in identifying 
when Atlantic sturgeon were present. 

Sampling was conducted under permits from the Virginia 
Marine Resources Commission, the Virginia Department 
of Wildlife Resources (formerly the Virginia Department 
of Game and Inland Fisheries), and the National Marine 
Fisheries Service (permit nos, 16547 and 19642). 


Results 


Both spawning and non-spawning fish regularly utilized 
Chesapeake Bay, starting as the water began to warm in 
the spring and ending as it cooled in the fall. The earli- 
est arrival of a tagged Atlantic sturgeon in Chesapeake 
Bay from the Atlantic Ocean was detected on 27 February 
for females and on 4 March for males, corresponding to 
daily mean water temperatures of 7.3°C and 6.5°C at the 
locations of first detection upon entry, respectively. Males 
and females generally did not arrive until water tempera- 
tures were above 6°C, except for 3 females, who arrived 
in water as cool as 2.9°C. The latest date of arrival of a 
fish implanted with an acoustic tag to the bay from the 
ocean, with the exception of the arrival of a male outlier 
on 22 August, was recorded on 17 July for males and on 
11 July for females, corresponding to daily mean water 
temperatures of 27.9°C and 26.0°C at their locations of 
first detection upon entry, respectively. In all years, at 
least 1 Atlantic sturgeon remained in the bay after its 
arrival, even during the spawning event, until departure 
back to the ocean. The latest departure (last detection) 
from the bay to the ocean after spawning was recorded 
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on 24 January for females and on 27 January for males, 
corresponding to daily mean water temperatures of 6.3°C 
and 6.5°C at the locations of last detection, respectively. 
Atlantic sturgeon generally left Chesapeake Bay for the 
ocean when water temperatures were approaching a mini- 
mum of 6°C, with the exception of a large male that left for 
the ocean when the daily mean water temperature at its 
location of last detection was 3.4°C. 

The majority of fish were detected moving through 
Chesapeake Bay on spawning runs. Spawning males were 
detected in the bay during 78.7% (244 detections of the 310 
times males are known to have moved through the bay to 
reach the river when they were detected) of their spawning 
runs, and spawning females were detected during 81.2% 
(108 detections of the 133 times females were known to 
be present) of theirs. Non-spawning fish did not always 
return to Chesapeake Bay, but when they did males were 
detected during 67.7% (23 detections of the 34 times males 
were known to be present) of their non-spawning years 
and females detected during 55.0% (61 detections of the 
111 times females were known to be present) of theirs. The 
detection rate between years varied (Table 1), but the loca- 
tions of the passive acoustic receivers we used to monitor 
migrations did not. 

Males and females were detected in Chesapeake Bay 
during January or February 1.5% (5 of 344 males) and 
1.2% (3 of 244 females) of the time, as both their arrival 
and departure between the bay and the ocean pushed into 
winter. There were no significant interannual differences 
in date of entry to the bay from the ocean in the spring 
or in date of departure from the bay to the ocean in the 
fall for males or females (Table 2). However, only 5 of 
77 (6.5%) returning fish (male or female) were detected 
in the bay before 1 April in 2014 (0 fish), in 2015 (1 fish, 
when the mean daily temperature near the detection loca- 
tion was 10.2°C), and in 2016 (4 fish, when temperatures 
were 11—12.6°C), but in 2018 and 2019, the first fish were 
detected entering the bay in February. 

Further illustrating this shift in arrival times, in 2018 
and 2019, 30.8% of all adults (33 of 107 fish) had entered 
the bay by 10 April. The first detections in February of 
2018 and 2019 were of the same female that arrived when 
the mean temperature in the bay was 2.9°C and 7.38°C. 
Likewise, in each of 2013, 2014, 2015, and 2016, cumula- 
tively only 3 fish (of 161 individuals, 1.9%) were detected 
still out-migrating through the bay toward the ocean 
during December but, by 2019, 10 fish (of 47 individu- 
als, 21.3%) were detected in the bay in December, with 
2 fish remaining in the bay until January. Of those fish 
detected in December 2019 and January 2020, 9 individ- 
uals were tagged with acoustic transmitters in 2013 and 
2014, with the other fish tagged with acoustic transmit- 
ters in 2015; therefore, the increase in sturgeon presence 
in December and January following the 2019 spawning 
season was not just a case of more Atlantic sturgeon being 
tagged with acoustic transmitters in later years. This 
apparent extension of the use of Chesapeake Bay was not 
significant, and only males had any significant differences 
in interannual behavior (Table 2). 
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Table 1 


Numbers of Atlantic sturgeon (Acipenser oxyrinchus oxyrinchus) of each sex tagged and the numbers and proportions of males 
and females that were detected within Chesapeake Bay at some point during a season, either before (pre) or after (post) spawning 
from 2013 through 2019. Atlantic sturgeon typically do not spawn every year; therefore, values are provided for both spawning 
and non-spawning fish. NA=not applicable because no tagged male or female Atlantic sturgeon were detected during the specified 


season. 


Spawning 


Males Females 


Detected in 


bay 
No. 


Season tagged 
Post 12 
Pre 9 
Post af 
Pre 24 
Post 24 
Pre 22 
Post 22 
Pre 28 
Post 29 
Pre aa 
Post val 
Pre 28 
Post 31 


The average date of arrival to Chesapeake Bay from 
the Atlantic Ocean for females was 12 April and for males 
was 22 April (P=0.005), correlating to mean temperatures 
at the locations of first detection upon entry into the bay 
of 11.8°C and 13.5°C (P=0.01), respectively. Fish of both 
sexes spent prolonged periods in Chesapeake Bay from 
entry to initiating a spawning run but took a quicker, 
more direct route out of the bay after spawning (Fig. 2). 
Females spent an average of 96.98 d (range: 27-184 
d) in the bay before spawning, compared with 82.86 d 
(range: 3-160 d) for males (P=0.051). Males entered their 
spawning river from the bay on average on 13 July (when 
the mean temperature in the bay was 25.6°C; range: 
11 April—3 September), and females entered on average 
on 15 July (24.8°C; range: 30 March—7 October). Females, 
on average, left the York River and reentered the bay on 
24 October (19.4°C; range: 17 September—19 December), 
and males reentered on 25 October (18.8°C; range: 
19 September—21 January). Males left the bay for the 
ocean after spawning on 13 November (15.4°C; range: 
25 September—27 January), on average, and the average 
departure date for females was 9 November (16.6°C; range: 
3 October—24 January). Arrival and departure dates rela- 
tive to average length are provided in Table 3. 

After spawning, males spent on average 20.51 d (range: 
1—89 d) in Chesapeake Bay, compared with 14.79 d (range: 
1-72 d) for females, as they out-migrated (P=0.022, for 
males versus females; P=0.11, when considering relative 
length of each sex). Results from assessment by period, 


Detected in 


Non-spawning 
Males Females 


Detected in Detected in 


bay bay bay 


% 


Z 
> 


0.33 
0.00 
1.00 
0.75 
0.66 
0.33 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
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either before or after spawning, indicate that females spent 
significantly longer in Chesapeake Bay only during the 
prespawning period in 2015 (P=0.02) and 2016 (P=0.038), 
spending 129.00 d (range: 116-142 d) and 115.30 d (range: 
74-155 d) in the bay when males spent an average of 
72.89 d (range: 26-134 d) and 71.92 d (range: 5-144 d) 
there. Differences between size groups and sexes were 
detected for bay entry day, river entry day, river entry tem- 
perature, and prespawning days (Table 4). 


Discussion 


The Atlantic sturgeon is an endangered species in Chesa- 
peake Bay, and the best mitigation to promote its conserva- 
tion is limiting potentially harmful interactions by identifying 
the times Atlantic sturgeon are present in Chesapeake Bay 
and by planning potentially harmful activities accordingly. 
Females were recorded in the bay during prespawning peri- 
ods from 27 February to 7 October. Males entered the bay as 
early as 4 March and stayed during the spring and summer 
until as late as 3 September. With the lone exception of the 
spawning initiation on 7 October, the departure of Atlan- 
tic sturgeon from Chesapeake Bay to rivers and reentry 
into Chesapeake Bay from rivers did not overlap. In every 
year of this study, females were in the bay from 28 March 
to 26 November, and males were in the bay from 3 April to 
23 November (Fig. 2). Females entered Chesapeake Bay from 
the Atlantic Ocean earlier than males, with the smallest 
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Table 2 


Comparison of annual averages of the days of arrival and departure of tagged 
Atlantic sturgeon (Acipenser oxyrinchus oxyrinchus) into or out of Chesapeake Bay 
or a river, temperatures in the bay and river at arrival and departure, and number 
of days tagged fish spent in the bay before and after spawning, by sex. Atlan- 
tic sturgeon were monitored with acoustic telemetry from August 2013 through 
January 2020. Letters denote statistically significant differences (P<0.05), with 
values having like letters not being significantly different from those within the 
compared group. Arrival and departure days are given as days of the year. 


Males before spawning 


Bay entry River 
day entry day 


189.107” 
187.15? 
186.22? 
194.092» 
212.09° 


115.85 
115.11 
104.00 
114.85 
107.41 


Males after spawning 


River exit Bay exit 
day day 


289.002 
296.252%° 
300.45°° 
295.537 
300.66"° 
306.44° 


309.14 
311.00 
322.36 
320.00 
319.72 
319.70 


Bay entry 
temp (°C) 


14.9 
13.8 
12.5 
12.4 
13.7 


River exit 
temp (°C) 


20.1? 
18.12” 
18.27» 
2022 
18.62 
siycatg 


River entry 


temp (°C) 


25.62” 
24,92» 
24.5? 
26.12” 
26.7° 


Bay exit 
temp (°C) 


15.4 
17.0 
14.9 
15.0 
13.6 
14.9 


Females before spawning 


Bay entry River 
day entry day 


104.67 222.00 
104.14 197.20 

88.47 176.25 
108.16 202.82 
103.15 192.59 


Bay entry 
temp (°C) 


13.2 
13.1 

9.4 
10.9 
13.1 


River entry 


temp (°C) 


Zyl 
22.9 
24.2 
26.8 
23.0 


Females after spawning 


River exit Bay exit 
day day 


297.25 
287.67 
289.67 
295.28 
292.93 
302.85 


313.71 
310.67 
307.18 
310.41 
311.35 
317.96 


females returning the earliest (Table 4). When spawning, 
females of below average size (<1880 mm FL) entered the 
York River the earliest, all males entered roughly at a simi- 
lar intermediary time, and then significantly later, females of 
above average size moved into the river (Table 4). The larg- 
est females with the most reproductive potential (Mitchell 


River exit 
temp (°C) 


18.0 
ig 
20.8 
20.1 
21.4 
18.0 


Bay exit 


temp (°C) 


14.5 
17.0 
16.9 
17.2 
16.7 
15.2 


No. of days 
in bay 


72.89? 
711927 
83.2170 
81.587” 
106.46° 


No. of days 
in bay 


20.75 
20.06 
21.09 
27.59 
19.44 
12.27 


No. of days 
in bay 


129.00 
115.25 
86.50 
89.18 
96.47 


No. of days 
in bay 


14.71 
20.50 
18.70 
16.69 
12.00 
12.13 


et al. 2020) spent significantly more time each summer in 
the bay (over 113 d on average) than smaller females or all 
males. Therefore, the threats faced by Atlantic sturgeon in 
Chesapeake Bay are disproportionately affecting the largest 
females of the population, those most critical to the recovery 
of this population of Atlantic sturgeon. 
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Day of the year 


Sep Oct Nov’ Dec 
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Figure 2 


Proportion of all tagged Atlantic sturgeon (Acipenser oxyrinchus oxyrinchus) present 
in Chesapeake Bay during the periods before and after spawning from August 2013 
through January 2020. Presence of females and males are also shown separately. No 
females were tagged or available for detection until the fall of 2014. Passive acoustic 
telemetry was used to monitor tagged fish. 


One of the threats in Chesapeake Bay is vessel strike, 
with the tendency of Atlantic sturgeon to drift along with 
water flow (Poletto et al., 2014) making them particu- 
larly vulnerable to being pulled through the propellers of 
vessels (Brown and Murphy, 2010; Balazik et al., 2012b; 
Demetras et al., 2020). In modeling efforts, this suction of 
propellers is great enough to pull a whale up from a 
depth equivalent to 2 full ship drafts and into contact 
with propellers (Silber et al., 2010). When fish are pulled 
up to propellers, there is a direct linear relationship 
between length and probability of being struck such that 
smaller fish have a better chance than longer fish of mov- 
ing through without contacting a blade. We recommend 
that vessel operation and movement be timed to avoid, to 
the extent possible, the periods from 1 April through 
31 August and from 15 October through 1 December. 
Although this timing would avoid interactions with most 
tagged adult individuals, adult presence between 2013 
and 2019, at its extremes, was documented between 
27 February and 24 January; therefore, Atlantic stur- 
geon will face threats in Chesapeake Bay even if those 
protection periods are established. 

Adult Atlantic sturgeon use Chesapeake Bay season- 
ally, as has been reported for this species in other sys- 
tems (Dovel and Berggren, 1983; Smith, 1985; Bain, 1997; 
Dunton et al., 2010; Ingram and Peterson, 2016; Ingram 
et al., 2019). Atlantic sturgeon overwinter in the Atlantic 
Ocean, and they move into Chesapeake Bay as waters 


begin to warm. Males and females generally do not arrive 
until water temperatures are above 6°C. In general, both 
males and females have been arriving to the bay earlier 
and earlier each year (Table 2, Fig. 2). This trend of earlier 
arrival times is consistent with temperature data (from 
buoys of the Chesapeake Bay Interpretive Buoy System) 
that indicate that the warming of Chesapeake Bay begins 
earlier each year, affecting migratory timing and species 
distributions in the bay and other areas along the Atlantic 
coast of the United States (Richardson et al., 2018; Crear 
et al., 2020; Pinsky et al., 2020). Our observations support 
the conclusion of Ingram et al. (2019) that use of offshore 
habitat was influenced primarily by inshore drivers. 
Departure of tagged fish from Chesapeake Bay and back 
into the Atlantic Ocean occurred relatively quickly after 
spawning (Table 4). On average, females took approxi- 
mately 2 weeks and males took approximately 3 weeks to 
move through the bay to the ocean. Females generally left 
the York River after spawning, and males remained until all 
females had left. Because the timing of the spawning season 
of Atlantic sturgeon is driven by temperature and photo- 
period (Hager et al., 2020), male out-migration into the bay 
was significantly later in 2016 and 2018, and the late move- 
ment from river to bay became even more extreme in 2019 
(Table 2). Despite the change in the day that males reentered 
the bay from rivers, there was no corresponding change in 
the day males left the bay for the ocean. Given that nearly 
every male spawns every year (Hager et al., 2020), acoustic 
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Table 3 


Range of ordinal dates of occupancy in Chesapeake Bay 
for each size class of Atlantic sturgeon (Acipenser oxyrin- 
chus oxyrinchus) monitored with acoustic telemetry from 
August 2013 through January 2020, by sex and season, 
before (pre) and after (post) spawning. The size classes are 
below or above average size, which was 1518 and 1880 mm 
in fork length [FL] for males and females, respectively. The 
year 2016 was a leap year. A day of the year exceeds 365 or 
366 for a particular year because at least 1 fish from that 
size class remained in Chesapeake Bay into the following 
calendar year. NA=not available because no females of 
that size class were detected in that year. 


Above average 
(mm F'L) 


Below average 
(mm FL) 


Males 


276-327 
106-224 
280-327 
101-242 
279-333 
69-241 
283-392 
71-245 
262-369 
91-246 
287-335 
85-239 
284-351 


275-312 
91-225 
271-326 
93-239 
271-370 
91-240 
285-359 
67-243 
276-356 
91-243 
279-374 
63-230 
285-373 


Above average 
(mm FL) 


Below average 
(mm FL) 


Season Females 

Post NA NA 
Pre NA NA 
Post 293-326 278 
Pre 85-200 103-244 
Post 278-349 273-287 
Pre 87-247 88-242 
Post 275-330 284-309 
Pre 68-210 90-205 
Post 270-347 260-320 
Pre 59-243 87—235 
Post 277-335 282-374 
Pre 58-252 64—280 
Post 272-365 270-389 


detections of tagged males in Chesapeake Bay in the fall 
were generally a result of them leaving the York River. How- 
ever, because fewer than half of all females spawn each year, 
detections of tagged females from September until they left 
the bay for the ocean were not necessarily a result of them 
leaving a river. Atlantic sturgeon continued to leave the bay 
at similar temperatures each year regardless of spawning 
status, size, or sex (Tables 2 and 4), indicating that the cue 
for exiting the bay is either day or temperature and not 


simply a rapid migration from spawning river to ocean once 
spawning is complete. 


Conclusions 


Atlantic sturgeon face the greatest risk of interacting with 
anthropogenic threats when in estuarine and fresh water. 
Virginia and Maryland both support a variety of fisheries 
in Chesapeake Bay that have bycatch of Atlantic sturgeon 
(Welsh et al., 2002; Stein et al., 2004; ASSRT, 2007; Trice’). 
Water quality in Chesapeake Bay has improved somewhat 
in the past 50 years, but large areas are still hypoxic or 
anoxic (Scavia et al., 2021; Tian et al., 2022). Vessel strikes 
are the most frequently observed source of mortality of 
Atlantic sturgeon in Chesapeake Bay (in records of Atlan- 
tic sturgeon salvaged and reported under National Marine 
Fisheries Service permit no. 21858), despite very low 
reporting rates (Brown and Murphy, 2010; Balazik et al., 
2012b; Fox et al.*). 

The results of our research indicate that Atlantic stur- 
geon use Chesapeake Bay primarily when water tempera- 
tures exceed 6°C. When Chesapeake Bay is cooler, Atlantic 
sturgeon in coastal waters off mid-Atlantic states seek 
thermal refuge in offshore waters (Ingram et al., 2019; 
Rothermel et al., 2020; Hager and Breault*). The amount 
of the year during which water temperatures are higher 
than 6°C in Chesapeake Bay is increasing relatively rap- 
idly. When this study began in 20138, the period of offshore 
overwintering by Atlantic sturgeon was 2—4 months, but 
the length of that period has fallen to approximately 35 
d in 2019. Waters along the Atlantic coast of the United 
States are expected to continue warming through the 21st 
century (IPCC, 2023), affecting the timing of arrivals and 
departures of Atlantic sturgeon into and from Chesapeake 
Bay. As waters in Chesapeake Bay provide more days of 
hospitable temperatures for Atlantic sturgeon, we antici- 
pate and have observed a corresponding increase in their 
use of Chesapeake Bay, particularly by the largest females 
(>1880 mm FL), who currently spend significantly longer 
in Chesapeake Bay than fish of any other sizes (Table 4). 
All adult Atlantic sturgeon in Chesapeake Bay are sub- 
ject to the threats of vessel strike, poor water quality, and 
bycatch. Vessel strikes appear to pose the greatest threat 


2 Trice, G. E., IV. 2014. Observing striped bass and Atlantic stur- 
geon bycatch in a striped bass fishery using raised footlines in 
the Chesapeake Bay, 9 p. Fishery Resource Grant FRG 2014-02. 
Va. Inst. Mar. Sci., William Mary, Gloucester Point, Va. [Avail- 
able from website.] 

3 Fox, D. A., E. A. Hale, and J. A. Sweka. 2020. Examination of 
Atlantic sturgeon vessel strikes in the Delaware River estuary, 
36 p. Final report for NMFS award no. NAI6NMF4720357. 
[Available from Agric. Annex Rm. 123, Del. State Univ., 1200 
N. DuPont Hwy., Dover, DE 19901.] 

* Hager, C., and D. K. Breault. 2023. Spatiotemporal distribu- 
tions of species detected within Virginia’s offshore lease areas. 
Volume 2: The Virginia Wind Energy Lease Area A-0483, 141 p. 
Report to Bureau of Ocean Energy Management for contract no. 
140M0122P0023. [Available from Chesapeake Scientific LLC, 
100 Sixpence Ct., Williamsburg, VA 23185.] 
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Table 4 


Comparison of average seasonal residency in Chesapeake Bay and its rivers for tagged 
Atlantic sturgeon (Acipenser oxyrinchus oxyrinchus) monitored with acoustic telemetry 
from August 2013 through January 2020, by sex and size class. Average days of arrival 
in the bay and then in a river, days of departure from that river and then from the 
bay, temperatures in the bay and river at entry and exit, and number of days tagged 
fish spent in the bay each season, before and after spawning, are provided. The size 
classes are below or above average size (BA or AA), which was 1518 and 1880 mm in 
fork length for males and females, respectively. Letters denote statistically significant 
differences (P<0.05), with values having like letters not being significantly different 
from values in each column for arrivals and departures separately. Days of entry and 
exit are given as days of the year. 


Before spawning 


River entry No. of 
temp (°C) days 


Bay entry 
temp (°C) 


Bay entry River 
Sex day entry day 


hi3.61" 


218.88” OOD Palo. 
179.35? 11.46 23.18 83.36" 
196.09° 13.78 26.10° 86.15? 
191.04°° 13518 24.82? 77.45" 


1e4ew Te” 
100.512 
Male 111.56” 
112.18” 


Female 


After spawning 


Bay exit 
temp (°C) 


River exit 
temp (°C) 


River exit Bay exit 
day day 


294.33 311.31 20.25 16.69 
298.07 313.44 18.86 16.01 
Agi Ao 316.04 18.96 15.08 
298.60 316.92 18.46 15.42 


to populations of this species, and the larger the fish, the 
less likely it is to pass through a propeller without being 
struck, disproportionately threatening large females 
who also have the greatest fecundity. The U.S. Congress 
established the Endangered Species Act to require federal 
agencies and encourage states to protect endangered spe- 
cies like the Atlantic sturgeon. The data provided herein 
should allow managers to establish appropriate avoidance 
periods to protect Atlantic sturgeon and prevent the extir- 
pation of populations of Chesapeake Bay. 


Resumen 


El esturién del Atlantico (Acipenser oxyrinchus oxyrin- 
chus) esta incluido la Ley de Especies Amenazadas desde 
2012, con una fraccion de la poblacién de la bahia de Ches- 
apeake listada como en peligro. Rastreamos los lapsos de 
ocupacion y movimiento de esturiones adultos en la bahia 
de Chesapeake del Atlantico procedentes del rio York para 
identificar cuando es probable que el esturién del Atlan- 
tico interactue con las amenazas antropogénicas. Moni- 
toreamos 84 esturiones del Atlantico adultos (40 machos 
y 44 hembras) desde agosto de 2013 hasta enero de 2020 
mediante telemetria actstica. Tanto los peces desovantes 
como los no-desovantes utilizaron regularmente la bahia de 


Chesapeake, con hembras y machos llegando desde el 27 
de febrero (cuando la temperatura media del agua fue de 
7.7°C) y el 4 de marzo (6.5°C) y partiendo el 24 de enero 
(6.3°C) y el 27 de enero (6.5°C), respectivamente. La mayor 
ocupacion en la bahia por el esturién del Atlantico fue del 
1 de abril al 31 de agosto y de nuevo del 15 de octubre al 
1 de diciembre. Antes de desovar, las hembras de tamano 
mayor al promedio (>1880 mm de longitud furcal) pasaron 
significativamente mas tiempo en la bahia (>113 d) que las 
hembras mas pequenas y todos los machos. Por lo tanto, 
las hembras capaces de producir la mayor cantidad de 
huevos estuvieron desproporcionadamente expuestas a las 
amenazas antropogénicas. Aunque los cambios en las fechas 
de llegada y salida no son estadisticamente significativos, 
durante los 7 anos de este estudio, tanto los machos como 
las hembras generalmente llegaron antes y se marcharon 
mas tarde cada ano que el anterior, aumentando en un mes 
la residencia en la Bahia de Chesapeake. 
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Abstract—Somatic growth parameters 
are used in age-structured stock assess- 
ment models, such as those used to 
assess the federally managed northern 
anchovy (Engraulis mordax); therefore, 
incorrect estimation of growth can cause 
errors that affect estimates from stock 
assessments. To our knowledge, we com- 
pleted the first comprehensive investi- 
gation to model somatic and otolith 
growth of the central subpopulation of 
northern anchovy (CSNA), which has a 
range from Northern California to Baja 
California, Mexico, by using fishery- 
dependent and fishery-independent 
data. Five growth models were fitted 
to length-at-age data, including 2 mod- 
els that account for seasonal growth 
oscillations, to determine the model 
that best described growth. Seasonal 
growth oscillations were found for the 
CSNA, and the best fitting models were 
the ones that accounted for seasonal- 
ity. The fastest growth occurred during 
summer, and growth decreased 90% in 
winter. Large variations in length at age 
and extensive overlap of lengths across 
ages were found for the CSNA, typical 
observations for fish with an opportu- 
nistic life history strategy that respond 
rapidly to changing environmental 
conditions in dynamic ecosystems. Tra- 
ditional growth models overestimated 
mean asymptotic length and underesti- 
mated the growth coefficient in compar- 
ison to results from the seasonal growth 
models. Accurately describing growth of 
the CSNA by accounting for seasonality 
is important to limiting biases that may 
propagate in stock assessments and 
management decisions. 
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Somatic growth rate is an import- 
ant determinant of population growth 
and recruitment in marine fish spe- 
cies (Kerns and Lombardi-Carlson, 
2017). Clupeoids have an opportunistic 
life history strategy, which is char- 
acterized by rapid somatic growth, 
high mortality rates, short life spans, 
and protracted spawning seasons in 
response to environmental conditions 
(Blaxter and Hunter, 1982; Armstrong 
and Shelton, 1990; Winemiller and 
Rose, 1992). Clupeoids opportunisti- 
cally utilize favorable environmental 
conditions that vary on different scales, 
such as the conditions of the Califor- 
nia Current Ecosystem (CCE), which 
vary on daily, interannual, and decadal 
scales (McClatchie, 2014), leading to 
high variation in growth rates of indi- 
viduals and large annual and seasonal 
fluctuations in growth (Spratt, 1975; 
Mallicoate and Parrish, 1981; Cubillos 
et al., 2001). The rapid response of fish 
to favorable environmental conditions 
leads to changes in growth rate and 


mortality that contribute to natural, 
large fluctuations in recruitment and 
population sizes of clupeoids (Blaxter 
and Hunter, 1982; Thayer et al., 2017; 
Sydeman et al., 2020). 

Estimating somatic growth rate 
correctly is important, as it is a pri- 
mary biological characteristic used in 
fisheries assessments, management, 
and research (Kerns and Lombardi- 
Carlson, 2017). Somatic growth char- 
acteristics are used as parameters in 
age-structured stock assessment mod- 
els, such as those used to assess bio- 
mass of coastal pelagic species (CPS) 
in the North Pacific Ocean along 
the West Coast of the United States 
(Methot, 1989; Crone et al., 2019; 
Kuriyama et al., 2020, 2022). Incor- 
rectly estimating growth characteris- 
tics, whether as a result of aging errors 
or use of improper growth models, can 
cause errors that affect estimates of 
other biological characteristics, such 
as age at maturity, length at age, 
and weight at age, and estimates of 
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fishery indices, such as catch at age and biomass, that 
depend on age data (Reeves, 2003). Incorrect estimation 
of growth characteristics can also lead to errors in esti- 
mates of recruitment and total allowable catch allocated 
to fisheries in assessment models (Reeves, 2003) and to 
errors that can significantly mask relationships between 
stocks and recruitment and between environmental fac- 
tors and year-class strength (Fournier and Archibald, 
1982; Richards et al., 1992). 

Many fishes have seasonal oscillations in growth, and 
not incorporating this seasonality in data analysis when 
it is present can lead to incorrect estimation of somatic 
growth characteristics (Liu et al., 2021). Misspecifica- 
tions of growth models can affect the computation of 
recruitment (age-0 fish) and estimation of biomass (age- 
1+ fish), information that is used to manage CPS stocks 
and set harvest guidelines (Reeves, 2003; Stawitz et al., 
2019; Kuriyama et al., 2020). Unless seasonal effects can 
be explicitly accounted for in growth models, only infor- 
mation for samples collected from the same season should 
be compared (Shoup and Michaletz, 2017), an important 
consideration when samples are collected year-round. 
When seasonality was ignored for small species of com- 
mercial importance in coastal areas of the China Seas, 
the fat greenling (Hexagrammos otakii) and yellow 
croaker (Larmichthys polyactis), Liu et al. (2021) found 
that biological reference points were underestimated, 
an outcome that was attributed to incorrect estimation 


of growth parameters resulting from ignoring seasonal - 


growth patterns. The underestimation of biological ref- 
erence points indicates that a species is less productive 
than it actually is and can result in the loss of potential 
yield (Liu et al., 2021). 

The most commonly used function to model fish growth 
is the von Bertalanffy growth function (VBGF) (von Berta- 
lanffy, 1938; Ogle et al., 2017), which does not account for 
seasonality. Several models that modify the VBGF have 
been proposed to account for seasonal growth oscillations 
(Ogle et al., 2017). Seasonal growth models, like the 2 com- 
monly used functions of Somers (1988) and Pauly et al. 
(1992), account for periods of faster and slower growth 
within a year. The Somers (1988) function allows modeled 
growth to be reduced, completely stopped, or decreased, 
whereas the function of Pauly et al. (1992) does not allow 
decreasing growth but can be used to model prolonged 
periods of no growth. 

The northern anchovy (Engraulis mordax) is one of 
the most abundant forage fish species in the CCE and 
is consumed by many different species (e.g., marine 
mammals, seabirds, salmon, tunas, and sharks), mak- 
ing it an important part of the food web (Koehn et al., 
2017; Thompson et al., 2022). Similar to other clupeoid 
species, the northern anchovy is characterized by high 
interannual variability in recruitment success and pop- 
ulation abundance in response to environmental condi- 
tions (Methot, 1983; MacCall et al., 2016; Thayer et al., 
2017; Sydeman et al., 2020). Additionally, the northern 
anchovy is a multiple batch spawner with a protracted 
spawning season; spawning may occur year-round, but 


peak spawning generally occurs from January through 
April (Hunter and Macewicz, 1980; Hunter and Leong, 
1981; Dorval et al., 2018). 

Three subpopulations, the northern, central, and south- 
ern subpopulations, are recognized for northern anchovy 
on the basis of morphometric and meristic characteristics 
(McHugh, 1951; Vrooman et al., 1981); however, Lecomte 
et al. (2004) did not detect genetic structure across its range. 
The central subpopulation of northern anchovy (CSNA) 
ranges between San Francisco, California, and the middle 
of the Baja California Peninsula in Mexico. The CSNA is 
managed in U.S. waters under the Coastal Pelagic Species 
Management Plan (PFMC, 2019) and was last assessed in 
2021 (Kuriyama et al., 2022), 26 years after the previous 
assessment. The current assumption is that the overall 
population size is at or around its highest level since the 
development of quantitative monitoring in the CCE around 
70 years ago (Thompson et al., 2019, 2022). 

No recent work has included examination of somatic 
growth of the CSNA, despite its ecological role in the CCE 
and the importance of growth estimates in informing stock 
assessments. The most recent work examining growth of 
northern anchovy through the adult stage was done by 
Parrish et al. (1985), and a few limited studies have pro- 
vided values of growth parameters (i.e., mean asymptotic 
length [Z,,] and the growth coefficient [K]) calculated by 
using traditional growth models (Spratt, 1975; Cisneros 
et al., 1990). In no work have seasonal growth oscillations 
of northern anchovy been modeled by using traditional 
growth models (e.g., the VBGF) with data on fish collected 
year-round, but seasonal variations in growth of northern 
anchovy have been generally observed (Mallicoate and 
Parrish, 1981; Methot, 1983; Parrish et al., 1985; Butler, 
1989; Takahashi et al., 2012). Seasonal growth oscilla- 
tions have been reported for other engraulids (Palomares 
et al., 1987; Cubillos and Arancibia, 1993; Bellido et al., 
2000; Bilgin et al., 2013; Cerna and Plaza, 2016). North- 
ern anchovy reportedly live up to 8 years (Mallicoate and 
Parrish, 1981), but older individuals (j.e., age-4+ fish) 
are thought to be relatively scarce (Collins, 1969; Spratt, 
1975; Mallicoate and Parrish, 1981; Parrish et al., 1985). 
Understanding the age structure and growth patterns of 
northern anchovy is important, as it is a model species 
for examining and projecting the effects of climate change 
on clupeoids and marine communities in the CCE and 
on fisheries management (Muhling et al., 2019, 2020; 
Tommasi et al., 2021). 

The goal of this study was to update the age and growth 
patterns of the CSNA to inform future stock assessments. 
The objectives were to model growth and determine 
whether the CSNA has seasonal growth oscillations 1) by 
comparing estimates from 3 growth models, the Gompertz, 
logistic, and von Bertalanffy growth functions, to integer 
age data; 2) by calculating fractional ages and compar- 
ing estimates between the VBGF and 2 models allowing 
seasonal growth oscillations (Somers, 1988; Pauly et al., 
1992); and 3) by examining seasonal variations in the rela- 
tionships of weight to body length and of body length to 
otolith length. 
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Materials and methods 
Data sources and collection methods 


Samples from the CSNA were collected during both fishery- 
dependent and fishery-independent surveys from 2015 
through 2021. Fishery-independent data were collected 
as part of acoustic-trawl surveys of CPS conducted aboard 
the NOAA Ships Reuben Lasker and Bell M. Shimada by 
the NOAA Southwest Fisheries Science Center (SWFSC) 
in summer and fall (June—September) for each year in 
2015-2021, as well as in spring (March and April) in 2017 
and 2021 (Dorval et al., 2018; Stierhoff et al., 2019; 
Zwolinski et al., 2019; Stierhoff et al., 2020, 2021a, 2021b). 
The primary objectives of the SWFSC trawl survey are 1) to 
survey the distributions and abundances of CPS and their 
abiotic environments in the CCE and 2) to measure size— 
frequency distributions and to gather life history informa- 
tion (Stierhoff et al., 2020; Renfree et al., 2022), in order to 
produce length and age composition data for use in acoustic 
biomass estimation and stock assessment models. 

Sagittal otolith samples were collected by using a 2-stage 
sampling design that follows best practices (Miranda and 
Colvin, 2017). For the first stage, a maximum of 75 indi- 
viduals of the CSNA were randomly collected from each 
nighttime trawl haul, and each individual fish was mea- 
sured for standard length (SL) to the nearest millimeter 
and for total weight (TW) to the nearest 0.5 g. For the 
second stage, 25 individuals that had lengths that encom- 
passed the size distribution of the first-stage samples were 
selected for otolith extraction and aging by the SWFSC. A 
complete description of the data collection methods of the 
SWFESC trawl survey can be found in Dorval et al. (2022) 
and Schwartzkopf et al. (2022). 

Fishery-dependent samples were collected by the Califor- 
nia Department of Fish and Wildlife (CDFW) during port 
sampling of commercial fishery landings. The CDFW sam- 
pling plan for CPS is based on a stratified-random design, 
where samples are collected during 12 random days within 
each month of the year for each of the port areas where 
the majority of landings occur: the areas of Monterey, Santa 
Barbara, and Los Angeles. The majority of fish were collected 
off Monterey, but in 2016 some additional samples were col- 
lected off Los Angeles. For each vessel from which landings 
were sampled, 25 fish were randomly selected throughout 
the offload period and brought back to a lab, where they 
were measured for SL to the nearest millimeter and for TW 
to the nearest 0.1 g. Otoliths were then extracted for aging 
to be completed by the CDF'W. A complete description of the 
design and methods of data collection by the CDFW can 
be found in Dorval et al. (2022). The combined data from 
fishery-dependent and fishery-independent surveys used 
in the CSNA stock assessment models were collected year- 
round, allowing seasonality to be explored. 


Age determination 


Whole (i.e., not sectioned) otoliths were immersed in dis- 
tilled water in a glass dish with a black background for 
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no longer than 3 min, counts were made of the number of 
annuli observed on the distal side of the otolith, and edge 
type was determined. Aging was conducted by using a dis- 
secting microscope with reflected light at a magnification 
of 25x. An annulus is defined as the interface between an 
inner translucent growth increment and the successive 
outer opaque growth increment (Collins and Spratt, 1969; 
Yaremko, 1996; Schwartzkopf et al., 2022). To our knowl- 
edge, annual age determination has not been validated 
for the northern anchovy, but first annulus formation and 
annual age determination has been validated for the Euro- 
pean anchovy (E. encrasicolus) (Aldanondo et al., 2016; 
Uriarte et al., 2016; Basilone et al., 2020); therefore, we 
assumed this validation held true for the northern anchovy. 

Integer age was assigned on the basis of a June birthdate 
by using the capture date and interpretation of the most dis- 
tal pair of increments (Collins and Spratt, 1969; Schwartz- 
kopf et al., 2022). Because samples were collected year-round, 
fractional age was also determined by dividing the integer 
age by 12 and by using a June birthdate. Age data from both 
fishery-independent and fishery-dependent surveys and 
from all years were combined because we were interested 
in modeling growth of the CSNA as a whole (similar to the 
modeling done in the stock assessment) and because effects 
of these factors (i.e., year and type of survey) are beyond the 
scope of this work. Sample sizes by age for each month of 
sampling can be found in Supplementary Table 1 (online only). 


Relationships of standard length to weight and to 
otolith length 


We fit a generalized linear model (GLM) incorporating 
season for all individuals measured for length and weight 
(sample size [n]=19,209) to examine the weight—length 
relationship of the CSNA, similar to the work done by 
Palance et al. (2019). The GLM was fitted by using the 
base function glm in R (vers. 4.2.1; R Core Team, 2022) and 
has the following structure: 


(TW) = By + B Season; + Balog, (SL,,;), (1) 


where / = the log link function; 
TW; =the total weight (in grams) of fish z in season j; 
Season = a categorical factor with fall, winter, spring, and 
summer; and 
SL; ; = the standard length (in millimeters) of fish 7 in 
season /. 


Note that TW;; has a gamma distribution, with the mean 
calculated as follows: 


TW: op alo* Br Season; +B» log. (SZ;,;)} (2) 
ye 


The variance of the gamma distribution was calculated as 
follows: 


POTENT (3) 


where v = a dispersion parameter estimated during fitting. 
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The season was determined on the basis of the month of 
capture, with months of capture from September through 
November assigned to the fall season, the months from 
December through February considered to be winter, the 
months from March through May assigned to spring, and 
the months from June through August considered to be 
summer. The model fit was determined with the adjusted 
amount of deviance accounted for by a GLM, D? (Guisan 
and Zimmermann, 2000), calculated by using the Dsquared 
function in the R package modEvA, vers. 3.8.4 (Barbosa 
et al., 2013). The effect of season among all season pairs was 
compared by running the GLM 4 times, each time with a dif- 
ferent season as the intercept, and evaluating the t-statistic 
for every season effect (e.g., if fall was the intercept, the 
t-statistics for the individual effects of winter, spring, and 
summer were evaluated relative to the fall intercept). 

We chose a random subset of otoliths (1=20386), from fish 
collected during both SWFSC trawl surveys and CDFW port 
sampling that had lengths that encompassed the size distri- 
bution within each season, to examine the relationship of SL 
to otolith length and to confirm that otolith growth reflects 
somatic growth. The right-side otolith was chosen when 
possible and was imaged by using a dissecting microscope 
with an AmScope camera (United Scope LLC, Irvine, CA) and 
analyzed with the R package shapeR (vers. 0.1-5; Libungan 
and Palsson, 2015) to calculate otolith length (Feret length) at 
the time of capture. We fitted a linear model that has season 
as a variable and has the following structure: 


(SL; ;] = B, + B, Season; + Bo (OL, ;) 
+ Bs (Season,OL, ;), (4) 


where OL ; ; = the otolith length Gin millimeters) of fish 7 
in season /. 


A model with an interaction between season and OL; ; 
was chosen, given that previous work on CPS found reduced 
effects of seasonality on otolith length as the length of fish 
increased (Thomas, 1983). A model with no interaction term 
was also fitted, but the model with the interaction term had 
a lower Akaike information criterion (AIC) than the model 
with no interaction term (AIC: 138632 and 13800, respec- 
tively) and was deemed a better fit (Akaike, 1973; Burnham 
and Anderson, 2002; Burnham et al., 2011). The effect of 
season among all season pairs was compared by running 
the linear model 4 times, each time with a different season 
as the intercept, and evaluating the t-statistic for every sea- 
son effect (e.g., if fall was the intercept, the ¢-statistics for 
the individual effects of winter, spring, and summer were 
evaluated relative to the fall intercept). Visual examination 
of the diagnostic plots of both of the models used to examine 
the relationships of weight to length and of length to oto- 
lith length revealed that the assumptions of normality and 
homogeneity of variance were met. 


Somatic growth modeling 


We fitted 5 growth models to length-at-age data to exam- 
ine which model best describes somatic growth of the 


CSNA. The Gompertz growth model (Gompertz, 1825) is 
as follows (Ricker, 1975; Ogle et al., 2017): 


-s"(t-1"] 
-e 


Eile (5) 


where L, = the mean length at age f; 
L,, = the mean asymptotic length (i.e., the maximum 
mean length for the population); 
g =the instantaneous growth rate at the inflection 
point ¢ ; and 
t = the maximum instantaneous growth rate. 


The logistic growth model (Campana and Jones, 1992; 
Ogle et al., 2017) is as follows: 


L, = L, ; £ sal (6) 


where g_,, = the positive instantaneous growth rate as age 
approaches negative infinity. 


The final 3 models are parameterizations or modifica- 
tions of the VBGF (von Bertalanffy, 1938), which has the 
following common foundation: 


L, =L,(1-e), (7) 


where g = at least a function of t. 


For the VBGF attributable to Beverton and Holt (1957), 
q takes this form: 


q=K(t-ty), (8) 


where K = the growth coefficient, a measure of the expo- 
nential rate at which L, approaches L,.; and 
ty = the theoretical age at which L; is zero. 


As CPS have seasonal variations in their somatic 
growth rates, 2 modifications to the VBGF were fitted 
to model seasonal oscillations in growth, similar to the 
adjustments to the VBGF done by Ogle (2017). The first 
modification of the VBGF is from Hoenig and Choudaray 
Hanumara (1982) and Somers (1988), with a clarification 
by Garcia-Berthou et al. (2012). In this model, hereafter 
referred to as the Somers model, q takes this form (Ogle, 
2017): 


q= K(t = ty + <= sin (2n(¢ “= t,)) 


-* sin (2n(t =e )) (9) 


where ¢, = the inflection point, or the amount of time 
between time zero and the start of the convex 
portion of the first sinusoidal growth oscilla- 
tion; and 
C =the proportional decrease in growth at the 
depth of the growth oscillation, normally the 
“winter” season. 
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If C=0 (i.e., if there is no seasonal oscillation in mean 
length), the equation reduces to the VBGF. If C>0 and 
C<1, the increase in mean length is reduced but does not 
stop; if C=1, the increase in mean length comes to a com- 
plete stop; and if C>1, mean length decreases during “win- 
ter.” The time of the year when growth is slowest, known 
as the winter point (WP), was calculated as follows: 


WP =t, + 0.5. (10) 


Pauly et al. (1992) argued that a C greater than 1 is bio- 
logically improbable for fish vertebrates whose skeletons 
largely preclude shrinkage in length and, therefore, pro- 
posed a modification to the Somers model in which mean 
length was not allowed to decrease. In this model, hereaf- 
ter referred to as the Pauly model, q takes this form: 


ae PL Gia 2 F 

oo aoe | on inl al =) 
K’\1- NGT 

side hal -+)} (11) 


where NGT = the no growth time, the length of the no 
growth period as a fraction of a year. 


The age ¢’ is found by “subtracting from the real age (¢) the 
total no-growth time occurring up to age zt” (Pauly et al., 
1992; Ogle, 2017). Pauly et al. (1992) suggested using K’ 
to minimize confusion with K in the VBGF as the units 
changed from year’ to (I-NGT)*. 

Model parameters and 95% confidence intervals were 
estimated by using nonlinear least squares regression 
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with the R package FSA (vers. 0.9.3; Ogle et al., 2022) 
and the nlsBoot function in the R package nlstools (vers. 
2.0-0; Baty et al., 2015), respectively. Starting values were 
obtained by visually fitting the models to the observed 
data, and alternative starting values were used to confirm 
a global minimum was obtained. The best fitting model 
was chosen by using the AIC values (Akaike, 1973; Ogle, 
2017), with the best fitting model having the lowest AIC 
and other models ranked according to the differences in 
their AIC values from that of the best model (Burnham 
and Anderson, 2002; Burnham et al., 2011). 


Results 
Length and age structure 


A total of 19,209 individuals were measured for length 
and weight, and 12,520 otoliths from fish collected during 
both fishery-dependent and fishery-independent surveys 
were aged from 2015 to 2021 and used in subsequent 
analyses. Lengths from all measured fish ranged from 23 
to 164 mm SL, lengths from aged fish ranged from 44 to 
164 mm SL, and ages ranged from 0 to 6 years (Fig. 1, A 
and B). Age-0 and age-1 fish accounted for 62% of the sam- 
ples, and a steady decline in the frequency of older aged 
fish was observed, with age-4+ fish accounting for 4% of 
the samples (Fig. 1B). All fish <69 mm SL were age 0. 
The GLM used to examine the relationship of weight to 
length was significant (P<0.0001) and fit the data well 
(D*=0.95). Winter did not differ from spring (t=0.36, 
P=0.71), but all other seasons were significantly different 
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Figure 1 


(A) Histogram of the standard lengths of all northern anchovy (Engraulis mordax) from 
the central subpopulation collected and measured for length (sample size [n]=19,209) 
and (B) histogram of the age distribution of sampled fish (n=12,520). This subpopu- 
lation ranges from Northern California to Baja California, Mexico. Data used in both 
histograms are from northern anchovy collected in 2015-2021 during sampling of ports 
by the California Department of Fish and Wildlife and acoustic-trawl surveys of coastal 
pelagic species conducted by the NOAA Southwest Fisheries Science Center. 
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from each other (P<0.0001), including fall from winter 
(t=4.11), spring (t=5.12), and summer (¢=36.58) and winter 
from summer (t=22.57) and spring (t=-31.30) (Table 1). 
Looking at the model predictions, fish from summer had a 
larger weight at a given length (Fig. 2A). Although fish 
from winter and spring statistically differed from fish 
from fall, the model predictions look similar among these 
3 seasons (Fig. 2A). Full results from the GLM can be 
found in Supplementary Table 2 (online only). 


Table 1 


Equations used to determine the relationship of total 
weight (TW) to standard length (SL) in the general lin- 
ear model for each season for the central subpopulation of 
northern anchovy (Engraulis mordax). This subpopulation 
ranges from Northern California to Baja California, Mexico. 
Parameter estimates for each season are given within the 
equation. For example, the estimates from the fall model 
are as follows: -—12.085 for the intercept and 3.122 for the 
slope. The model was fitted to a combination of fishery- 
independent and fishery-dependent data from 2015 through 
2021, and the sample size (7) for each season is provided. 


Season Model equation 

TW = e(-12.085)SL?1”2 
TW = e(-12.071)SL? 7 
TW = e(-12.073)SL?1 
PW = 1 2001SE 


Fall 
Winter 
Spring 
Summer 
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A total of 2036 otoliths were imaged to calculate oto- 
lith length. The linear model of the relationship of SL to 
otolith length was significant (P<0.05) and fit the data 
well (coefficient of multiple determination [R*]=0.90). 
Fall, winter, and spring did not differ from each other, 
including fall from winter (t=1.48, P=0.15) and spring 
(t=-0.31, P=0.75) and winter from spring (¢=-1.58, 
P=0.11), but summer differed from all the other seasons 
(from fall: t=9.74; from winter: t=8.99; and from spring: 
t=7.77; P<0.0001) (Table 2). Looking at the model pre- 
dictions, otolith length was greatest in winter, followed 
by fall, spring, and summer for small fish (less than 
~120 mm SL) (Fig. 2B). As fish length increased, the 
effects of seasonality were reduced, and otolith length 
was similar between seasons for fish of the same size 
(Fig. 2B). Full results from the linear model can be found 
in Supplementary Table 3 (online only). 


Somatic growth modeling 


Large variation in length at age was found among indi- 
viduals from the CSNA, with higher variation in length 
at age for young fish (<2 years) and smaller variation for 
old fish (>5 years) (Fig. 3, A and B). The von Bertalanffy, 
Gompertz, and logistic growth functions yielded asymp- 
totic growth patterns when fitted to both integer and 
fractional age data (Fig. 3, A and B). Curves for modeled 
length at age are nearly identical for all 3 models at all 
ages (Fig. 3, A and B), with similar parameter estimates 
within each age data set (Table 3). Notably, parameter 


Otolith length (mm) 


Season -= Fall -- Winter -- Spring ~~ Summer 


Figure 2 


Relationships (A) between total weight and standard length and (B) between standard length and otolith length by 
season for the central subpopulation of northern anchovy (Engraulis mordax) that ranges from Northern California 
to Baja California, Mexico. The lines represent the predictions from the growth models used in this study, and the 
points represent the raw data for fish collected during fishery-independent and fishery-dependent surveys from 2015 


through 2021. 
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Table 2 


Equations used to determine the relationship of standard 
length (SL) to otolith length (OL) in the linear model for each 
season for the central subpopulation of northern anchovy 
(Engraulis mordax). This subpopulation ranges from North- 
ern California to Baja California, Mexico. Parameter esti- 
mates for each season are given within the equation. For 
example, the estimates from the fall model are as follows: 
—2.316 for the intercept and 32.134 for the slope. The model 
was fitted to a combination of fishery-independent and 
fishery-dependent data from 2015 through 2021, and the 
sample size (n) for each season is provided. 

Season Model equation 
Fall 
Winter 
Spring 
Summer 


SL = -2.316 + 32.134(OL) 
SL = -8.584 + 33.460(OL) 
SL = 0.042 + 31.864(OL) 

SL = 24.491 + 26.243(0L) 


estimates differed between the integer age and frac- 
tional age data sets for the same model, with estimates 
of L,, being greater and estimates of the growth coeffi- 
cients (i.e., K, g , g_..) being lower with a given model 
fitted to the integer age data set than with that model 
fitted to the fractional age data set (Table 3). For 


Standard length (mm) 


Age (years) 


Model 


Standard length (mm) 
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example, with the VBGH, L,, was estimated at 130.6 mm 
SL and K was estimated at 0.52 year’’ when the inte- 
ger age data set was used, whereas L.,, was estimated 
at 125.1 mm SL and K was estimated at 0.72 year’ 
when the fractional age data set was used. The VBGF 
was the best fitting model for both the integer age and 
fractional age data sets, followed by the Gompertz and 
logistic models (Table 4). 

When modeling seasonal oscillations of growth, the 
Somers model and Pauly model were indistinguishable 
(Fig. 4, Table 5) and fit the data better than the VBGF, 
although the Somers model had the lowest AIC value 
and was deemed the best fit. With the Somers model, L,, 
was estimated at 126.2 mm SL and K was estimated at 
0.68 year’. Growth was fastest during the first year of 
growth. There was evidence of a 90% reduction in growth 
during the winter season (C=0.9), and the month of slow- 
est estimated growth was December (WP=0.62). We fit a 
VBGEF to integer age data from August and September to 
compare parameter outputs when seasonality was 
removed. Data from August and September were chosen 
because data for these months included information for 
individuals with the full range of ages in sufficient num- 
bers (Suppl. Table 1) (online only). When seasonality was 
removed from this model fitted with integer age data, L,. 
was estimated at 126.9 mm SL and K was estimated at 
0.78 year’, values similar to the parameter estimates 
from the Somers model. 
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Figure 3 


Growth curves from von Bertalanffy, Gompertz, and logistic growth functions fitted to (A) length-at-integer-age 
data and (B) length-at-fractional-age data for the central subpopulation of northern anchovy (Engraulis mordax) 
that ranges from Northern California to Baja California, Mexico. The solid lines indicate model predictions based 
on observed data, and the shaded areas indicate the 95% confidence intervals for these predictions. The black circles 
(with transparency) indicate the raw data points for individuals. The data used in the models are for fish collected 
from 2015 through 2021 during fishery-independent and fishery-dependent surveys. Parameter estimates from the 


model fits are provided in Table 3. 
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Table 3 


Equations used in 3 growth models to estimate the mean length at age ¢ (L,) for 
the central subpopulation of northern anchovy (Engraulis mordax) that ranges 
from Northern California to Baja California, Mexico. Parameter estimates from 
each model are given within the equation; for example, estimates from the first 
equation for the von Bertalanffy growth function are as follows: 130.6 mm in 
standard length for the mean asymptotic length, 0.52 year! for the growth 
coefficient, and 2.51 years for the theoretical age at which L, is zero. The von 
Bertalanffy, Gompertz, and logistic growth functions were fitted to both integer 
and fractional age data for fish collected from 2015 through 2021 during fishery- 
independent and fishery-dependent surveys. 
Model df 


Data set Model equation 


Integer ages von Bertalanffy 12,517 


L, = 130.6 E .: ae 


ra 
Je 0.60(t+1.95) 
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Gompertz P2514 Ee 120 5e=70e0 


ot 
Logistic 12,517 L, = 128.6| 1 A speed 


Fractional ages von Bertalanffy 12 


Gompertz 
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Table 4 


Statistical measures of fit for von Bertalanffy, Gompertz, 
and logistic growth functions fitted to data sets of lengths 
at integer and fractional ages for northern anchovy 
(Engraulis mordax) of the central subpopulation col- 
lected during 2015-2021. This subpopulation ranges from 
Northern California to Baja California, Mexico. For each 
model, number of parameters (K), the Akaike informa- 
tion criterion (AIC) and the difference in its AIC value 
from that of the best model (the one with the lowest AIC) 
(AAIC) are provided. 


Data set Model K AIC 

97,208.86 
OF 22 1°82 
97,235.67 
96,428.34 
96,477.62 
96,526.79 


von Bertalanffy 
Gompertz 
Logistic 
Fractional ages von Bertalanffy 
Gompertz 
Logistic 


Integer ages 


Discussion 


The results of this study provide an update of the under- 
standing of CSNA growth patterns, given that the last 
work examining growth of the CSNA through the adult 
stage was completed 38 years ago (i.e., Parrish et al., 
1985). We found patterns in somatic growth of the CSNA 
that are typical of r-selected and opportunistic species, 
which have small body sizes, low maximum ages, and 
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high growth rates (Adams, 1980; Winemiller and Rose, 
1992; King and McFarlane, 2003). Additionally, fast 
growth during the first year of life, high individual varia- 
tion in growth rates, and seasonal fluctuations in growth, 
as found in our study, are characteristic of CPS in upwell- 
ing environments (Cubillos et al., 2001; Takahashi et al., 
2012). High variability in growth rates among individuals 
from the CSNA has been postulated to lead to increased 
population growth rates due to size-dependent mortal- 
ity and, therefore, to lead to increased population size 
(Lo et al., 1995). The growth characteristics of the CSNA 
contribute to high interannual variability in recruitment 
success and population abundance observed for this pop- 
ulation of northern anchovy (Lasker, 1981; MacCall et al., 
2016; Thayer et al., 2017). 

Large variations in length within each age class, exten- 
sive overlap of lengths across ages, rapid somatic growth 
during the first year of life, and high proportions of young 
fish (<4 years old) with scarcity of older individuals have 
been found in studies investigating the CSNA (Collins and 
Spratt, 1969; Spratt, 1975; Mais, 1981; Mallicoate and 
Parrish, 1981; Parrish et al., 1985) and other species of 
Engraulis (Bacha et al., 2010; Cerna and Plaza, 2016). 
Spratt (1975) modeled growth of individuals from the 
CSNA collected during both fishery-dependent and 
fishery-independent surveys and found a higher L,, 
(165 mm SL) and lower K (0.29 year”) than the estimates 
for those parameters from our study, and Cisneros et al. 
(1990) found an L., of 153 mm SL and a K of 0.7 year? for 
northern anchovy collected from the Gulf of California in 
Mexico (the southern subpopulation). The higher L,, and 
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Figure 4 


Growth curves from the von Bertalanffy growth function, Somers model, and 
Pauly model fitted to length-at-fractional-age data for the central subpopu- 
lation of northern anchovy (Engraulis mordax) that ranges from Northern 
California to Baja California, Mexico. The solid lines indicate model predic- 
tions based on observed data, and the shaded areas indicate the 95% confi- 
dence intervals for these predictions. The black circles (with transparency) 
indicate the raw data points for individuals. The data used in the models are 
for fish collected from 2015 through 2021 during fishery-independent and 
fishery-dependent surveys. Parameter estimates and their 95% confidence 


intervals from the model fits are provided in Table 5. 


lower K found by Spratt (1975) reflect their sample distri- 
bution, with lengths of sampled fish ranging from 88 to 
160 mm SL and no age-0 fish. Without large and old fish, 
L, is underestimated and K is overestimated (Bolser et al., 
2018); therefore, it follows that the absence of small and 
young fish would result in overestimation of L,, and under- 
estimation of K, as is revealed when comparing estimates 
from Spratt (1975) and those from our study. Cisneros 
et al. (1990) collected fish with size and age ranges similar 
to the fish collected for our study (~50—-140 mm SL, 0-5 
years), likely contributing to similar K values. However, 
they modeled growth using length—frequency analysis 
instead of ages calculated from otolith readings, possibly 
driving the differences in L,.. 

Seasonal CSNA growth oscillations and large varia- 
tions of length within each age class were not unexpected 
and reflect the life history characteristics of this species in 
response to environmental conditions of the CCE. Somatic 
growth of the CSNA was most rapid during the summer 
months in our study, with increased length at a given 
weight in summer compared to that in other seasons, a 
result that also has been found for northern anchovy col- 
lected in the CCE (Palance et al., 2019). This rapid sum- 
mer growth is likely in response to warmer sea-surface 
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temperatures and the increase in pro- 
ductivity from upwelled, nutrient-rich 
waters that stimulates growth of phy- 
toplankton (Checkley and Barth, 2009; 
Jacox et al., 2018). Although productiv- 
ity is increased during summer months, 
food aggregations are patchy and may 
be composed of prey of different levels 
of quality (Lasker and Zweifel, 1978; 
Litz et al., 2010), potentially resulting 
in large, opportunistic feeding events 
of variable nutritional quality for some 
individuals of the CSNA. Differences in 
the quality of prey encountered between 
individuals from the CSNA can contrib- 
ute to high variability in their growth 
rates. _ 

Additionally, the protracted spawning 
season of the CSNA may also contrib- 
ute to large variations of length within 
each age class, especially at younger 
ages (Methot, 1983). For example, an 
age-0 fish born early in the year would 
be larger in length than an age-0 fish 
born later in the year when collected 
at the same time. The use of fractional 
ages in our work helped to address this 
phenomenon, but the birthdates of indi- 
vidual fish are still unknown. The lack 
of birthdates may bias any estimate 
of actual somatic growth; therefore, 
caution should be taken when making 
interpretations. 

Somatic growth of the CSNA was 
slowest in winter, a result that is linked 
to cool sea-surface temperatures and reproductive strat- 
egies (Cubillos and Arancibia, 1993; Cubillos et al., 2002; 
Bilgin et al., 2013). The peak CSNA spawning season is 
from January through April (Hunter and Leong, 1981; 
Reiss et al., 2008); therefore, the slowest growth in Decem- 
ber aligns with a shift from somatic growth to gametic 
growth in preparation for spawning. The duration of 
spawning season increases with female length, weight, 
and age, and although a large majority of individuals are 
mature by age 0, older females contribute disproportion- 
ately more to total annual egg production than first-time 
spawners (Hunter and Macewicz, 1985; Parrish et al., 
1986; Schwartzkopf et al., 2022). The reduction in season- 
ality in somatic growth of older individuals is likely due to 
the increased investment in reproduction over more of the 
year and reduced somatic growth at older ages (>3 years) 
(Armstrong and Shelton, 1990; Kozlowski, 1996). If fish 
mature at different ages, individual variability in repro- 
ductive investment may also contribute to high variability 
in growth rates among individuals (Brosset et al., 2016; 
Claireaux et al., 2021). 

Large variations and extensive overlaps in length at 
age present challenges in estimating other life history 
characteristics and consequently in conducting stock 
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Table 5 


Estimates and 95% confidence intervals for parameters of the von Bertalanffy growth function (VBGF) fitted to integer age data, 
the VBGF fitted to fractional age data, the Somers model, and the Pauly model used to examine the growth of the central sub- 
population of northern anchovy (Engraulis mordax) that ranges from Northern California to Baja California, Mexico. The model 
parameters are mean asymptotic length (Z..); growth coefficients (K and Kk’); the theoretical age at which length is zero (¢,), the 
inflection point, or the amount of time between time zero and the start of the convex portion of the first sinusoidal growth oscilla- 
tion (¢,); the proportional decrease in growth at the depth of the growth oscillation (C); and no growth time (NGT), or the length of 
the period of no growth as a fraction of a year. For each model, the Akaike information criterion (AIC) and the difference in its AIC 
value from that of the best model (the one with the lowest AIC) (AAIC) are provided. No AIC or AAIC values are given for the VBGF 
fitted to integer age data as that data set is different from the data set of fractional ages used with the other models. All models 
were fitted to length-at-age data for fish collected during fishery-independent and fishery-dependent surveys from 2015 through 
2021. SL=standard length. 


Parameter estimates (95% confidence intervals) 


K 


Model L..(mmSL) K(year) ((4-NGT)} t, (year) NGT (year) AIC AAIC 


ty (years) 


VBGF 130.6 0.52 —2.51 
(integer) (129.1-132.3) (0.47—0.56) (-2.68——2.35) 
VBGF 125.1 0.73 -1.53 
(fractional) (124.2-125.9) (0.69-0.77) (-1.63—-1.45) 
Somers 126.2 0.68 —1.74 0.12 0.90 
model (125.4—-127.1) (0.64—0.71) (-1.79--1.67) (0.10-0.15) (0.77—1.00) 
Pauly model 126.4 — 0.67 -1.76 0.13 — 0.00 


96,428.3 172.10 
96,256.2 0.00 


96,257.9 1:69 


(125.5-127.4) (0.63—0.73) (-1.81--1.61) (0.11-0.15) (0.00—0.08) 


assessments. When time or costs do not allow all samples 
that were measured for length to be aged, age—length keys 
(ALKs) are used to draw inferences about the age compo- 
sition of a stock and are used in stock assessments to esti- 
mate growth or mortality rates (Coggins et al., 2013) and 
to estimate abundance and catch at age (Kimura, 1977). 
When ages and growth characteristics are estimated with 
extrapolation-based methods (e.g., ALKs), more errors are 
generated for populations with extensive overlap of length 
distributions across ages and large variability in length at 
age (Miranda and Colvin, 2017). These errors could prop- 
agate throughout an assessment when an ALK is used to 
calculate weight at age, as has been done in the CSNA 
stock assessment (Kuriyama et al., 2022), which is ulti- 
mately used to calculate recruitment and biomass quanti- 
ties relevant for management of this population. 
Although the VBGF is the most typical model used to 
describe growth of CPS on the U.S. Pacific coast (Arm- 
strong and Shelton, 1990; Crone et al., 2019; Kuriyama 
et al., 2020, 2022), no other work has examined seasonal 
modifications to the VBGF for northern anchovy. We found 
that the best fitting growth model for the CSNA was the 
Somers model, which allows growth to oscillate seasonally. 
Seasonal growth oscillations have been reported for small 
CPS around the world: the European anchovy (E. encra- 
sicolus) (Bellido et al., 2000; Bilgin et al., 2013), Peruvian 
anchoveta (E. ringens) (Palomares et al., 1987; Cubillos 
and Arancibia, 1993; Cerna and Plaza, 2016), European 
pilchard (Sardina pilchardus) (Silva et al., 2008), and 
Araucanian herring (Strangomera bentincki) (Cubillos 
et al., 2002; Feltrim and Ernst, 2010). However, in many 


of the studies that generated these reports, length-based 
methods were used to estimate growth parameters with 
seasonal growth models, possibly yielding less reliable 
estimations than those from the use of age-based methods 
(Liu et al., 2021). 

We found that the VBGF overestimated L,, and under- 
estimated K in comparison to estimation with the Somers 
model, and Liu et al. (2021) had the same observation 
when seasonal growth oscillations of small species were 
ignored. The misestimation of growth parameters when 
seasonality is ignored can propagate biases in stock 
assessments because underestimating K results in under- 
estimation of natural mortality and possibly in charac- 
terization of a stock as less productive than it actually is 
(Scherrer et al., 2021). Additionally, when seasonality was 
included with the estimation of fractional ages, overlaps in 
length at age were reduced, and that reduction may result 
in more accurate ALKs and, therefore, reduce bias in stock 
assessments. 

Validating that otolith growth has a relationship with 
somatic growth is imperative, as it permits the use of 
data from otolith studies (e.g., increment widths) to test 
for relationships between fish growth and climate indices 
or regional environmental conditions, and information 
about such relationships is important for projections of 
climate- and ecosystem-level change on marine communi- 
ties (Muhling et al., 2018; Erisman et al., 2021; Tommasi 
et al., 2021). The linear relationship between otolith 
length and fish length indicates that otolith growth does 
reflect somatic growth to a degree for the CSNA; how- 
ever, it is not uncommon for otolith and somatic growth 
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to be decoupled because somatic growth is tied to the 
environment and otolith growth is tied to metabolic rate 
(Buckmeier et al., 2017). 

The differences in otolith growth among seasons 
observed in our study may also provide examples of the 
uncoupling of otolith and somatic growth reported for 
some other fish species that have seasonally different 
growth rates (e.g., Francis et al., 1993). In the summer 
season, when the fastest growth was found, mean otolith 
length for fish of a given length was smaller than that for 
fish in other seasons (and growth rates were lower), with 
seasonal differences diminished at larger fish lengths 
(greater than ~120 mm SL). The larger otolith lengths in 
winter may also indicate that otoliths continue to grow 
in seasons during which somatic growth is reduced, a 
common phenomenon (Secor and Dean, 1992; Buckmeier 
et al., 2017). The negative relationship between otolith 
size and somatic growth rate is predicted for individuals 
erowing according to the VBGF (Xiao, 1996), but the pres- 
ence of a strong growth effect can create substantial errors 
for studies in which otoliths are used to back-calculate fish 
length (Campana, 1990). Although we did not directly test 
the effect of growth rates on the relationship between oto- 
lith length and fish length, the results of our study provide 
a starting point for examination of this relationship in the 
future. 

One potential source of uncertainty in this study is the 
sample distribution. Only 162 fish (0.008%) in our study 
were <50 mm SL, possibly a reflection of the selectivity 
of the commercial fishery, size selectivity of trawl nets, 
locality of trawl tows, or a combination of these factors. 
The commercial fishery is known to avoid small fish (gen- 
erally <90 mm SL) because small individuals may clog 
nets (Huppert et al., 1980). More importantly, the fishery 
is extremely limited by market demand because large 
anchovies have a higher fat content and are, therefore, 
more desirable for their main uses as feed for the finfish 
mariculture industry, pet food, oil for pharmaceuticals, 
live bait, and fresh or canned food for human consumption 
(CDFW, 2022; Kuriyama et al., 2022). Therefore, fishery- 
dependent data will consistently have this source of uncer- 
tainty, and fishery-independent data can help fill this gap. 

The net used in SWFSC acoustic-trawl surveys has 
an 8-mm mesh liner to retain juvenile and adult fish as 
small as 50 mm SL and potentially as large as 800 mm 
SL (Dotson and Griffith, 1996; Dotson et al., 2010; Dorval 
et al., 2022), but small fish are known to be poorly rep- 
resented in catches made with nets with similar sizes of 
mesh liner in contrast to their representation in purse 
seine samples (Dotson and Griffith, 1996). Small juveniles 
of the CSNA tend to remain closer to shore (Parrish et al., 
1985), but NOAA research vessels are unable to operate 
in waters with depths less than 20-30 m. Therefore, any 
CPS aggregations in shallow, nearshore areas are missed; 
it should be noted that the SWFSC is addressing this 
data gap through the inclusion of sampling of nearshore 
areas by commercial fishing vessels with the sampling 
conducted during its fishery-independent surveys (Dorval 
et al., 2022; Renfree et al., 2022). Except for 1 individual, 
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all fish <50 mm SL were caught in 2015, the year in 
which abundances of young of the year in the CSNA were 
found to have increased drastically (Zwolinski et al., 2017; 
Thompson et al., 2019). The presence of small individuals 
in 2015 indicates that small fish (<50 mm SL) will be col- 
lected when they are present in high numbers with the 
current sampling design and trawl! net. 

The disproportionate distribution of samples from both 
the fishery-independent and fishery-dependent surveys, 
with a decline in older age classes (>4 years), may also be 
due to the sampling designs used. The commercial fishery 
fleet fishes anchovy close to shore, where younger fish are 
common (Parrish et al., 1985), and close to ports because 
fish do not survive well in the holds over the time it takes 
to cover long distances from shore (Kuriyama et al., 2022). 
The SWESC trawl survey is designed to sample the entire 
populations of CPS to assess their biomass (Stierhoff et al., 
2020), a different objective than that of a survey focused 
purely on age and growth in which ideally all age classes 
are sampled uniformly regardless of relative abundance 
(Bolser et al., 2018). 

It is important to ensure that growth data accurately 
represent the variability of the lengths at each age, as sam- 
pling that is not representative of the underlying popula- 
tion may produce biased estimates of growth characteristics 
and mortality (Sainsbury, 1980; Goodyear, 2019). Further 
sampling of older age classes could cause the current esti- 
mates of growth characteristics to change, but attempting 
to add samples in size bins for which catch is more rare in 
size-stratified sampling can substantially increase bias in 
mean size estimates (Goodyear, 2019). The current 2-stage 
sampling design of the SWFSC trawl survey meets both of 
these important needs through minimization of bias due 
to sampling distribution by first selecting samples at ran- 
dom so that they are representative of what is caught and 
then by making sure a sufficient number of otoliths from all 
length (and age) classes of a species or population are aged 
to accurately approximate the true variation in length at 
age (Schwartzkopf et al., 2022). 

The high percentage of fish from age 0 to age 2 in sam- 
pling and the low percentage of fish at age 4+ reflect the 
current belief that the maximum age of northern anchovy 
has declined over time with catches being heavily dom- 
inated by age-O and age-1 fish and having diminished 
levels of older age groups (Mais, 1981; MacCall, 2009). 
This result is in line with a high mortality rate typical of 
clupeoids (Armstrong and Shelton, 1990) and the CSNA 
(Kuriyama et al., 2022), assuming that the survey samples 
are representative of the adult population. Given the low 
harvest rate by the commercial fishery (PFMC, 2022), it is 
reasonable to conclude that the high mortality rates reflect 
natural mortality (e.g., mortality that results from preda- 
tion or starvation). Additionally, true length at age may 
be underestimated because slower-growing marine fish 
can have higher predation mortality rates (Fennie et al., 
2020). The northern subpopulation of northern anchovy is 
thought to have larger, older, and faster-growing fish than 
those in the CSNA (Parrish et al., 1985; Stierhoff et al., 
2020; senior author, unpubl. data). Future work could 
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examine growth characteristics for the entire species with 
the addition of older individuals (>4 years old) as well as 
could compare characteristic estimates between subpopu- 
lations and stratify by latitude. 


Conclusions 


This work is the first comprehensive investigation into 
modeling non-seasonal and seasonal growth oscillations 
in northern anchovy by using an age-based method. The 
results indicate the importance of incorporating seasonal- 
ity when modeling somatic growth of CPS with data from 
samples collected year-round. If seasonality is ignored, 
L,, may be overestimated and growth coefficients (i.e., K) 
may be underestimated. The similarity in parameter esti- 
mates from the VBGF fitted to integer age data when 
only 2 months of data were included and from the Somers 
model fitted to fractional age data when data from year- 
round sampling were included indicates that 1) the sea- 
sonal model appears to account well for seasonal growth 
oscillations of the CSNA and 2) a VBGF fitted with inte- 
ger ages may be appropriate for use with data collected 
over a short time period. Accurately describing growth is 
important to decrease biases that may propagate in stock 
assessments and subsequently in management decisions. 


Resumen 


Los pardmetros de crecimiento somatico se utilizan en 
modelos de evaluaciédn de poblaciones estructurados 
por edades, como los que se emplean para evaluar y mane- 
jar federalmente la anchoa del norte (Engraulis mordax); 
por lo tanto, un estimado incorrecto del crecimiento puede 
causar errores que afecten a los estimados de las evalua- 
ciones de poblaciones. Hasta donde sabemos, hemos com- 
pletado la primera investigacion exhaustiva para modelar 
el crecimiento somatico y de los otolitos de la subpoblaci6n 
central de la anchoa del norte (CSNA), que tiene su distri- 
bucidn desde el norte de California hasta Baja California, 
México, utilizando datos dependientes e independientes de 
la pesca. Para determinar el modelo que mejor describia el 
crecimiento, se ajustaron cinco modelos de crecimiento a 
los datos de talla por edad, incluidos 2 modelos que tienen 
en cuenta las oscilaciones estacionales del crecimiento. Se 
observaron oscilaciones estacionales del crecimiento en el 
CSNA, y los modelos que mejor se ajustaban fueron los que 
tenian en cuenta la estacion. El crecimiento mas rapido se 
produjo durante el verano, y disminuy6 90% en invierno. 
Se observaron grandes variaciones en la longitud a la edad 
y un amplio solapamiento de longitudes entre edades en 
el CSNA, observaciones tipicas de peces con un modo de 
vida oportunista que responden rapidamente a las condi- 
ciones ambientales cambiantes en ecosistemas dinamicos. 
Los modelos de crecimiento tradicionales sobrestimaron 
la longitud asint6tica media y subestimaron el coeficiente 
de crecimiento en comparacién con los resultados de los 
modelos de crecimiento estacional. Describir con precisién 


el crecimiento de la CSNA teniendo en cuenta la estacio- 
nalidad es importante para limitar los sesgos que pueden 
propagarse en las evaluaciones de las poblaciones y en las 
decisiones de manejo. 
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Abstract—In this study, we examined 
whether otolith microconstituents are 
deposited seasonally in a manner sim- 
ilar to optical annulus formation and 
thereby can be used to validate age 
interpretations. In temperate species, 
seasonal temperature changes drive 
the formation of optical annuli, and 
we hypothesized that they similarly 
caused oscillations in microconstituent 
deposition. Using laser ablation induc- 
tively coupled plasma mass spectrom- 
etry analysis, we tested for periodicity 
in the deposition of barium (Ba), cal- 
cium, copper, magnesium (Mg), manga- 
nese (Mn), phosphorus (P), strontium, 
and zinc in otoliths and compared that 
periodicity to the periodicity of the 
annulus zonation (optical opaque and 
translucent zones). For this investiga- 
tion of chemical annulus periodicity, 
we used black sea bass (Centropristis 
striata), a species with yearly optical 
annulus formation that has been val- 
idated. Periodicities in elemental pro- 
files in otoliths from black sea bass 
were detected by using Lomb—Scargle 
periodogram analyses. Optical annu- 
lus formation aligned with Mg profiles, 
but periodicity in annular deposition of 
other elements—Ba, Mn, and P—was 
observed only after the first optical 
annulus, an outcome that is indicative 
of ontogenetic changes in habitat (from 
coastal to near-shelf waters) known to 
occur during the juvenile life stage of 
this species. Certain elements, such 
as Mg, identified through this otolith 
analysis and the periodogram analyses 
could be applied to species for which no 
validated aging procedure exists. 
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The lack of validated methods for age 
determination directly affects many 
fisheries assessments and manage- 
ment efforts. The use of unvalidated 
aging techniques can affect estimation 
of catch at age, selectivity at age, num- 
bers at age, maturity at age, and weight 
at age, and errors in estimates of these 
factors can contribute to major bias and 
errors in age-based assessments (Brad- 
ford, 1991; Wilhelm et al., 2008; Elvarsson 
et al., 2018). An error in aging can result 
in overexploitation and the eventual col- 
lapse of a stock and its associated fishery 
(Mills and Beamish, 1980; Chale-Matsau 
et al., 2001; Yule et al., 2008). Despite the 
importance of accurate age determina- 
tions, validation studies on annulus 
formation are rare for many import- 
ant commercial stocks, including the 
majority of species harvested in waters 
along the Atlantic coast from Maine to 
South Carolina (Suppl. Table) (online 
only). Validation studies have been done 
for only slightly more than half (10 of 19; 
Suppl. Table) (online only) of the species 
routinely aged by the NOAA Northeast 
Fisheries Science Center. 


3 Jackson School of Geosciences 
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Interpretation of optical zones 
(annuli) in fish hard parts is consid- 
ered a direct method of aging and is the 
conventional approach for age determi- 
nation, which is needed in age-based 
assessments (Casselman, 1983; Vitale 
et al., 2019). Hard parts are calcified 
structures, which include otoliths, verte- 
brae, opercula, fin rays, and scales. They 
often have patterns of seasonal zonation, 
which can be optically distinguished 
through microscopy as alternating dark 
and light (translucent and opaque) zones 
that vary in their organic concentrations 
(Jolivet et al., 2008) (Fig. 1). Opaque 
zones typically have higher protein con- 
tent, are relatively narrow in comparison 
to translucent zones, and are interpreted 
as representing periods of lower sea- 
sonal growth rates such as temperate 
winters (Thomas, 1983; Machias et al., 
1998; Zlokovitz et al., 2003). The 2 zones 
together constitute an annulus, but it 
remains the task of the investigator to 
determine whether each annulus forms 
annually and therefore whether it is 
valid to use the structure for determin- 
ing ages of fish (Campana, 2001). 
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Figure 1 


Age annotations for each optical annulus identified on a (A) sectioned and (B) whole otolith from 
a black sea bass (Centropristis striata) caught on 15 July 2018 off Ocean City, Maryland. (C) Oto- 
lith sections, like the one shown in the bottom panel, were annotated following analysis through 
laser ablation inductively coupled plasma mass spectrometry. Distances were measured from the 
ablated point within the core region of an otolith to the distal edge of each opaque zone found in 
optical annulus assignments. Images were captured with reflected light on a black background. In 
the top panels, the asterisk indicates the central opaque core, which was not counted according to 
the otolith aging approach used by the NOAA Northeast Fisheries Science Center. 


Efforts to validate age determination can be hampered 
by the requirement to test age interpretations on fish of 
known ages. The approaches used to validate the annual 
formation of presumed annuli in hard parts can be direct or 
indirect (Vitale et al., 2019). In direct validation methods, 
individuals are tracked over time as annular zones form 
over seasons and years, by using temporal reference points 
(markers) on hard parts (Vitale et al., 2019). In indirect val- 
idation methods, annular structures are examined in rela- 
tion to changes in season, year-class strength, or size modes 
and include length—frequency analysis, marginal incre- 
ment analysis, and microchemistry (Vitale et al., 2019). 

In a promising approach for age validation, periodicity 
in annular deposition of microconstituents in hard parts is 
examined by using microchemistry. Otolith microconstit- 
uent uptake can vary seasonally and can be regarded as 
a “chemical calendar-clock” with time-keeping properties 
similar to the optical annuli on hard parts (Limburg et 
al., 2018). The results of the study by Hiissy et al. (2021a) 
reveal that the elemental composition of otoliths can 
vary seasonally and that these variations can be used 


to validate the accuracy of otolith chemical chronologies 
in determining the age of fish. Siskey et al. (2016) and 
Brophy’, employing periodogram analysis, have shown 
that the periodicity of some elements in hard parts is 
annual and synchronizes with the optical patterns that 
form annuli. 

The underlying premise of microconstituent validation 
is that optical banding patterns (annuli) in hard parts coin- 
cide with changes in chemical composition, both of which 
represent seasonal changes in otolith organic content and 
environmental temperature. Fish that reside or migrate 
through coastal and shelf environments are exposed 
to mixtures of terrestrial, anthropogenic, and marine 
sources of aqueous metals and isotopes. Found in higher 
abundance in hard parts than calcium (Ca), the elements 
strontium (Sr), magnesium (Mg), manganese (Mn), and 


1 Brophy, D. (coord.). 2019. Validating age-determination [sic] 


of anglerfish and hake, 97 p. EASME/EMFF/2016/1.3.2.7/ 
S12.762036. Final report. Publ. Off. Eur. Union, Luxembourg, 
Luxembourg. [Available from website.] 
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barium (Ba) substitute for Ca within the crystalline lat- 
tice of aragonitic otoliths in proportion to environmental 
exposure but at rates modified by temperature and growth 
(Hiissy et al., 2021b). Trace constituents, such as phospho- 
rus (P), cobalt (Co), copper (Cu), and zinc (Zn), are bound 
in the protein matrix of otoliths and it has been hypoth- 
esized that they are physiologically controlled (Htissy et 
al., 2021b), albeit some evidence indicates a role for envi- 
ronmental exposure (Arslan and Secor, 2005). Therefore, 
across the elements that we considered candidates for use 
in age determination, the kinetics of tracer incorporation 
are driven variously by exposure, temperature, and physi- 
ology, the relative contribution of which remains an unset- 
tled question (e.g., Hussy et al., 2020). 

In temperate shelf environments, concentrations of can- 


didate elements do not vary as substantially as they would 


in freshwater, estuarine, and nearshore environments. In 
this study, we assumed that exposure to seasonally chang- 
ing temperatures drives variation in microconstituent 
concentrations. We tested this assumption by comparing 
known seasonal patterns of optical annulus zonation with 
those of microconstituent zonation in otoliths. 

To evaluate chemical annulus periodicity, it’s important 
to test a model species for which age determination has 
been validated previously. The black sea bass (Centropris- 
tis striata) is a ubiquitous reef fish species found in the 
Gulf of Mexico and along the Atlantic coast of the United 
States, from Florida to the Gulf of Maine. This species is 
a target for important commercial and recreational fish- 
eries along the Atlantic coast that are assessed by using 
statistical catch-at-age models (Musick and Mercer, 1977; 
NEFSC, 2020; ASMFC?). Adult black sea bass are most 
abundant in temperate shelf environments, although 
their young (<2 years in age) often occur in estuarine 
and coastal habitats. The aging method for black sea bass 
has been validated through marginal increment analysis 
(Dery and Mayo, 1988; Penttila and Dery, 1988; Robillard 
et al.®: Koob et al., 2021). 

We hypothesized that cycles of microconstituent depo- 
sition (i.e., chemical zonation) would relate to age with 
a dominant periodicity of 1 cycle per optical annulus. By 
linking the optical interpretation to the chemical mea- 
sures, the periodicity of each cycle can be tested through 
time-series analysis. We first identified candidate micro- 
constituents through elemental profiles and 2-dimensional 
(2D) maps. We then conducted Lomb-—Scargle periodogram 
analyses on candidate microconstituent profiles to test the 
hypothesis that the dominant periodicity is 1.0 cycle per 
optical annulus. Under the premise that black sea bass 


2 ASMFC (Atlantic States Marine Fisheries Commission). 2021. 
Review of the interstate fishery management plan for black sea 
bass (Centropristis striata): 2020 fishing year, 15 p. ASMFC, 
Arlington, VA. [Available from website.] 

3 Robillard, E., J. W. Gregg, J. Dayton, and J. Gartland. 2016. 
Validation of black sea bass, Centropristis striata, ages using 
oxytetracycline marking and scale margin increments, 17 p. 
Stock Assess. Rev. Comm., SARC 62 working paper. [Available 
from Northeast Fish. Sci. Cent., Natl. Mar. Fish. Serv., NOAA, 
166 Water St., Woods Hole, MA.] 
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undertake a habitat shift from inshore to offshore waters 
during their second year of life, 2 alternate periodogram 
analyses were conducted with data sets containing con- 
centrations of microconstituents along a line from the core 
to the edge of otoliths for 1) the entire series of annuli on 
each otolith or 2) the series of annuli after the first one. 


Materials and methods 
Analysis of otoliths 


Sagittal otoliths were collected from black sea bass 
(sample size [n]=125) ranging in total length from 19.5 
to 36.3 cm and captured through hook and line at wreck 
and artificial reef sites (depths: 20-30 m) off Ocean City, 
Maryland, during the summer and fall of 2016 and 2018. 
Annuli in whole and sectioned otoliths were defined as 
continuous opaque zones with no breaks (Fig. 1). Aging 
conventions were adopted from Robillard et al.® and Koob 
et al. (2021), who discounted initial opaque material 
deposited around the core of an otolith as material formed 
during age 0 and identified a second opaque region as the 
completion of the first annulus (Fig. 1) (Penttila and Dery, 
1988; Robillard et al.°). | 

By convention, a birth date of 1 January was assumed 
(given that spawning occurs in the spring for this spe- 
cies) (Vanderkooy et al.*). The seasonality of optical and 
chemical annuli could vary; therefore, fish collected before 
30 June were assumed to have completed less than half of 
their most recent year of life, and counts of opaque zones 
were assumed to represent the estimated age. For those 
captured after 30 June, the edge was counted as an integer 
because these fish had completed >0.5 years of life and 
so that we could consider the indication of growth pat- 
terns beyond the last opaque zone (Dery and Mayo, 1988; 
Vanderkooy et al.*). Whole otoliths from black sea bass 
were photographed under reflected light at 1.5x magnifi- 
cation by using a microscope camera and Infinity Analyze 
and Capture” software (vers. 6.5.6; Teledyne Lumenera, 
Ottawa, Canada). Following age assignment with whole 
otoliths, otoliths were then embedded in Struers epoxy 
resin (Struers A/S, Ballerup, Denmark) for sectioning. 
A 2.0-mm-thick section containing the core region (the 
innermost portion of the otolith, pertaining to the first 
year of life) was cut along a transverse plane with an 
IsoMet Low Speed Saw (Beuhler Ltd., Lake Bluff, IL). The 
otolith section was then attached to a glass slide by using 
Crystalbond 509 mounting adhesive (Aremco Products 
Inc., Valley Cottage, NY) in preparation for microchemical 
analysis and further interpretations of annuli. 


* Vanderkooy, S., J. Carroll, S. Elzey, J. Gilmore, and J. Kipp. 2020. 
A practical handbook for determining the ages of Gulf of Mexico 
and Atlantic coast fishes, 3rd ed., 268 p. Gulf States Mar. Fish. 
Comm., Ocean Springs, MS. [Available from website. ] 

° Mention of trade names or commercial companies is for identi- 
fication purposes only and does not imply endorsement by the 
National Marine Fisheries Service, NOAA. 
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Microchemical analysis of elemental variations from the 
core region to the otolith’s periphery was done by using 
laser ablation inductively coupled plasma mass spectrome- 
try (LA-ICP-MS) at the Department of Geological Sciences 
of the University of Texas at Austin. Analysis was limited 
to 34 of the 125 otoliths collected because of the cost of 
the analysis; otoliths were randomly selected from among 
otoliths taken from fish assigned ages >1 year. The instru- 
mentation used was an NWR193 excimer laser ablation 
system (193-nm wavelength, 4-ns pulse width; Elemental 
Scientific Inc., Omaha, NE) coupled to an Agilent 7500ce 
ICP-MS (Agilent, Santa Clara, CA). This LA-ICP-MS sys- 
tem was equipped with a large-format, 2-volume sam- 
ple cell with fast washout (<1 s) that accommodated all 
samples and standards (USGS MACS-3, NIST 612, and 
ECRM-752-NP) in 8 laser cell loadings. 

The LA-ICP-MS system was optimized daily for sensitiv- 
ity across the atomic mass unit range and low oxide produc- 
tion (ratio of thorium monoxide to thorium: 0.58 [standard 
deviation (SD) 0.06]) by tuning with a standard glass (NIST 
612), and sensitivity was checked with trial transects on 
representative specimens. Following pre-ablation (75-ym 
spot, 50-ym/s scan rate, 3.3-J/cm? energy density [fluence]) 
to remove shallow surface contaminants, core-to-edge 
profiles were collected from each otolith along the sulcal 
ridge for consistency, by using a spot with a 50-um diam- 
eter, a 5-um/s scan rate, a fluence of 2.92 J/em? (SD 0.15), 
a 15-Hz repetition rate, and carrier gas flows of 0.8 L/min 
for argon and 0.80—0.85 L/min for helium. The quadrupole 
time-resolved method was used to measure 12 masses 
with integration times of 10 ms (“Mg, **“‘Ca, and ®°Sr), 
20 ms Meg, 21P and °°Mn), and 25.ms @°Co, Cu, Zn, 
and '°"!8®Ba). The quadrupole duty cycle of 0.3732 s corre- 
sponds to 94% of total measurement time, with a resultant 
linear sampling rate of 1.867 m/s, equivalent to 13.4 mea- 
surements captured within the footprint of either aperture. 
Measured intensities were converted to elemental concen- 
trations (parts per million) by using iolite 4 software (Paton 
et al., 2011), with *’Ca as the internal standard and a Ca 
index weight percentage of 38.3% for unknowns. 

For the otoliths, USGS MACS-3 (synthetic aragonite) was 
used as the primary calibration standard, and NIST 612 
(glass) and ECRM-752-NP (nanoparticulate pressed lime- 
stone powder pellet made from BCS-CRM 393 limestone 
from Bureau of Analysed Samples Ltd., Middlesbrough, 
UK) were used as external reference standards. The grand 
average of secondary standard recovery fractions for all ele- 
ments was typically within 5% of preferred values in the 
GeoReM database (Jochum et al., 2005). 

To create 2D maps of otolith microconstituents, areas 
were scanned by using contiguous line traverses with 
a 17-um” aperture, spaced 17 pm apart. The scan rate 
(105 pm/s) was set to provide 1 quadrupole measurement 
cycle per aperture footprint. Using the quadrupole method 
for mapping, we surveyed elements with integration times 
of 10 ms (*° “Ca and ®Sr), 20 ms (°°Mn and '°Ba), and 
25 ms (*4Mg, © Cu, and °°Zn) and with a significantly shorter 
duty cycle of 0.1612 s corresponding to 90% of total mea- 
surement time. Calibration and secondary standards were 


as described previously. External reference standard recov- 
eries were typically within 10% of certified reference values. 
Image conversions, from contiguous-line traverses along the 
otolith to 2D maps of element deposition, were performed by 
using iolite 4 software (Paton et al., 2011). 


Analysis of time series 


To determine if optical annuli coincide with the micro- 
constituent profiles, LA-ICP-MS data were aligned with 
micrographs of transverse otolith sections through image 
analysis. Optical annuli were assigned and measured for 
each cycle of a paired translucent zone and opaque zone. 
Distances along the trench formed through ablation on the 
surface of an otolith were measured from the ablated point 
within the core region to the distal edge of each opaque 
zone according to optical annulus assignments (Fig. 1C), 
with a minimum of 39 measurements per year of growth. 
Each opaque zone was considered the terminus of each 
annulus (1 year of growth), and zones measured through 
LA-ICP-MS within a cycle (pair of translucent and opaque 
zones) were converted to fractional ages within each annu- 
lus cycle, assuming that otolith growth was linear. 
Typical concentrations of elements over transects on 
otoliths from representative specimens were 10—10,000 
times higher than the limits of detection for most elements 
(Table 1). A rule-of-thumb ratio of signal to noise of 2 to 3 
was employed as a threshold for inclusion (mean transect 
concentration divided by minimum detection limit) 
(Armbruster and Pry, 2008; Rajakovié et al., 2012). Con- 
secutive moving average filters were used to smooth 
derived elemental time series (profiles) by using a 9-point 
boxcar width (13- and 26-ym equivalent distances), and 
the profiles were converted to Z-scores, resulting in 
smooth, locally weighted signals that were free from high- 
frequency outliers. Signals were converted to distance (in 
microns) from the core on the basis of scan rate and duty 
cycle. Because signals for elements measured with 


Table 1 


Average concentrations (in parts per million) of elements 
in otoliths of black sea bass (Centropristis striata) (n=34) 
caught off Ocean City, Maryland, during the summer and 
fall of 2016 and 2018. Also provided for each element are the 
standard deviation (SD) of the average concentration and 
the ratio of signal to noise (S:N), which is the mean tran- 
sect concentration divided by the estimated mean detection 
limit among the analyzed otoliths. The elements are mag- 


nesium (Mg), phosphorus (P), manganese (Mn), cobalt (Co), 
copper (Cu), zinc (Zn), strontium (Sr), and barium (Ba). An 
asterisk (*) indicates that the average concentration of the 
element was less than twice as high as detection limits. 


Value Mo ocPy cM rCo-s GurZn Sr Ba 
Average 13.6 78.6 2.3 0.2 0.3 0.4 2148.9 10.1 
SD 3.0828) 08 0202-404 T68i4es: 3:5 
S:N 23:7 . 2.7.46 3.2°:0.8* 0.3* 45,000.0 370.0 
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multiple isotopes were identical, we used the isotopes with 
the highest abundance (e.g., “Mg and '°°Ba) for plotting 
and statistical analysis. The detrended time series for 
each individual fish and element were subjected to Lomb— 
Scargle periodogram analysis, a flexible approach that 
accounts for unequally spaced measurements across the 
annuli of multiple otoliths (Lomb, 1976; Scargle, 1982; Ruf, 
1999). Normalized power was used to evaluate significant 
frequencies. Alternate periodogram analyses were per- 
formed with data for the entire series of annuli on each 
otolith and with data for the series of annuli after the first 
one on each otolith. 


Results 
Analysis with mass spectrometry 
The elements Mg, Mn, Co, Cu, Zn, Sr, and Ba were all mea- 


sured at detectable levels though LA-ICP-MS (Table 1). 
Copper and Zn were below the signal, a noise threshold 


of 2.0—3.0 (both at a signal-to-noise ratio of 0.3), and were 
not considered further in analyses. Concentration ampli- 
tudes were within the range expected for otoliths from 
marine teleost fish: <1 ppm for transition metals Mn, Co, 
Cu, and Zn; <8 ppm for Mg, P, and Ba; and <2000 ppm for 
Sr (Campana, 1999; Hiissy et al., 2021b). 

Among measured elements, Mg had the sharpest con- 
trast across presumed chemical annuli as revealed in 2D 
maps of element deposition for the 2 analyzed otoliths 
(Fig. 2). The chemical zonation for 4 of the 6 elements ana- 
lyzed, Ba, Mg, Mn, and Sr, aligned with conventional light 
microscopy. For Ba, chemical zonation was most apparent 
for annuli after the second one, with weak contrast prior to 
the third annulus. Strong zonation of Mg was observed for 
the entire series of annuli. For Mn, concentrations were ele- 
vated in the core region for both scanned otoliths, and 
chemical zonation was faint but distinct across the series of 
annuli. Zonation of Sr was inconsistent between the 2 oto- 
liths (Fig. 2, “Sr” panel in both sets of panels), with a faint 
series of annuli in one otolith and a narrow chemical zona- 
tion for annuli after the second one in the other otolith. 


Figure 2 


Two-dimensional maps of element deposition, with a light micrograph (Raw) image, for each of the sectioned otoliths from 2 
black sea bass (Centropristis striata) caught on 15 July 2018 off Ocean City, Maryland. The elements are barium (Ba), calcium 
(Ca), magnesium (Mg), manganese (Mn), and strontium (Sr). Colors in the maps indicate the concentration of elements from 
blue (low) to red (high). Note that contrast is optimized and depicted through pseudocolor for each map, and it is not common 
across elements. White circles indicate optical annuli, and the associated numerals indicate assigned ages; the central opaque 


core is not included because it was not counted in otolith aging. 
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Zonation of Ca was homogeneous in both 2D maps as 
expected. 

Despite notable zonation in the 2D maps, oscillating pat- 
terns in standardized elemental profiles were not aligned 
across individuals nor clearly associated with age (Fig. 3). 
Note that the majority of profiles correspond to estimated 
ages <4 years, with 2 profiles continuing to age 6. Profiles 
for the transition metal Co were especially noisy, corre- 
sponding to the lack of patterning in the 2D scan. For Ba, 
Mg, and Mn, peaks and nadirs of oscillations give some 
indication of association with optical annuli, but visual 
classifications were deemed too subjective. 


Analysis with periodograms 


The results of the principal periodogram analysis of the 
entire elemental time series, inclusive of the first annu- 
lus, indicate a significant periodicity in microchemi- 
cal transects, albeit at frequencies that mostly differed 
from the expected 1 year (Fig. 4, top set of panels). Only 
annular deposition of Mg had significant periodicity of 
1 cycle per annulus. Significant periodicity in deposition 
is evident for the elements Ba and Mn but at a rate of 
~1.5 cycles per annulus, each at a lower normalized power 


Standardized concentration (ppm) 


than that for deposition of Mg (Fig. 4, top set of panels). 
Strontium had no defined peak in oscillations and, there- 
fore, no determinable primary chemical cycle. In the peri- 
odograms for oscillations of P and Co, the frequencies are 
1.5 cycles and 1 cycle per annulus, respectively, albeit at 
low normalized power (i.e., reduced significance). 

In the alternative periodogram analysis of the time 
series excluding the first annulus, significant periodic- 
ity was detected for oscillations of Mg, Ba, Mn, and P at 
the hypothesized 1 cycle per annulus at high normalized 
power (Fig. 4, bottom set of panels). Magnesium deposition 
periodicity remained at 1 cycle per annulus after trunca- 
tion of the time series. The periodogram for oscillations 
of Sr again failed to indicate a dominant period. In the 
truncated series, oscillations of Co had multiple peaks, all 
of which deviated from the expected frequency of 1 year. 


Discussion 


The hypothesis of aligned optical and chemical oscillations 
across all annuli (not excluding the first annulus) was 
confirmed for the element Mg and alignment across age 
estimates >1 year was confirmed for the microconstituents 


Age (years) 


Figure 3 


Profiles of microconstituent concentrations, smoothed by using consecutive moving median and average filters, in 5 of 
the 34 otoliths examined through laser ablation inductively coupled plasma mass spectrometry. Otoliths were taken 
from black sea bass (Centropristis striata) caught off Ocean City, Maryland, in October and July 2018. The microcon- 
stituents are barium (Ba), cobalt (Co), manganese (Mn), magnesium (Mg), phosphorus (P), and strontium (Sr). The 
lines for the individual otoliths are depicted in different colors for readability. 
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Period (years) 


Figure 4 


Lomb—Scargle periodograms based on elemental profiles for all otoliths examined through micro- 
chemical analysis (n=34) combined, with the top set of panels including data for the full series of all 
annuli on each otolith and the bottom set of panels including data for only the series of annuli after 
the first one. Otoliths were taken from black sea bass (Centropristis striata) caught off Ocean City, 
Maryland, during the summer and fall of 2016 and 2018. The horizontal dashed lines indicate the 
significance level of 0.001. The vertical dashed line indicates the hypothesized frequency of 1 cycle 
per annulus. Note the difference in scales between the top and bottom rows in each set of panels. 


Ba, Mg, P, and Mn through the Lomb-Scargle periodogram 
analyses. The periodicity we observed for deposition of Mg 
in otoliths through the analysis that included data for all 
annuli supports recent work by Koob et al. (2021) and the 
conventions used by the Northeast Fisheries Science Cen- 
ter for first annulus assignments. Because the dominant 
period was 1 year, rather than for instance 1.2 years, the 


influence of the assumption that otolith growth is linear 
on the estimated period was negligible. The 2D maps show 
that chemical zonation corresponded with optical zonation 
for all annuli in each otolith (in the case of Mg) or for the 
annuli after the first one in each otolith (in the cases of Ba, 
P, and Mn). Concentrations of other elements were below 
detection thresholds (Cu and Zn), had no zonation (Co), or 
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had patterns that were inconsistent with optical annuli 
(Sr). Patterns of Sr oscillations are weak in 2D maps and 
in standardized profiles. This lack of annual periodicity in 
Sr explains why X-ray dispersive spectroscopy of Sr depo- 
sition in otoliths was inconclusive in past work (Richards, 
2016; Frey, 2022). 

Both periodogram analysis results and 2D maps indi- 
cate a lack of consistent chemical zonation associated with 
the first optical annulus. The inclusion and removal of the 
first annulus in periodogram analysis had a large influ- 
ence on the periodicity of Sr, Ba, and Mn oscillations. The 
shifts in Sr, Ba, and Mn concentrations could be explained 
by habitat changes occurring during the juvenile life 
stage. Nursery habitats in estuarine and nearshore envi- 
ronments for black sea bass (Dery and Mayo, 1988) are 
heterogenous in ambient concentrations of Sr, Ba, and Mn. 
Strontium is conservatively mixed with and is often pos- 
itively associated with salinity (Kraus and Secor, 2004), 
observations that are consistent with the lower otolith 
concentrations of Sr deposited during the first year of life 
than the levels deposited during later years in otoliths 
analyzed in our study (Fig. 2). Barium becomes less avail- 
able in eutrophic waters of estuaries and coastal waters 
(Hiissy et al., 2021b), a finding that aligns with the lower 
otolith concentrations of Ba around the first optical annu- 
lus compared to the levels around the rest of the annuli 
beyond the first annulus in otoliths analyzed in our study. 
How Mn is associated with ambient exposure remains 
unsettled in the literature, although higher than normal 
concentrations of Mn in otoliths have been associated with 
hypoxic waters (Limburg and Casini, 2018). Early elemen- 
tal oscillations could also reflect a change in reproductive 
state for black sea bass (maturation is achieved at ages =2 
years) or early habitat transitions (Drohan et al., 2007). 
Understanding the lack of expected periodicity for age-1 
black sea bass requires a deeper consideration of otolith 
composition during the first year of life. 

Seasonal changes in otolith composition with regard to 
the elements Mg, Sr, and Mn have been found to conform to 
expectations of environmental and physiological influences 
(Grammer et al., 2017; Htissy et al., 2021a). In their scheme, 
Hiissy et al. (2021a) proposed that otolith uptake of Ba, K, 
and Sr, versus that of P, Cu, and Zn, was respectively under 
environmental versus physiological controls. Further, Mg 
and Mn uptake was influenced jointly by the environment 
and physiology (Hiissy et al., 2021a). Concentration of Mg 
in otoliths is known to be linked to growth (Grammer et al., 
2017), possibly explaining the strength of the periodogram 
analyses and the visibility of chemical annuli in the 2D 
maps created in our study. Results of recent work by Brophy 
et al. (2021) indicate that seasonal variation in deposition 
of Sr in otoliths is detectable in anglerfish (Lophius piscato- 
rius), but in our analysis, Sr oscillations were weak and no 
yearly periodic fluctuations were apparent. Concentrations 
of Ba in otoliths, like those of Sr, may reflect seasonal tem- 
perature cycles (Barnes and Gillanders, 2013; Reis-Santos 
et al., 2013; Stanley et al., 2015). 

A central premise of aligning otolith microconstituent 
cycles with annulus zone formation is that black sea bass 


are exposed to seasonal changes in temperatures. Moser 
and Shepherd (2008) placed archival tags in black sea bass 
released off coastal states from Massachusetts to Virginia. 
These and conventionally tagged black sea bass made off- 
shore migrations during winter months, cued by cooling 
temperatures in near-shelf regions. Off the coast of Mary- 
land, black sea bass migrate to inner shelf regions in spring 
and reside in bottom waters that have stable thermal con- 
ditions until fall, when turnover or seasonal mixing causes 
a rapid increase in water column temperatures (Wiernicki 
et al., 2020). These conditions are then followed by cooling 
in late fall, cueing winter migrations of black sea bass to 
offshore waters (Moser and Shepherd, 2008). Off the coast 
of Maryland, where this study’s samples originated, this 
pattern of residency in the spring and fall and offshore 
migration in winter would be associated with exposure toa 
seasonal amplitude in temperature of 12—17°C, with most 
of that change in temperature occurring during August— 
November (based on observed bottom temperatures in 
near- and mid-shelf waters off Deleware, Maryland, and 
Virginia; Rothermel et al. 2020). During a single, long-term 
(more than 3 years) mark-recapture study in which archi- 
val tags were used, Moser and Shepherd (2008) observed 
temperatures experienced by tagged fish that were con- 
sistent with this seasonal cycle and that ranged between 
6°C (in winter) and 25°C (from late summer to fall). As 
reviewed in the “Introduction” section, seasonal tempera- 
ture was expected to influence concentrations of Mg, Mn, 
Sr, and Ba—the elements that are incorporated into an 
otolith’s aragonitic structure. 

For temperate species, such as the black sea bass, 
changes in composition of an otolith are thought to relate 
to the formation of organic-rich (winter) and Ca-rich 
(spring, summer, and fall) zones in an otolith (Secor and 
Piccoli, 1996; Hitissy et al., 2021b). These optical zones are 
also associated with seasonal changes in temperature, 
with the translucent zone forming during high growth 
seasons of summer and fall and the narrow opaque zone 
forming during winter and early spring (Zlokovitz et al., 
2003). For example, variation in the level of Sr in hard 
parts can serve as a proxy for oscillations in temperature: 
in cooler temperatures, higher levels of Sr are substituted 
for Ca in the lattice structure of a hard part (Engstedt 
et al., 2012). Such oscillations across the optical zones of 
otoliths have been observed for the temperate species red 
emperor (Lutjanus sebae) (Seyama et al., 1991) and bluefin 
tuna (Thunnus thynnus) (Siskey et al., 2016). Therefore, 
the variations in microchemical composition of otoliths 
across optical zones can validate the seasonal frequency of 
these chemical annuli. 


Conclusions 


In addition to their use in validation of microconstituent 
composition in otoliths, one could directly age individuals 
on the basis of chemical annuli rather than optical annuli. 
Such an approach would have application for species like 
the goosefish (Lophius americanus), for which optical 
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annuli are not easily interpreted (Bank et al., 2020). For 
black sea bass, optical annuli are well aligned with chem- 
ical annuli in the 2D maps created in our study (Fig. 2), 
especially those for Mg, indicating that age can be assigned 
through direct interpretation of the 2D maps of element 
deposition in otoliths. The high-contrast visualizations of 
such analysis of otoliths from black sea bass can provide 
improved precision in age validation for optical annuli 
interpretation (Limburg et al., 2018; Heimbrand et al., 
2020; Hiissy et al., 2021b). On the other hand, for even 
moderate sample sizes (n>50), reliance on 2D maps could 
be cost-prohibitive and demand greater investment in time 
than standard aging approaches. Further, 2D map inter- 
pretations are subjective and not supported by objective 
time-series analyses, such as the periodogram approach 
applied in our study. 

Validation of direct aging through hard-part analysis 
provides the foundation for strong stock assessments. The 
results from our use of a chemical aging approach for black 
sea bass support the findings of recent marginal incre- 
ment validation studies and optical annulus interpreta- 
tions. Still, the variable interpretations of the first optical 
annulus in otoliths merit additional examination given 
the 1) varying results across the candidate elements Mg, 
Ba, P, and Mn and 2) past uncertainty in assignment of 
the first annulus. Chemical validation of optical annulus 
formation in otoliths of black sea bass further strength- 
ens the confidence in accurate age data for age-structured 
black sea bass stock assessments. 


Resumen 


En este estudio, examinamos si los microconstituyentes de 
los otolitos se depositan estacionalmente de forma similar 
a la formacion del anillo 6ptico y, por tanto, pueden uti- 
lizarse para validar las interpretaciones de la edad. En las 
especies templadas, los cambios estacionales de tempera- 
tura impulsan la formacién de anillos 6pticos, y supone- 
mos que de la misma manera causan oscilaciones en la 
deposicién de microconstituyentes. Utilizando andAlisis 
de espectrometria de masas con plasma de acoplamiento 
inductivo por ablacién laser, comprobamos la periodicidad 
en la deposicién de bario (Ba), calcio, cobre, magnesio (Mg), 
manganeso (Mn), fosforo (P), estroncio y zinc en los otolitos 
y comparamos esa periodicidad con la periodicidad de la 
zona anular (zonas 6pticas opacas y transltcidas). Para 
esta investigacion de la periodicidad de los anillos quimi- 
cos, utilizamos serrano estriado (Centropristis striata), 
una especie para la que se ha validado la formaci6n anual 
de anillos 6pticos. Se detectaron periodicidades en los per- 
files de los elementos de otolitos de lobina negra mediante 
analisis de periodogramas de Lomb-Scargle. La formacion 
de anillos 6pticos se alineé con los perfiles de Mg, pero la 
periodicidad en la deposicién de anillos de otros elemen- 
tos, Ba, Mn y P, se observo solo después del primer anillo 
6ptico, un resultado que es indicativo de los cambios onto- 
genéticos en el habitat (de aguas costeras a aguas cerca- 
nas a la plataforma) que se sabe ocurren durante la etapa 


de vida juvenil de esta especie. Ciertos elementos, como 
el Mg, identificados mediante este andlisis de otolitos y 
el andlisis del periodograma podrian aplicarse a especies 
para las que no existe un procedimiento de validacién de 
la edad. 
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Abstract—The distribution and abun- 
dance of demersal young of the year 
(YOY) of rockfish species (Sebastes spp.) 
were investigated over a 17-year period 
along the central coast of California at 
depths of 19-326 m. More than 87,000 
YOY rockfish were counted during sur- 
veys of 1292 transects on the seafloor. 
The 3 most abundant species were the 
pygmy (S. wilsoni), halfbanded (S. semi- 
cinctus), and shortbelly (S. jordani) 
rockfish. Young of the year of many 
rockfishes were found at greater depths 
than have been previously reported. The 
highest densities of all species combined 
occurred in shallow and intermediate 
depths (<100 m). Most YOY rockfish 
occupied rock or mixed substrata, except 
the greenstriped (S. elongatus) and 
stripetail (S. saxicola) rockfish, which 
were associated primarily with soft sub- 
strata. Statistically significant (P<0.01) 
hot spots of YOY rockfish densities were 
located throughout the region, but only 
Soquel Canyon, a known refuge for large 
adult rockfish (>50 cm in total length), 
was a significant (P<0.01) cold spot. Bio- 
logical processes (e.g., ontogenetic move- 
ment and predation) and environmental 
characteristics (e.g., seafloor habitat and 
oceanographic conditions) are potential 
influences on the distribution of YOY 
of rockfishes. The results of our study 
provide insights into the distributions of 
YOY of rockfishes in deep, complex sea- 
floor habitats and can help improve esti- 
mation of year-class strength through 
targeted surveys of preferred habitats. 
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Over 70 rockfish species (Sebastes spp.) 
occur in the northeast Pacific Ocean 
from Alaska to Baja California, Mexico, 
and many of them are important to 
commercial and recreational fisheries. 
Rockfish species have long lifespans 
(from decades to more than 100 years 
for certain species), are slow growing, 
and mature late (Love et al., 2002), 
making them slow to recover from high 
levels of fishing mortality. In studies 
of the distribution and abundance of 
young of the year (YOY) of rockfishes, 
large interannual fluctuations in abun- 
dance have been found (Woodbury and 
Ralston, 1991; Sakuma et al., 2006; 
Ralston et al., 2013). Such variation in 
abundance is a contributor to variable 
year-class strength, which is largely 
determined within the first year of life 
(Ralston and Howard, 1995; Hobson 
et al., 2001; Laidig et al., 2007). Annual 
indices of YOY abundance are used in 
stock assessments for several rockfish 
species (He et al., 2011; Dick et al., 2017; 
Field et al., 2017; He and Field, 2018). 
Young-of-the-year rockfish begin as 
pelagic larvae, transform to pelagic juve- 
niles, and settle to a variety of benthic 
habitats generally after 3-6 months, 
typically at shallower depths than 
those occupied by conspecific adults 


(Boehlert, 1977; Love et al., 2002). In 
nearshore areas (depths <30 m), species- 
specific habitat associations of settled 
YOY rockfish have been identified with 
kelp beds, estuaries, and man-made 
structures, such as piers and jetties 
(Anderson, 1983; Hallacher and Roberts, 
1985; Carr, 1991; Gallagher and Hep- 
pell, 2010). In offshore areas, large num- 
bers of YOY rockfish have been found 
in high-relief rocky habitats at depths 
of 30-100 m off central California and 
Oregon (Stein et al., 1992; Yoklavich 
et al., 2002), and newly settled YOY of 
3 rockfish species were commonly found 
in a study on soft sediments at depths 
of 40-100 m in Monterey Bay, located 
south of the San Francisco Bay Area in 
California (Johnson et al., 2001). 
Relatively few surveys of settled YOY 
rockfish have been conducted in high- 
relief, rocky habitats at depths >380 m 
compared to the number of surveys 
conducted in shallow nearshore areas 
(depths <30 m). The latter are often 
conducted by scuba divers or with small 
trawl nets deployed from small boats 
or from shore, and the former require 
large, specialized survey vehicles (such 
as human-occupied and remotely oper- 
ated vehicles) and support vessels, as 
well as skilled equipment operators 
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(Yoklavich et al., 2015). The researchers that have con- 
ducted the few surveys that have been carried out in rocky 
habitats at depths >30 m have reported on YOY rockfish 
in general rather than on individual species, because of 
difficulties in identifying young fish to species. As a result, 
our understanding of the abundance and habitat use of 
YOY rockfish is lacking for those species that inhabit 
depths deeper than 30 m. 

In our study, a manned submersible was used to conduct 
visual surveys of YOY of deep-living, demersal rockfishes 
along a 140-km section of the central coast of California 
from the seaward edge of kelp beds (depths ~30 m) to 
depths >300 m. Survey data were used 1) to character- 
ize species composition, 2) to determine associations of 
YOY rockfish with seafloor substrata by species and for 
all species combined, and 3) to examine spatial patterns of 
density, as a measure of abundance, relative to depth. The 
resultant data on YOY rockfish distributions will help to 
unravel the complex life history strategies of many com- 
mercially important rockfish species and to discover vital 
rockfish nursery areas in deep waters that are not often 
surveyed. 


Material and methods 


Visual surveys of demersal fish were conducted along 
the central coast of California from Big Creek (36°00’N, 
121°35°W) to Ascension Canyon (37°02’N, 122°25’W) 
(Fig. 1) by using the 4.6-m-long Delta manned submers- 
ible (Delta Oceanographics’, Torrance, CA). A total of 11 
research expeditions were completed over the 17-year 
period from 1992 through 2009 (Table 1) during late sum- 
mer and fall (typically the time of the calmest waters 
along this coastline) at depths from 19 to 326 m (the max- 
imum depth to which the submersible can travel safely is 
365 m). All surveys were conducted during daylight hours 
from 0700 to 1700. 

The submersible had space for 1 pilot and 1 scientific 
observer and was equipped with an external color video 
camera (fabricated by Delta Oceanographics) mounted 
on the starboard side above the observer viewports. 
Before 2002, tapes in the Hi8 format were used in the 
external camera, and in 2002, the tape format changed 
to MiniDV to improve video resolution. Also, in 2002, a 
handheld MiniDV video camera (DCR-VX1000, Sony 
Group Corp., Tokyo, Japan) was placed in the lowest 
viewport inside the submersible and angled down to 
document the fish that could be seen close to the sub- 
mersible, along the inside edge of a transect. Two parallel 
lasers were mounted on either side of the external color 
camera about 20 cm apart (except during the research 
expedition in 1992, when the lasers were spaced 39.5 cm 
apart). These lasers were used to estimate fish length and 
as a tool to help in determining the width of transects. 


' Mention of trade names or commercial companies is for identi- 
fication purposes only and does not imply endorsement by the 
National Marine Fisheries Service, NOAA. 
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A conductivity, temperature, and depth instrument (SBE 
19plus SeaCAT, Sea-Bird Scientific, Bellevue, WA) was 
attached to the submersible to accurately measure depth. 
From 2007 through 2009, a ring laser gyro and Doppler 
velocity log were affixed to the submersible to measure 
the distance travelled. 

Strip-transect surveys were conducted to quantify the 
fish within a 2-m-wide swath along the seafloor and from 
the seafloor to a height 2 m above it. During transect sur- 
veys, the submersible transited approximately 1 m above 
the bottom at rates of 0.4-0.9 kt (depending on bottom 
topography). The surveys were conducted for either 10 
or 15 min (depending on the year), with 2-6 transects 
surveyed during each submersible dive. The start posi- 
tion and direction of travel on transects were directed to 
the submersible pilot by a scientist aboard the support 
vessel who tracked the submersible’s position using an 
ultra-short baseline acoustic system (ORE Trackpoint II, 
EdgeTech, West Wareham, MA) and WinFrog software 
(vers. 3.1; Fugro Pelagos Inc., San Diego, CA). During tran- 
sect surveys, the submersible followed the depth contour 
with minor course corrections (e.g., for obstacle avoidance) 
left to the submersible pilot’s discretion. 

The scientific observer (all of whom were biologists 
experienced in fish identification) in the submersible 
identified (to the lowest possible taxon), enumerated, 
and estimated the length (to the nearest 5 cm in total 
length [TL] using the lasers spaced 20 cm apart) of all 
fish observed within the 2-m-wide transects and <2 m 
above the seafloor, verbally recording this information 
on an audio channel of each video tape (for more details, 
see Yoklavich et al., 2002, 2007). Few YOY rockfish were 
observed >2 m above the seafloor, and most were within 
0.5 m of the seafloor. The 2-m width of transects was 
determined by using a combination of the distance of the 
laser spots from the submersible and measurements with 
a handheld sonar. 

Young-of-the-year rockfish were identified to species 
by their body shape, length-specific pigment patterns, 
and physical characteristics (e.g., the bocaccio, S. pau- 
cispinis, has a large mouth; Butler et al., 2012). How- 
ever, many of these patterns and characteristics were not 
easily observed from the submersible or in video footage 
because of the small size of some fish and their orienta- 
tion, and these fish remained classified as unidentified 
YOY rockfish. All rockfish that were <5 cm TL were con- 
sidered YOY, as were fish <10 cm TL that were identified 
as bocaccio, olive rockfish (S. serranoides), and shortbelly 
rockfish (S. jordani). These 3 species are relatively fast 
growing during their first year and can attain lengths 
>10 cm TL by the end of summer and fall of their first 
year (Anderson, 1983; Love et al., 2002). 

All video footage was reviewed back in the laboratory. 
Data from the observer’s verbal recordings were entered 
into a database, and fish identifications and counts were 
confirmed. Sections of seafloor observed in video footage 
for at least 5 s in duration along transects were charac- 
terized as a combination of primary and secondary sub- 
stratum types. Primary substratum types covered at 
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Map of the study area along the central coast of California showing locations where dives (black 
dots) of the Delta manned submersible were made to conduct surveys of transects during 1992-2009. 
Depth contours are shown at intervals of 100 m. 


least 50% of the seafloor, and secondary types repre- 
sented at least 20% of the remaining seafloor. We used 
these substratum types, following Greene et al. (1999): 
rock outcrop (bedrock outcrops), rock pinnacle (vertical 
bedrock pinnacles), continuous flat rock (flat-topped 
bedrock), boulder (individual rocks >25 cm in diameter), 
cobble (stones 6-25 cm in diameter), pebble (stones 
<6 cm in diameter), sand, and mud. Combinations of 


primary and secondary substrata were further classified 
as hard, mixed, or soft. Any combination of cobble, peb- 
ble, boulder, flat rock, rock outcrop, or rock pinnacle was 
considered hard, and a combination of mud and sand 
was deemed soft. For the mixed category, hard and soft 
substratum types each composed either the primary 
(>50%) or the secondary (=>20%) portion of the observed 
substratum. 
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Table 1 


Year, month or period, location, number of transects surveyed, total duration in minutes of all transect surveys, area 
surveyed, depth range surveyed, and total number of young-of-the year (YOY) rockfish (Sebastes spp.) observed for each 
of the 11 research expeditions conducted with the Delta manned submersible off the central coast of California from 1992 


through 2009. 


Period Location 


Aug Soquel Canyon 

Oct Monterey Bay 

Sep Ascension and Carmel Canyons 
Sep—Oct Big Creek 

Sep Big Creek 

Oct Soquel Canyon 

Sep Soquel Canyon 
Aug—Sep Davenport and Monterey Bay 
Sep—Nov Big Creek to Soquel Canyon 
Sep—Nov Big Creek to Soquel Canyon 
Sep-Oct Soquel Canyon 


Transect length was estimated by using 1 of 3 meth- 
ods: 1) counting intervals of paired laser spots end-to-end 
in the video footage as it advanced (in 1992-1994, 1997, 
and 1998), 2) plotting transect lines in ArcGIS (vers. 10.3; 
Esri, Redlands, CA) with edited and smoothed latitudes 
and longitudes from the ultra-short baseline acoustic 
system (in 2002-2009), or 3) using distances travelled 
measured with a ring laser gyro and Doppler velocity 
log (in 2007-2009). Transect area was estimated by mul- 
tiplying the transect length by the 2-m transect width. 
Densities of YOY rockfish on each transect were calcu- 
lated by dividing the total number of YOY rockfish by the 
total transect area. 

For analyses, data from all years were combined. We 
analyzed all YOY rockfish grouped together and YOY of 
the most common species that had overall densities 20.01 
fish/100 m? (with a minimum of 64 individuals counted 
on all transects) and that occurred on at least 15 tran- 
sects (Table 2). We excluded fish identified to species in 
the subgenus Sebastomus from analyses because mul- 
tiple species (e.g., the rosy rockfish, S. rosaceus, and the 
egreenspotted rockfish, S. chlorostictus) in this group had 
varied distributions that would confound these analyses. 
A Kruskal—Wallis (KW) test, a nonparametric method 
of one-way analysis of variance, was conducted for each 
species and for all rockfishes combined to determine if 
there was a statistically significant difference in over- 
all YOY density with depth (average depths of transects 
were grouped in 10-m depth bins). A canonical correlation 
analysis was used to examine the relationships between 
densities of the most common species and environmental 
factors (i.e., average depth of transects [grouped in 10-m 
bins] and substratum type). We examined which substrata 
were used by YOY rockfish by using a chi-square test to 
test the null hypothesis that YOY of all species combined 


Number 


172 1880 
317 3170 
391 3910 


1292 13,485 


Transects 


Area Depth 
Minutes (m7) 


No. of YOY 
range (m) observed 
42 420 
90 900 
78 1170 
64 640 25,768 
81 810 29,567 
10 “Lales 5066 
30 300 12,794 
104,630 
162,694 
176,856 
17 170 1972 
638,133 


14,142 
48,772 
49,872 


95-305 
62-296 
79-316 
19-222 
22-250 
110-292 
113-296 
66—260 
23-325 
19-326 
99-288 


were distributed proportionately to the relative availabil- 
ity of the substratum types. 

We examined the spatial distribution of YOY rockfish 
densities on transects, using ArcMap in ArcGIS. We con- 
ducted hot-spot analyses (with the Spatial Analyst exten- 
sion in ArcGIS), employing the Getis—Ord Gi* statistic 
(Getis and Ord, 1992) to identify statistically significant 
spatial clusters of transects with high (hot spots) or low 
(cold spots) densities of YOY rockfish. Hot-spot analyses 
of YOY densities were conducted for each rockfish species, 
for all rockfishes and depths combined, and for all rock- 
fishes combined within 50-m depth strata. 

To calculate the Getis—Ord Gi* statistic for each fea- 
ture (i.e., transect) relative to its neighbors, inputs that 
define the spatial relationships among features and the 
size of the neighborhood are required. For each analysis, 
we chose the input of 3000 m as a fixed distance band 
that provided a reasonable spatial scale relative to the 
seafloor topography and the input of 8 or more neighbors 
for most transects (or a minimum of 6 neighbors for some 
transects). The input of a false discovery rate correction 
was also applied to adjust critical thresholds of P-values 
from analysis with the Gi* statistic to account for multiple 
testing and spatial dependency issues (Caldas de Castro 
and Singer, 2006). As a local statistic, the Gi* tests each 
feature in a data set (i.e., allows multiple testing), and 
the identification of some statistically significant hot or 
cold spots may actually be driven by random spatial pro- 
cesses (i.e., some results may falsely reject the spatial ran- 
domness null hypothesis). Spatial dependency is caused 
by the tendency of features near each other to be simi- 
lar and can artificially inflate statistical significance. We 
considered Z-scores with P-values, adjusted with a false 
discovery rate procedure, at the 99% confidence level to 
be statistically significant; at this level, we rejected the 
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Table 2 


Total number, overall density (the total number of fish observed per total area surveyed), and minimum, maximum, and 
mean depths of observed young of the year for all rockfish species (Sebastes spp.) identified during or in video footage from 
transect surveys conducted with the Delta manned submersible off the central coast of California from 1992 through 2009. 
An asterisk (*) after a common name indicates that the species is one of the most abundant species that were used in anal- 
yses. A dagger ({) after a common name indicates that the species belongs to the subgenus Sebastomus. Some of the fish 
identified as blue rockfish (S. mystinus) may have been deacon rockfish (S. diaconus) because these species have not been 
differentiated at small sizes. Mean depths are rounded to the nearest 10 because 10-m depth bins were used in analyses. 


SD=standard deviation. N/A=not applicable because only 1 individual was observed. 


Common name 


Aurora rockfish 
Pinkrose rockfish’ 
Swordspine rockfish’ 
Yelloweye rockfish 
Chameleon rockfish 
Gopher rockfish 
Black rockfish 

Copper rockfish 
Vermilion rockfish 
Bank rockfish 

Canary rockfish 
Rosethorn rockfish’ 
Bocaccio 

Olive rockfish 
Darkblotched rockfish 
Splitnose rockfish 
Greenspotted rockfish ' 
Stripetail rockfish’ 
Yellowtail rockfish’ 
Starry rockfish ' 

Rosy rockfish ' 
Greenstriped rockfish’ 
Squarespot rockfish’ 
Widow rockfish’ 

Blue rockfish’ 
Unidentified Sebastomus 
Shortbelly rockfish’ 
Halfbanded rockfish’ 
Pygmy rockfish’ 
Unidentified YOY 
Total 


Scientific name 


Sebastes aurora 
Sebastes simulator 
Sebastes ensifer 
Sebastes ruberrimus 
Sebastes phillipsi 
Sebastes carnatus 
Sebastes melanops 
Sebastes caurinus 
Sebastes miniatus 
Sebastes rufus 
Sebastes pinniger 
Sebastes heluomaculatus 
Sebastes paucispinis 
Sebastes serranoides 
Sebastes crameri 
Sebastes diploproa 
Sebastes chlorostictus 
Sebastes saxicola 
Sebastes flavidus 
Sebastes constellatus 
Sebastes rosaceus 
Sebastes elongatus 
Sebastes hopkinsi 
Sebastes entomelas 
Sebastes mystinus/diaconus 
Sebastes spp. 
Sebastes jordani 
Sebastes semicinctus 
Sebastes wilsoni 
Sebastes spp. 


Density 
(fish/100 m7) 


Number 


<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
0.01 
0.01 
0.02 
0.02 
0.02 
0.03 
0.04 
0.06 
0.08 
0.08 
0.21 
Os 
0.65 
Lone 
10.75 
13.76 


OPRWWONN RR FFE 


Min 


241 
209 
90 
72 


36 
39 
44, 
29 


tt 


28 
30 
89 


81 
41 
28 
AA, 
33 
82 
44 
30 
24 
32 
38 
39 
43 
22 


Depth 


Max 


241 
209 
90 
72 
275 
38 
45 
51 
48 


52 


78 
83 


Mean 


240 
210 
90 
70 
260 
30 
40 
50 
40 
210 
40 
170 
50 
40 
220 
200 
120 
130 
50 
80 
60 


70 
40 
40 
90 
70 
60 
70 
70 
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null hypothesis that the density of YOY rockfish is com- 
pletely spatially random. 


Results 


Surveys of 1292 transects were conducted, resulting in 
>220 h of observations (Table 1). The depths of transects 
ranged from 19 to 326 m (from roughly the outer edge of 
the kelp bed to near the maximum depth to which the sub- 
mersible can travel safely). The total area surveyed was 
638,133 m”, and the mean transect area (with all transects 
included regardless of YOY presence) was 494 m? (standard 
deviation [SD] 151.6). The mean transect area for 10-min 


transects (sample size [n]=1203) was 480 m? (SD 137.9) and 
for 15-min transects (n=89) was 647 m? (SD 137.1). 

A total of 87,687 YOY rockfish from at least 28 different 
species were counted, with an overall density of 138.8 
fish/100 m? (Table 2). Young-of-the-year rockfish were 
observed on nearly two-thirds of all transects (822 of 1292 
transects). Approximately one-quarter (26%) of YOY were 
identified to a taxon lower than the genus Sebastes 
(n=22,967), with 94% of them identified to species. The most 
abundant identified species were as follows: the pygmy 
rockfish (S. wilsont), with 11,292 individuals observed at a 
density of 1.8 fish/100 m”; the halfbanded rockfish (S. semi- 
cinctus), with 4143 individuals counted at a density of 0.7 
fish/100 m?; and the shortbelly rockfish, with 3656 
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individuals observed at a density of 0.6 fish/100 m7”. Individ- 
uals were also identified to the subgenus Sebastomus, a 
group of 11 similar-looking rockfishes characterized by 
5 large white spots along the dorsal half of their body (Love 
et al., 2002; Butler et al., 2012). Species in this subgenus are 
difficult to identify as adults and even more so as YOY; how- 
ever, in our study, rockfish were identified to 6 species in 
this group, and 3 of these species are among the most com- 
monly identified (Table 2). 

Young-of-the-year rockfish were found at depths from 
22 to 314 m. Overall, the highest densities of all taxa 
combined within 10-m depth bins (i.e., the total number 
of fish divided by the total area sampled, for each 10-m 
depth bin) were found on transects at depths between 40 
and 90 m (Figs. 2 and 3), and the densities in this depth 
range were significantly higher than they would have been 
if the depth distribution had been uniform as expected 
(KW test: P<0.0001). Most species had depth distributions 
of YOY that were significantly different from a uniform 
depth distribution (KW test: P<0.0001), except the strip- 
etail rockfish (S. saxicola), which had a broad distribu- 
tion ranging from 41 to 257 m. For the blue (S. mystinus), 
widow (S. entomelas), and yellowtail (S. flavidus) rockfish, 
YOY were found at the relatively shallow depths of 40 
or 50 m. Young-of-the-year rosy, halfbanded, squarespot 
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(S. hopkinst), shortbelly, starry (S. constellatus), and 
pygmy rockfish were observed at mean depths ranging 
between 60 and 100 m, and greenspotted and greenstriped 
(S. elongatus) rockfish were found at mean depths >100 m 
(Table 2). The species for which YOY occurred at the great- 
est depths, such as the splitnose rockfish (S. diploproa) 
and the chameleon rockfish (S. phillipsi), were too scarce 
to include in our analyses. Sampling effort (measured as 
the area surveyed) for each 10-m depth bin varied from 
>30,000 m? at a depth of 90 m to <1000 m’? at the depths 
of 10-20 m (the submersible could not sample kelp beds 
safely; therefore, few transects within this depth stratum 
were surveyed) and depths >330 m (Fig. 2). 

Results from the canonical correlation analysis indicate 
that YOY are separated by depth in 3 groups of rockfish 
species (Fig. 4). The group with the negative second canon- 
ical coefficient includes the greenstriped rockfish, green- 
spotted rockfish, and stripetail rockfish. These species 
occurred at the deepest average depths (120-130 m) and 
had generally high densities (0.3-0.4 fish/100 m7) on soft 
substrata. The group with high second canonical and low 
first canonical coefficients includes the fishes that occu- 
pied average mean depths between 60 and 100 m (ie., 
pygmy, starry, rosy, squarespot, and halfbanded rockfish), 
and YOY of all rockfishes combined (mean depth for all 


(,wW) peAesins ealy 
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Figure 2 


Densities of young-of-the-year (YOY) rockfish (Sebastes spp.) by depth (grouped in10-m bins), for (A) all YOY rockfish 
combined, (B) shortbelly rockfish (S. jordani), (C) pygmy rockfish (S. wilson), and (D) blue rockfish (S. mystinus), from 
surveys conducted with the Delta manned submersible off central California from 1992 through 2009. The black line 
in each graph indicates the total area surveyed for each depth. 
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Big Creek Pt.Lobos Port.Ledge Pt. Sur Soquel Davenport 
Location 


Figure 3 


Box plots of the mean depth at which all young-of-the-year (YOY) rockfish 
(Sebastes spp.) were observed at the 6 locations along the coast of central Cal- 
ifornia that were surveyed the most during 1992—2009. Each box shows the 
mean depth (diamond), the distribution between the first and third quartile 
(top and bottom lines), and the median depth (horizontal line). The whiskers 
indicate the maximum and minimum depths at which YOY rockfish were 
observed. Port. Ledge=Portuguese Ledge, a marine protected area located at 
the southern end of Monterey Bay. 
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Figure 4 


Comparison of first and second canonical correlation coefficients for den- 
sities of young of the year (YOY) of the 12 most abundant rockfish spe- 
cies (Sebastes spp.) and of all rockfishes combined from surveys conducted 
with the Delta manned submersible along the central coast of California 
during 1992-2009. The species with abbreviated common names are the 
greenstriped rockfish (Sebastes elongatus) (gstrip), greenspotted rockfish 
(S. chlorostictus) (grspot), stripetail rockfish (S. saxicola) (stripe), short- 
belly rockfish (S. jordani) (shortb), yellowtail rockfish (S. flavidus) (yellow), 
squarespot (S. hopkinsi) (square), and halfbanded rockfish (S. semicinctus) 
(halfba). 
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YOY was 67 m) had similar results in 
this analysis. The remaining species 
grouped in the middle of Figure 4 (i.e., 
blue, widow, yellowtail, and shortbelly 
rockfish) occurred at mostly shallow 
depths in high densities on hard sub- 
strata. Although shortbelly rockfish were 
found at a mean depth equal to that of 
the mid-depth group, over 2000 individu- 
als of this species were found at shallow 
depths of 50 m or less (Fig. 2B). 

The breakdown of seafloor substrata 
types for all transects combined (over 
the entire survey area of 638,133 m’), 
is as follows: 48% hard, 27% mixed, 
and 30% soft (Fig. 5). Most YOY rock- 
fish occupied hard and mixed substrata 
(Fig. 5). The frequency of occurrence of 
YOY for each of 7 species was dispro- 
portionally greater on hard substrata 
relative to its availability (chi-square 
test: P<0.0001), with frequencies vary- 
ing from 58% (greenspotted rockfish) 
to 96% (squarespot rockfish). For 4 of 
these species, the squarespot, widow, 
yellowtail, and blue rockfish, YOY were 
never observed on soft substrata. Each 
of 3 species, the halfbanded, pygmy, and 
rosy rockfish, occurred disproportionally 
more on mixed substrata relative to its 
availability (chi-square test: P<0.0001), 
with 75% of halfbanded rockfish occur- 
ring on this substrata. Each of 2 species, 
the greenstriped and stripetail rockfish, 
occurred disproportionally more on soft 
substrata relative to its availability (chi- 
square test: P<0.0001). 

Although YOY of rockfish species 
were present throughout the study area 
at densities ranging from 0.1/100 m’ to 
3236/100 m7, density hot spots for all 
but greenstriped rockfish were located 
within the southern region from Point 
Joe to Big Creek. For all taxa combined, 
there was a statistically significant 
(P<0.01) hot spot at Big Creek (Table 38, 
Fig. 6), where hot spots of YOY densities 
were found for individual species, 
including the stripetail, blue, shortbelly, 
and halfbanded rockfish, and for uniden- 
tified rockfish. Point Sur was a statisti- 
cally significant hot spot for 7 taxa: 
greenspotted, yellowtail, starry, rosy, 
widow, and pygmy rockfish and uniden- 
tified rockfish (hot spots for pygmy rock- 
fish are shown in Figure 6, as an 
example). Three species, the starry, 
squarespot, and pygmy rockfish, had a 
hot spot at Point Lobos. Point Joe, 
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Davenport, and Monterey Canyon were 
m Hard CMixed eSoft statistically significant hot spots for 
YOY pygmy, greenstriped, and stripetail 
rockfish, respectively. The greenstriped 
rockfish was the only species that had a 
cold spot, which was a statistically sig- 
nificant (P<0.01) cluster of transects 
with low densities located in Soquel 
Canyon inside Monterey Bay. 

In the spatial examination of the den- 
sity at depth for YOY of all rockfishes 
combined, no statistically significant hot 
spots were found at depths <50 m, depths 
of 150-200 m, or depths of 250-300 m. 
Statistically significant (P<0.01; Table 4) 
hot spots were located at Big Creek at 
depths of 50-100 m, Point Sur at depths 
of 100-150 m, and Point Lobos at depths 
of 200-250 m. A statistically significant 
Figure 5 cold spot was located at Soquel Canyon 


87,687 402 503 150 3656 147 541 86 11,292 174 4143 236 112 


Percentage of YOY on each substraum type 


Frequency of occurrence on substratum types occupied for young of the year at depths of 200-250 m. 
(YOY) of each rockfish species (Sebastes spp.) and for YOY of all rockfishes 

combined during surveys conducted with the Delta manned submersible off 

the central coast of California from 1992 through 2009. Percentages of sur- Discussion 

veyed area classified as each substratum type, combined over the entire 
study period, are given in the first bar. The number at the top of each bar is 
the number of individuals observed for that species. The total area surveyed : , 
during this study was 638,133 m”. An asterisk (*) indicates that YOY were rockfish BPECies have been categorized 
observed significantly more on this substratum type relative to what was into a single group. In this study, we 
expected. were able to identify individuals of 


28 species, because of recent progress in - 


In most previous studies, YOY of all 


Table 3 


Results from hot-spot analysis with the Gi* statistic for young of the year 
(YOY) of rockfish species (Sebastes spp.) along the central coast of California 
from 1992 through 2009. This analysis was used to identify locations of statis- 
tically significant (P<0.01) clusters of survey transects with high (hot spots) 
and low (cold spot) densities of YOY rockfish, by species. 


Location 


Species Cold spot Hot spots 


All rockfishes None Big Creek 

Blue rockfish None Big Creek 

Greenspotted rockfish None Point Sur 

Greenstriped rockfish Soquel Canyon Davenport 

Halfbanded rockfish None Big Creek 

Pygmy rockfish None Point Sur, Point Lobos, Point Joe 
Rosy rockfish None Point Sur 

Shortbelly rockfish None Big Creek 

Squarespot rockfish None Point Lobos 

Starry rockfish None Point Sur, Point Lobos 
Stripetail rockfish None Big Creek, Monterey Canyon 
Widow rockfish None Point Sur 

Yellowtail rockfish None Point Sur 

Unidentified rockfish None Point Sur, Big Creek 
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Figure 6 


Map showing the results from hot-spot analysis with the Gi* statistic for young-of-the-year (YOY) 
pygmy rockfish (Sebastes wilsoni) and for YOY of all rockfish species combined from surveys con- 
ducted with the Delta manned submersible off the central coast of California from 1992 through 
2009. This analysis was used to identify hot spots, or statistically significant (P<0.01) clusters of 


survey transects, with high YOY densities. 


the identification of rockfishes for individuals at the YOY 
stage along the West Coast of the United States (Butler 
et al., 2012). In 2002, YOY were identifiable for less than 
30 species (Love et al., 2002); currently, about 40 of the 
over 70 rockfish species have detailed, descriptive develop- 
mental series of the larval and juvenile stages, allowing 


for more accurate identifications (Anderson, 1983; Love 
et al., 2002; Butler et al., 2012). 

Most of the rockfishes that occur along the central 
coast of California have had their juvenile pigment pat- 
terns identified and documented, and many of these spe- 
cies (i.e., pygmy, halfbanded, shortbelly, and squarespot 
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Table 4 


Results from hot-spot analysis with the Gi* statistic for 
young of the year (YOY) rockfish (Sebastes spp.) along the 
central coast of California from 1992 through 2009. This 
analysis was used to identify locations of statistically sig- 
nificant (P<0.01) clusters of survey transects with high 
(hot spots) and low (cold spot) densities of YOY rockfish, 
by 50-m depth bin. The number of transects in the >300-m 
depth bin was insufficient for analysis. 


Location 


Depth (m) 


Cold spot Hot spots 
<50 None None 

>50 and <100 None Big Creek 
>100 and <150 None Point Sur 
>150 and <200 None None 

>200 and <250 Soquel Canyon Point Lobos 
>250 and <300 None None 


rockfish) are among the most numerous that we observed. 
Many of them, such as the halfbanded rockfish with its 
large, dark, diamond-shaped pigment spot on the side of 
its body and the blue rockfish with its blue coloration, are 
quite distinctive. Young-of-the-year blue rockfish have not 
been positively differentiated from the more northerly 
distributed deacon rockfish (S. diaconus) (Butler et al., 
2012). Therefore, some of the YOY rockfish identified as 
blue rockfish in our study may be deacon rockfish. Biz- 
zarro et al. (2020) found that, as adults, blue rockfish were 
more frequently observed from Monterey Bay south with 
from two-thirds to three-quarters of the catch identified 
as blue rockfish. Further studies are needed to determine 
the identification and distribution of YOY blue and dea- 
con rockfish. It is more difficult to differentiate between 
species in a few groups, including species in the subgenus 
Sebastomus, the yellowtail and olive rockfish, and a group 
consisting of the widow, pygmy, squarespot, and shortbelly 
rockfish. The species in the latter group are all tan-colored 
and, even though their pigment patterns are distinctive 
as large individuals, they have only subtle differences at 
small sizes (e.g., shortbelly rockfish are more elongated 
than the other species in this group). 

Advances in high-definition camera systems that take 
high-resolution images allow researchers to see the defin- 
itive characters needed for identification of YOY in video 
and still images. Despite these advances, species identi- 
fication for YOY rockfish during visual surveys remains 
difficult, and we were able to identify to species only one- 
quarter of the total number of individuals observed on 
transects in our study. Most of the unidentified rockfish 
were small individuals (<5 cm TL), partially obstructed 
from view, or not observed with enough resolution in 
video images to get a positive identification. However, 
on the basis of the number of individuals of each species 
that we identified, their tan color, and lack of obvious pig- 
ment patterns, we believe that most of the unknown YOY 
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rockfish were probably pygmy, squarespot, or shortbelly 
rockfish. The number of YOY rockfish that can be iden- 
tified to species can only be expected to increase in the 
future, along with increases in camera quality and res- 
olution. With these improvements in identification and 
image quality and with the advent of emerging technolo- 
gies (i.e., environmental DNA and light-field cameras), it 
may be possible in the future to identify most YOY rock- 
fish to species during visual surveys. 

Rocky substrata have been shown to be important hab- 
itats for YOY rockfishes living in shallow water (depths 
<30 m) (Anderson, 1983; Carr, 1991; Laidig et al., 2007). 
In our study, we observed the highest densities of YOY 
rockfish on hard or mixed substrata at depths between 40 
and 90 m and the lowest densities at depths >160 m or 
over soft bottom habitat. Love et al. (2006) found the high- 
est densities of YOY bocaccio on oil platforms (hard man- 
made substrate) and on shell mounds near platform bases 
in Southern California at depths between 25 and 80 m. 
High abundances of YOY rockfish have also been seen in 
rocky substrata at depths <100 m in southeast Alaska and 
in Oregon, with low numbers of YOY rockfish observed 
in soft sediments (Carlson and Straty, 1981; Stein et al., 
1992). Although densities were not available from these 
studies, these results indicate that rocky substrata at 
these depths may be important to YOY of rockfishes along 
the West Coast. 

For many rockfish species, YOY had statistically sig- 
nificant associations with rocky areas, but soft sediments 
served as a habitat for YOY greenstriped and stripetail 
rockfish. Not surprisingly, the highest densities of these 
2 species were located in areas with high amounts of 
muddy seafloor (at Monterey Canyon and Davenport), at 
the deep average depths of 120 m for greenstriped rock- 
fish and 130 m for stripetail rockfish. Soft sediments 
are also the common habitats of conspecific adults, but 
adults have been found at deeper depths than YOY 
(Love et al., 2002). Young-of-the-year stripetail rockfish 
have been observed in large numbers on soft sediments 
off Big Creek in central California at depths of about 
70 m and in Southern California at depths of 18-150 m 
(Mearns et al., 1980; Yoklavich et al., 2002), and, for 
YOY, both the greenstriped rockfish and the stripetail 
rockfish have been found to be the most abundant rock- 
fishes surveyed on the soft sediments of Monterey Bay 
from the depth of 60 m to the 100-m limit of their sam- 
pling (Johnson et al., 2001). 

Most settled YOY of rockfishes live at shallower depths 
than their adult counterparts. Ontogenetic shifts in depth 
are most obvious in the deepwater species that recruit ini- 
tially to kelp beds (e.g., the shortbelly and yellowtail rock- 
fish; Anderson, 1983; Love et al., 1991, 2002, 2009; Butler 
et al., 2012). The YOY ofa handful of kelp-bed-settling spe- 
cies remain in shallow depths (e.g., kelp rockfish, S. atro- 
virens, and black-and-yellow rockfish, S. chrysomelas), and 
YOY of most species either remain for only a short time 
(days to a month; e.g., halfoanded and shortbelly rockfish) 
or for many months until the first big storms of winter hit 
in October or November (e.g., copper rockfish, S. caurinus; 
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stripetail rockfish; and yellowtail rockfish; Anderson, 
1983). Most YOY blue and olive rockfish move to deeper 
water in fall, but some individuals remain in kelp beds 
throughout the winter. 

Although these ontogenetic habitat movements have 
been recorded in nearshore areas, the movements of only 
a few species that recruit to deeper waters have been 
observed. Johnson et al. (2001) used a size-frequency analy- 
sis to document the ontogenetic movement to deeper waters 
for cowcod (S. levis), greenstriped rockfish, and stripetail 
rockfish in Monterey Bay. Love et al. (1991) found that 
bocaccio move successively deeper as they mature, with a 
shallow-water-dwelling juvenile taking around 5 years to 
move to deeper adult habitats. In our study, YOY rockfish 
were observed living in habitats similar to, but shallower 
than, habitats utilized by conspecific adults, according to 
data on depth distributions from Love et al. (2002). 

The average depth at which YOY of most rockfishes were 
found was over 50 m shallower than the average depth at 
which adults have been found. For example, YOY widow 
and yellowtail rockfish in our study were found at aver- 
age depths of 40 and 50 m, respectively, and in contrast 
the average depth of adults has been reported as well over 
100 m (Love et al., 2002). However, some species make an 
ontogenetic depth change of well over 100 m; for example, 
stripetail rockfish can settle in depths as shallow as 15 m 
but as adults live at 200 m (Love et al., 2002). 

Some rockfish species recruit to nearshore kelp beds, but 
many of these same species also recruit to deeper depths. 
Young-of-the-year blue, widow, and yellowtail rockfish 
have been recorded in large numbers in nearshore kelp 
beds (Anderson, 1983; Carr, 1991; Laidig et al., 2007). 
These beds reach maximum depths of 20-30 m depend- 
ing on water clarity (Schiel and Foster, 2015). The highest 
densities of blue rockfish in our study occurred at depths 
around 30—40 m, and widow and yellowtail rockfish were 
most abundant at depths of 40-50 m. These species utilize 
both shallow nearshore areas (depths <30 m) and deeper 
rocky habitats. Carr (1991) counted YOY blue rockfish in 
kelp beds and estimated densities of 180-400 fish/100 m? 
over a 2-year period. Laidig et al. (2007) surveyed YOY 
of rockfishes in the kelp beds off Northern California and 
found large interannual fluctuations in abundance for 
YOY blue and yellowtail rockfish with densities varying 
from 0.07 to 655 fish/100 m7”. In this study, the density of 
blue rockfish at the depth of 30 m was 2 fish/100m?, and 
the density of yellowtail rockfish at the same depth was 
0.1 fish/100 m?. The submersible could not survey within 
the kelp bed; therefore, a substantial proportion of each of 
the populations of blue and yellowtail rockfish may have 
been missed. 

The differences in densities may also be explained 
through the interannual variation in recruitment strength 
realized in all rockfishes (Carr, 1991; Ralston et al., 2013). 
Some species that recruit to the kelp bed, including the 
bocaccio and the black (S. melanops), gopher (S. carnatus), 
olive, and vermilion (S. miniatus) rockfish, were rarely 
seen at deeper depths in this study. The YOY of these spe- 
cies may be more restricted to shallow depths than YOY of 


other rockfish species. Other species have more protracted 
distributions that stretch from the kelp bed to depths 
>100 m. Stripetail, halfbanded, and shortbelly rockfish 
all occasionally recruit to kelp beds (Anderson, 1983), but 
distributions of their YOY occur at much deeper depths. 
Love et al. (2009) surveyed rockfishes in Southern Califor- 
nia and found YOY halfbanded rockfish in areas ranging 
from the kelp bed (<30 m) to depths over 200 m. Therefore, 
recruitment numbers for many species in kelp bed envi- 
ronments provide only a relative measure of their abun- 
dance, and the numbers of these species living in deeper 
waters should be considered. 

The abundance of YOY pygmy, halfbanded, and short- 
belly rockfish in our study reflects the dominance of adults 
of these dwarf species in the fish assemblages of deepwa- 
ter rocky habitats off California (Yoklavich et al., 2002; 
Baskett et al., 2006). Love et al. (2009) found that the 4 
most abundant species in Southern California at depths 
between 19 and 365 m were the squarespot, halfbanded, 
shortbelly, and pygmy rockfish. In a study of marine pro- 
tected areas (MPAs) along the central coast of California, 
Starr and Yoklavich (2008) observed that at most depths 
between 100 and 200 m, the dominant species group was 
dwarf rockfishes. In another study at the Big Creek EKco- 
logical Reserve also along the central coast of California, 
Yoklavich et al. (2002) observed that the halfbanded and 
shortbelly rockfish were the 2 dominant fish species as 
adults within the reserve. 

Although dwarf species dominate many reefs, some pro- 
tected areas harbor many old fish of larger species. An 
example of a relatively pristine refuge is Soquel Canyon. 
There, considerable numbers of large (>50 cm TL), old 
individuals of piscivorous rockfishes have been reported, 
but, not surprisingly, few adults of dwarf rockfishes have 
been observed there (Yoklavich et al., 2000). As expected, 
with so many large predators there, few YOY rockfish were 
observed in Soquel Canyon during this study. Much of the 
deep rocky habitat off the U.S. West Coast has been heavily 
fished, and fishing has disproportionately removed rock- 
fish of large, slow-growing, piscivorous species that prey 
on small, fast-growing dwarf rockfish species (O’Farrell 
et al., 2009). With the decline in abundance of larger 
species, populations of dwarf species have increased in 
abundance. The recovery of large rockfish species may be 
slowed by competition with these highly abundant dwarf 
rockfish species (an outcome referred to as the cultivation 
effect; Walters and Kitchell, 2001). 

Initial settlement of pelagic juvenile rockfish to benthic 
areas is dependent on many factors, including upwell- 
ing, fronts, depth, available habitat, and temperature 
(Anderson, 1983; Love et al., 2002; Laidig, 2010). Upwell- 
ing has been suggested as a settlement mechanism for 
winter-spawned rockfish off California, with pelagic juve- 
niles actively descending through the water column into 
the deeper shoreward currents that transport them to the 
appropriate shallow habitats (Lenarz et al. 1991; Larson 
et al., 1994; Johnson et al., 2001). Markel and Shurin 
(2020) observed high recruitment of black rockfish off Can- 
ada during a year with strong upwelling. Conversely, they 
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ascertained that periods of downwelling were correlated 
to higher settlement of other species: the copper, quillback 
(S. maliger), and brown (S. auriculatus) rockfish. Woodson 
et al. (2012) found that ocean fronts in the California 
Current Large Marine Ecosystem were positively cor- 
related with recent settlement of winter-spawned fish of 
Sebastes species. Larson et al. (1994) found an unusually 
shallow distribution of late-stage pelagic juvenile rock- 
fish above the thermocline in Monterey Bay and sug- 
gested that this behavior may contribute to shoreward 
transport of these large juveniles. They also described a 
multistage settlement process, with moving to shore as 
the first stage and finding suitable settlement habitat as 
the second stage. 

The distribution of settled YOY rockfish is determined 
by the availability of food and shelter, mortality, the 
number of rockfish that initially settled there, or the 
redistribution of YOY rockfish on different substrata 
(Love et al., 1991). Once settled, predator avoidance 
(West et al., 1994; Hobson et al., 2001; Love et al., 2002) 
and food availability (Hallacher and Roberts, 1985; Love 
et al., 1991) are likely factors affecting distribution. 
Stein et al. (1992) observed clouds of YOY rockfish near 
the top of a rocky bank (at depths of 60—80 m) and larger 
adults in deeper areas (at depths of 100-150 m), indi- 
cating that predation is a potential cause for this dis- 
tribution. The tops of rocky banks have more ambient 
light than deeper areas, aiding in predator detection and 
avoidance. In a laboratory study, predator avoidance by 
YOY black rockfish was higher when the predator was 
visible to the YOY rockfish than when they had only 
chemical cues (Giger, 2015). 

In a study of artificial reefs, West et al. (1994) determined 
that abundance of demersal YOY rockfish was greatest on 
the most complex reefs, with aggregations near crevices. 
Small rockfish near the top of reefs may also be feeding 
on planktonic organisms in open waters above the reef 
more than from the reef itself (Carlson and Haight, 1976; 
Singer, 1985). Love et al. (1991) found evidence that food 
resources were less abundant in deeper, offshore areas 
than in shallow reefs. Besides food resources, temperature 
may be a factor in determining YOY rockfish distributions. 
Boehlert and Yoklavich (1983) studied YOY black rock- 
fish in the laboratory and determined that growth rates 
increased at warmer temperatures. The YOY rockfish that 
settle in shallow zones may have an advantage by growing 
faster than deeper-living conspecifics. 

The spatial distribution of settled YOY rockfish varied 
by species within our study area. Big Creek was the hot 
spot for YOY of all rockfishes combined, with the identi- 
fication of this hot spot driven by high densities of half- 
banded and shortbelly rockfish and of unidentified YOY 
rockfish there during 1997 and 1998. These taxa along 
with the blue and stripetail rockfish also each had sta- 
tistically significant hot spots at this location. For several 
species (i.e., the greenspotted, pygmy, rosy, starry, widow, 
yellowtail rockfish) and unidentified rockfish, hot spots of 
high densities were found at Point Sur, a major upwelling 
center (Traganza et al., 1981). 
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Four other locations were also hot spots for 1 or 2 spe- 
cies: Davenport, Point Joe, Point Lobos, and Monterey 
Canyon. None of these locations was a hot spot for all rock- 
fish species combined (Table 3), and the oceanographic 
currents differ between locations. Therefore, it is unlikely 
that ocean current transport alone is the basis for these 
distributions because many of these species aggregate 
together as pelagic juveniles before settling to the seafloor 
(Lenarz et al., 1991; Larson et al., 1994; Field et al., 2021). 
If they had all settled together and not migrated after set- 
tlement, we would expect hot spots to be populated by the 
species that aggregated together in the planktonic envi- 
ronment, and that is not the case. In a study of pelagic 
juveniles of rockfishes off central California during 2001-— 
2019, only the shortbelly rockfish was in the top ten of 
the most abundant rockfish species over the 19-year study 
(Field et al., 2021). If ocean currents brought fish to these 
locations, we would expect the most abundant species of 
settled YOY rockfish to be proportional to the abundances 
of pelagic juvenile rockfish in the same areas, and they 
are not. It is more likely that the fish migrated after set- 
tlement to specific regions of high food availability or of 
low predation at appropriate species-specific depths and 
seafloor habitats (Love et al., 1991; Hobson et al., 2001; 
Auster et al., 2003). Those individual YOY rockfish that 
remain in the less desirable areas may perish or experi- 
ence reduced growth rates. 

Some of the areas of highest abundance of YOY rock- 
fish were located within or near an MPA. The fish living 
in these areas receive increased protection from fishing 
and other human activities, leading to increased num- 
bers of rockfish that can translate into greater larval 
production and greater numbers of YOY rockfish settling 
there. Most of the MPAs in our study area were estab- 
lished in 2007 Gnformation on the MPAs in California 
are available on the website of the California Depart- 
ment of Fish and Wildlife); therefore, they were too new 
to have a significant effect on larval production. However, 
2 MPAs were established prior to 2007, the Point Lobos 
and Big Creek Ecological Reserves in 1973 and 1994, 
respectively (McArdle, 1997). There may have been some 
additional contribution of larval rockfish from within 
these 2 reserves, but these reserves are small and their 
contribution may have been minimal. Also, high numbers 
of YOY were observed at Big Creek in 1997, only 3 years 
after the reserve was established, likely not enough time 
for rockfish abundance to increase enough to result in a 
substantial increase in the number of larvae produced. In 
the future, the added larval production from these MPAs 
should increase the numbers of YOY rockfish recruiting 
along the coast of central California. 

A better understanding of the habitats and areas of 
importance to YOY rockfish can improve estimates of year- 
class strength and management measures to protect them. 
Marine protected areas, gear restrictions, and other pro- 
tective measures can be employed to reduce human effects 
in these areas, which may prove to be nursery areas where 
juveniles may be able to thrive and survive to maturity, 
increasing their contribution to the population. Year-class 
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strength is a valuable component of assessing the status of 
fish populations and forecasting future growth. Indices of 
year-class strength for YOY rockfishes have traditionally 
been estimated by using catch-per-unit-of-effort data from 
midwater trawl and scuba surveys (Laidig et al., 2007; 
Ralston et al., 2013). Visual surveys conducted routinely 
in deep, complex rock habitats can enhance these tradi- 
tional indices and provide indices of year-class strength 
for YOY rockfish species not represented in traditional 
surveys. Habitat-specific estimates of YOY abundance 
improve our ability to manage rockfish populations on an 
ecosystem basis. 


Conclusions 


Young-of-the-year rockfish were distributed throughout 
the central coast of California at depths of 22-314 m. In 
the visual surveys of YOY in this study, 28 species of rock- 
fish were identified, and for many this study was the first 
time they were identified in situ, leading to new depth dis- 
tributions and habitat associations for them. The highest 
YOY densities for rockfishes were found in hard or mixed 
habitats at depths of 40-90 m. For YOY, rockfishes fell into 
3 groups, the group found on soft substrata in deep water, 
the mid-depth group, and the shallow-water group. The 
shallow-water group consisted of many species that also 
recruit to kelp beds. For 2 of the species that recruit to 
shallow kelp beds, the blue and yellowtail rockfish, den- 
sities were generally higher in these beds than at deeper 
depths. However, the maximum depths at which these 
2 species were observed in this study were 81 and 50 m, 
respectively, depths that are much deeper than the edge of 
the kelp beds (depths of ~30 m). Therefore, estimating spe- 
cies densities from surveys of kelp beds alone may result 
in underestimation of population density for YOY of spe- 
cies that recruit to kelp beds. 

The hot spots of YOY rockfish densities coincided with 
3 MPAs at Big Creek, Point Sur, and Point Lobos. These 
MPAs offer protection for these young fish, potentially 
allowing more individuals to reach maturity. Understand- 
ing of the species-specific distributions and habitat pref- 
erences of rockfish in these early life stages will improve 
fisheries management by enhancing 1) the reliability of 
density estimates with sampling of correct depths and 
habitats, 2) the protection of potential nursery areas and 
placement of future MPAs, and 3) estimates of year-class 
strength for stock assessment models. 


Resumen 


Se investigé la distribucién y abundancia de juveniles 
demersales del afio de especies de peces de roca (Sebastes 
spp.) durante un periodo de 17 anos a lo largo de la costa 
central de California a profundidades de 19-326 m. Se 
contaron mas de 87,000 juveniles de peces de roca durante 
los estudios de 1292 transectos en el fondo marino. Las 
3 especies mas abundantes fueron el rocote pigmeo 


(S. wilsoni), el rocote inspector (S. semicinctus) y rocote 
pancita (S. jordani). Los juveniles del afio de muchas 
especies de peces de roca se hallaron a profundidades 
mayores que las registradas anteriormente. Las mayores 
densidades de todas las especies combinadas se dieron en 
profundidades someras e intermedias. La mayoria de los 
juveniles ocuparon sustratos rocosos 0 mixtos, excepto el 
rocote reina (S. elongatus) y el rocote cola listada (S. saxi- 
cola), que se asociaron principalmente con sustratos blan- 
dos. En toda la region se localizaron puntos calientes de 
densidades de juveniles de peces de roca estadisticamente 
significativos (P<0.01), pero sdlo el canon Soquel, un con- 
ocido refugio de grandes peces de roca adultos (>50 cm de 
longitud total), fue un punto frio significativo (P<0.01). 
Los procesos biolégicos (ej., el movimiento ontogénico y 
la depredacion) y las caracteristicas ambientales (ej., el 
habitat del fondo marino y las condiciones oceanograficas) 
pueden influir en la distribucién de juveniles del afio de 
peces de roca. Los resultados de nuestro estudio propor- 
cionan informacioén sobre la distribucién de juveniles del 
ano de peces de roca en habitats profundos y complejos del 
fondo marino, y pueden ayudar a mejorar la estimacién de 
la fuerza de la cohorte anual mediante estudios especificos 
de los habitats preferidos. 
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Guidelines for authors 


Contributions published in Fishery Bulletin describe 
original research in marine fishery science, fishery engi- 
neering and economics, and marine environmental and 
ecological sciences (including modeling). Preference will 
be given to manuscripts that examine processes and 
underlying patterns. Descriptive reports, surveys, and 
observational papers may occasionally be published but 
should appeal to an audience outside the locale in which 
the study was conducted. 

Although all contributions are subject to peer review, 
responsibility for the contents of papers rests upon the 
authors and not on the editor or publisher. Submission 
of an article implies that the article is original and is not 
being considered for publication elsewhere. 

All submissions are subject to a double-blind review 
to remove the identity of author and reviewer during the 
review process. 

Plagiarism and double publication are considered 
serious breaches of publication ethics. To verify the orig- 
inality of the research in papers and to identify possible 
previous publication, manuscripts may be screened with 
plagiarism-detection software. 

Manuscripts must be written in English; authors 
whose native language is not English are strongly advised 
to have their manuscripts checked by English-speaking 
colleagues before submission. 

Artificial intelligence (AI), machine learning, or similar 
algorithmic tools should not be used to write or generate 
scientific content, including text, tables, and figures, for 
manuscripts submitted to Fishery Bulletin. Although some 
authors may use such tools to improve the readability or 
syntax of their manuscript, it must be done with extreme 
caution, and authors must review and edit the result. If 
such a tool is used, authors must declare the name of the 
tool, where it was used, and the purpose for which it was 
used in an “Artificial intelligence declaration” section after 
the “Acknowledgments” section. Authors remain account- 
able for all elements of their manuscript and should state 
that they take full responsibility for the content of their 
manuscript in the AI declaration. 

Once a paper has been accepted for publication, online 
publication takes approximately 2 weeks. 

There is no cost for publication in Fishery Bulletin. 


Types of manuscripts accepted by the journal 


Articles generally range from 20 to 30 double-spaced 
typed pages (12-point font) and describe an original 
contribution to fisheries science, engineering, or econom- 
ics. Tables and figures are not included in this page count, 
but the number of figures should not exceed 1 figure for 


every 4 pages of text. Articles contain the following divi- 
sions: abstract, introduction, methods, results, discussion, 
and conclusions. 


Notes are generally less than 10 double-spaced typed 
pages (12-point font), including the Literature cited sec- 
tion. Like articles, notes describe an original contribu- 
tion to fisheries science. They follow a format similar to 
that for articles: abstract, introduction, methods, results, 
and discussion, but the results and discussion sections 
may be combined and a conclusions section should not 
be included. They should include no more than 2 figures 
or tables (2 of each would be too many). They are distin- 
guished from full articles in that they report a noteworthy 
new observation or discovery—such as the first report of 
a new species, a unique finding, condition, or event that 
expands our knowledge of fisheries science, engineering, 
or economics—and do not require a lengthy discussion. 
Manuscripts on range extensions will not be considered. 


Companion articles should be submitted together and 
are published together as a scientific contribution. Both 
articles should address a closely related topic and may be 
articles that result from a workshop or conference. 


Review articles of exceptional quality may be considered 
occasionally for publication on a case-by-case basis. They 
address a timely topic that is relevant to aspects of fisheries 
science. They should include an abstract, but the format of 
the article, per se, will be up to the author. Please contact 
the scientific editor to discuss your ideas regarding a poten- 
tial review article before embarking on such a project. 


Perspective essays generally are less than 5 double- 
spaced typed pages (12-point font) and provide clarifica- 
tion of misused terms or misunderstood topics in fisheries 
science. Essays are accepted at the discretion of the 
scientific editor. Opinion pieces on policy will not be con- 
sidered. They should include an abstract and a list of the 
literature cited, but no other sections should divide the 
article. Please contact the scientific editor to discuss your 
idea before making a submission. 


Preparation of manuscript 


Title page should include authors’ full names, affiliations, 
mailing addresses, and the senior author’s email address. 


Abstract should be limited to 200 words (one-half typed 
page), state the main scope of the research, and empha- 
size the authors’ conclusions and relevant findings. Do 
not review the methods of the study or list the contents of 
the paper. Because abstracts are circulated by abstract- 
ing agencies, it is important that they represent the 
research clearly and concisely. 


General text must be typed in 12-point Times New Roman 
font throughout. A brief introduction should convey the 
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broad significance of the paper; the remainder of the paper 
should be divided into the following sections: Materials and 
methods, Results, Discussion, Conclusions, and Acknowl- 
edgments. Headings within each section must be short, 
reflect a logical sequence, and follow the rules of subdivision 
(i.e., there can be no subdivision without at least 2 subhead- 
ings). The entire text should be intelligible to interdisciplin- 
ary readers; therefore, all acronyms, abbreviations, and 
technical terms should be written out in full and defined 
the first time they are mentioned. Abbreviations should be 
used sparingly because they are not carried over to index- 
ing databases and slow readability for those readers out- 
side a discipline. They should never be used for the main 
subject (e.g., species or method) of a paper. 

For general style, follow the U.S. Government Publish- 
ing Office Style Manual (2016, available at website) and 
Scientific Style and Format: the CSE Manual for Authors, 
Editors, and Publishers (2014, 8th ed.) published by the 
Council of Science Editors. For scientific nomenclature, 
use the current edition of the American Fisheries Soci- 
ety’s (AFS) Common and Scientific Names of Fishes from 
the United States, Canada, and Mexico and its companion 
volumes (Crustaceans, Mollusks, Cnidaria and Ctenophora, 
and World Fishes Important to North Americans). For spe- 
cies not found in the previously mentioned AFS publica- 
tions and for more recent changes in nomenclature, use the 
Integrated Taxonomic Information System (ITIS, available 
at website), or, secondarily, the California Academy of Sci- 
ences Catalog of Fishes (available at website) for species 
names not included in ITIS. Common (vernacular) names of 
species should be lowercase. Citations must be given for the 
identification of specimens. For example, “Fish species were 
identified according to Collette and Klein-MacPhee (2002); 
sponges were identified according to Stone et al. (2011).” 

Dates should be written as follows: 11 November 2018. 
Measurements should be expressed in metric units, for 
example, “58 metric tons (t);” if other units of measurement 
are used, please make this fact explicit to the reader. Use 
numerals, not words, to express whole and decimal num- 
bers in the general text, tables, and figure captions (except 
at the beginning of a sentence). For example, “We consid- 
ered 3 hypotheses. We collected 7 samples in this location.” 
Use American spelling. Refrain from using the shorthand 
slash (/), an ambiguous symbol, in the general text. 

Cite all software, special equipment, and chemical solu- 
tions used in the study within parentheses in the general 
text, including the version number, company name, and the 
city and state (or nation) of the company headquarters, for 
example, “SAS, vers. 9.4 (SAS Institute Inc., Cary, NC).” 


Word usage and grammar that may be useful are the 
following: 


e Aging 
For our journal, the word aging is used to mean both 
age determination and the aging process (senescence). 
Authors should make clear which meaning is intended 
where ambiguity may arise. 


e Fish and fishes 
The plural of the word fish (a collective noun that 
implies individuals without regard to species) is fish. 
Example: The fish were collected by trawl net. 
Example: The numbers of fish collected that season 
were less than the numbers from previous years. 


The plural for fish species is fishes (a contrived plural 
used by taxonomists to mean several or more fish spe- 
cles) or one can use fish species (which is preferred in 
this journal for clarity across disciplines). 
Example: The fishes of Puget Sound [biodiversity is 
implied] or 
Example: The fish species of Puget Sound [preferred 
plural for clarity across disciplines]. 


e Crab and crabs, squid and squids, ete. 
The plural of the word crab (i.e., many individuals 
without regard to species) is crab. 
Example: The crab were sorted by weight. 
Example: Many red king crab were dying [Many 
individuals of one species of crab.] 


The plural of crab species is crabs (a word used by tax- 
onomists) or crab species (the latter is preferred in this 
journal for clarity). 
Example: These crabs were selected for treatment. 
[Different crab species are implied.] 
Example: These crab species were selected for 
treatment. [Preferred word choice for clarity.] 
Example: Snow crabs are found throughout the 
North Pacific Ocean and Bering Sea. [There are 
2 species of snow crab; therefore the word crabs can 
be used here. | 
Example: Two species of snow crab are found through- 
out the North Pacific Ocean and Bering Sea. [Pre- 
ferred usage for clarity.] 
Example: Three crabs were selected for treatment. 
[3 species of crab are implied.] 
Example: ‘Three crab species were selected for treat- 
ment. [Preferred word choice for clarity.] 


e We use fisherman and fishermen, not fisher and fishers, 
in this journal. One can always use crew member, vessel 
operator, and angler (the latter for recreational fishing). 


e The definite article with common names of species 
When the singular common name of a species rep- 
resents the entire class or group to which it belongs, 
use the definite article. 

Example: Only one species of the genus Salmo is 
found in the Atlantic Ocean—the Atlantic salmon 
(Salmo salar). 

Example: The sonic emissions of the bottlenose dol- 
phin are complex. 


For plural common names, this rule does not apply. 
Example: Chinook salmon are found throughout the 
Pacific Ocean. 
Example: Bottlenose dolphins are found in temper- 
ate and tropical waters. 


Guidelines for authors 


e Sex 
For the meaning of male and female, use the word sex, 
not gender. Do not write, “fish were sexed.” Write, “sex 
was determined.” 


e Participles 
As adjectives, participles must modify a specific noun 
or pronoun. 

Example: Using mark-recapture methods, these sci- 
entists determined the size of the population. [Correct. 
The participle using modifies the word scientists.| 

Example: These scientists, based on the collected 
data, concluded that the mortality rate of these fish 
had increased. [Incorrect. The scientists were not 
based on the collected data.] 

Example: ‘These scientists concluded, on the basis of 
collected data, that the mortality rate of these fish 
had increased. [Correct. The offending participle has 
been eliminated and an adverbial phrase modifies 
the verb concluded.| 


Equations and mathematical symbols should be set froma 
standard mathematical program (MathType or Equation 
Editor). Equations formatted in LaTex are not accept- 
able. For mathematical symbols in the general text (a, 7, 
T, +, etc.), use the symbols provided by the MS Word pro- 
gram and italicize all variables, except those variables 
represented by Greek letters and the superscript and 
subscript parts of variables and expressions. Do not use 
photo mode when creating these symbols in the general 
text, and do not cut and paste equations, letters, or sym- 
bols from a different software program. 

Number equations (if there are more than one) for 
future reference by scientists; place the number within 
parentheses at the end of the first line of the equation. 


Literature cited section comprises published works and 
those accepted for publication (in press) in peer-reviewed 
journals. Follow the name and year system for citation for- 
mat in this section (i.e., citations should be listed alpha- 
betically by the authors’ last names, and then by year if 
there is more than one citation by the same author. A list 
of abbreviations for citing journal titles can be found on 
our website. 

Authors are responsible for the accuracy and com- 
pleteness of all citations. Avoid the use of multiple 
citations when a single citation sufficiently supports a 
statement; cite the work that first reported the informa- 
tion that supports a statement, not all of the subsequent 
works. Literature citation format: Authors (last name, 
followed by initials for first name and, if given, mid- 
dle name of first author; then list names of additional 
authors with initials before last names). Year. Title of 
article. Abbreviated title of the journal in which it was 
published. Always include either the range of page num- 
bers (for a journal article) or a total number of pages (for 
a book or other type of publication). List a sequence of 
citations in the general text chronologically, for example, 
“(Smith, 1932; Green, 1998; Smith and Jones, 2015).” 
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Acknowledgments should be no more than 6 lines of 
text. Only those who have contributed in an outstanding 
way should be acknowledged by name. For recognition of 
other persons or groups, use a general term, such as crew, 
observers, or research coordinators, and do not include 
names with these terms. 


Digital object identifier (doi) code ensures that a publica- 
tion has a permanent location online. A doi link (which 
may include a doi code) should be included at the end of 
citations of published literature. Authors are responsi- 
ble for submitting accurate doi links. Faulty links will be 
deleted at the page-proof stage. 


Footnotes are used for all documents that have not been 
formally peer reviewed and for observations and personal 
communications, but these types of references should be 
cited sparingly in manuscripts submitted to the journal. 

All reference documents, administrative reports, inter- 
nal reports, progress reports, project reports, contract 
reports, personal observations, personal communications, 
unpublished data, manuscripts in review, and council meet- 
ing notes are footnoted in 10-point font and placed at the 
bottom of the page on which they are first cited. Footnote 
format is the same as that for formal literature citations. A 
link to the online source (e.g., [Available from http://www... , 
accessed July 2017.]), or the mailing address of the agency 
or department holding the document, should be provided 
so that readers may obtain a copy of the document. 


Tables are often overused in scientific papers; it is seldom 
necessary to present all the data associated with a study. 
Tables should not be excessive in size and must be cited 
in numerical order in the text. Headings should be short 
but ample enough to allow the table to be intelligible on 
its own. 

All abbreviations and unusual symbols must be 
explained in the table legend. Other incidental com- 
ments may be footnoted with numeral footnote markers. 
Use asterisks only to indicate significance in statistical 
data. Do not put a table legend on a page separate from 
the table; place the legend above the table. Do not submit 
tables in photo mode. 


e Note probability with a capital, italic P. 


e Provide a zero before all decimal points for values less 
than one (e.g., 0.07). 


e Round all values to 2 decimal points. 


e Use a comma in numbers of 5 digits or more (e.g., 
13,000 but 3000). 


Figures must be cited in numerical order in the text. 
Graphics should aid in the comprehension of the text, but 
they should be limited to presenting patterns rather than 
raw data. The number of figures should not exceed 1 figure 
for every 4 pages of text. 
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Figure legends should explain all symbols and abbrevi- 
ations seen in the figure and should be double spaced on a 
separate page at the end of the manuscript. 

Line art and halftone figures should be saved at res- 
olutions >600 dots per inch (dpi) and >300 dpi, respec- 
tively. Color is allowed in figures to show morphological 
differences among species (for species identification), to 
show stain reactions, to show gradations (such as those of 
temperature and salinity within maps), and to distinguish 
between numerous lines and symbols in graphs. Figures 
approved for color should be saved in CMYK format. 

All figures must be submitted as PDF, TIFF, or EPS 
files. 


e Capitalize only the first letter of the first word and 
proper nouns in all labels within figures. 


e Do not use overly large font sizes for labels in maps and 
for axis labels in graphs. 


e Use the same point size for all labels, except for panel 
labels (e.g., A), which should be slightly larger than 
other labels (e.g., 11 point versus 8 point). 


e Use a sans serif font for all labels. Panel labels (e.g., A), 
however, should be in Times New Roman font. 


e Do not use bold fonts or bold lines in figures. Do not 
use italic fonts. Exceptions include use of italic fonts 
for labels of bodies of water in maps and a bold font for 
panel labels (e.g., A). 


e Do not place outline rules around graphs. 


e Do not include vertical and horizontal lines in the back- 
ground of graphs. Ticks for values on the x-axis and 
y-axis will suffice. 


e Place a north arrow and label degrees latitude and lon- 
gitude (e.g., 170°E) in all maps. If scale of map requires 
more than degrees, use degrees minutes, not decimal 
degrees. 


e Place panel labels (e.g., A, B, C) within the upper-left 
area of each graph or photo in a multi-panel figure, 
from left to right, then top to bottom. If the letter is 
not visible against a dark background, put a white box 
behind it. Do not use white labels. 


e Avoid placing labels vertically or diagonally. Y-axis 
labels can be vertical. Words in horizontal labels can be 
stacked vertically to fit. 


e Use symbols, shading, or patterns (not clip art) in maps 
and graphs. 


e For scale bars in maps, use kilometers. Use the label 
km or kilometers (lowercased). 


Supplementary materials that are considered essential, 
but are too large or impractical for inclusion in a paper 
(e.g., metadata, figures, tables, videos, and websites), may 


be provided at the end of an article. These materials are 
subject to the editorial standards of the journal. A URL to 
the supplementary material and a brief explanation for 
including such material should be sent at the time of ini- 


_ tial submission of the paper to the journal. 


¢ Metadata, figures, and tables should be submitted in 
standard digital format (MS Word or PDF file) and 
should be cited in the general text, for example, as 
“,.. was determined (Suppl. Table 3, Suppl. Fig. 1).” 


e Websites should be cited with a URL in the general 
Ext: 


e Videos must not be larger than 30 MB to allow a swift 
technical response for viewing the video. Authors should 
consider whether a short video uniquely captures what 
text alone cannot capture for the understanding of a 
process or behavior under examination in the article. 
Supply an online link to the location of the video. 


Copyright law does not apply to Fishery Bulletin, which 
falls within the public domain. However, if an author 
reproduces any part of an article from Fishery Bulletin, 
reference to source is considered correct form (e.g., Source: 
Fish. Bull. 117:105). | 


Failure to follow these guidelines 
and failure to correspond with editors 
in a timely manner will delay 
publication of a manuscript. 


Submission of manuscript 


Submit a manuscript online at the ScholarOne Manu- 
scripts website for Fishery Bulletin. Commerce Depart- 
ment authors must provide proof of internal clearance 
of their manuscripts with either a completed and signed 
NOAA Form 25-700 or a copy of the clearance email from 
the Research Publication Tracking System. For further 
details on electronic submission, please contact the asso- 
ciate editor, Cara Mayo, at 


cara.mayo@noaa.gov. 


When requested, the text and tables should be submitted 
in MS Word format. Each figure should be sent as a sep- 
arate PDF, TIFF, or EPS file. Send a copy of a figure in 
the original software if conversion to any of these formats 
yields a degraded version of the figure. 


Questions? If you have questions regarding these guide- 
lines, please contact the managing editor, Kathryn 
Dennis, at 


kathryn.dennis@noaa.gov. 


Questions regarding manuscripts under review should be 
addressed to Associate Editor Cara Mayo. 
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